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The response and stability of a rotor system are critically dependent on bearing 
characteristics. Problems involving nonlinear characteristics of bearings - both, rolling 
element and fluid film, arise quite frequently in rotors (Ragulskis et al, 1974; Harris, 
1984; Dimentberg, 1988; Bendat, 1990; Stolarski, 1990; Choy et al, 1992; Childs, 1993; 
Soong and Grigoriu, 1993; Lin and Cai, 1995; Zhou and Hashimoto, 1995). It has 
become evident in recent years that an important class of rotor bearing phenomenon 
cannot be studied without adequately accounting for the nonlinear forces produced by the 
bearings. Bearing nonlinearity assumes a greater role for high speed and low weight rotor 
applications, where the vibration amplitudes tend to be relatively large. In comparison to 
the amount of work done in recent years, in the area of characterisation of linear elastic 
parameters of bearings, research in identification and estimation of the nonlinear 
parameters has been relatively scarce. This is true for fluid film bearings as well as 
rolling element bearings, which, despite their mechanical simplicity, are known to display 
highly nonlinear behaviour and present some very complex rotor problems. 

Parameter estimation is an inverse problem of determination of the elements of the 
system from measures of the forcing function induced upon it and the resultant response. 
Inverse problems in nonlinear analysis require techniques with rigorous theoretical base, 
which can provide valid routes to parameter estimation. The structure of Volterra series 
(Volterra; 1889, 1930, 1959) and Wiener series (Wiener, 1958), which model the 
relationship between the system response and the input in terms of a series of first and 
higher order convolution integrals, provide analytical platforms which can be utilised for 
parameter estimation. 
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The present study attempts to employ the kernel form of response representation of 
nonlinear systems, developed by Volterra/ and Wiener, for estimation of nonlinear 

li 

stiffness parameters of rotor-bearing systems. These theories have been in engineering 
literature for quite some time and have been used by researchers, mainly in electrical 
control systems area, for kernel identification and synthesis. However, Volterra and 
Wiener theories have been scarcely used in mechanical engineering applications. An 
engineering approach is adopted in the present work, whereby these theories have been 
utilised beyond first and higher order kernel identification, to the problem of estimation 
of both linear and nonlinear stiffness parameters of rotor-bearing systems. 

The Volterra series, recognised as a powerful tool for nonlinear analysis of systems 
(Billings, 1980; Rugh, 1981; Korenburg and Hunter, 1990), employs multidimensional 
kernels, which upon convolution with the applied excitation, express the response in the 
form of a power series. The kernels of the system are understood as multi-dimensional 
unit impulse response functions. Two basic difficulties associated with practical 
application of Volterra series are the convergence of the series and measurement of 
Volterra kernels of the given system. Their measurement is possible only if the 
contributions of each of the system’s Volterra operators can be separated from the total 
respon.se. These difficulties are circumvented in the Wiener theory. Wiener functionals 
are formed by a set of orthogonal functionals from the Volterra functionals for an input, 
which is a white Gaussian time function. The Wiener kernels can be determined through 
cross-correlation techniques, which can be carried out in time domain or the frequency 
domain. Reference can be made to the works of Lee and Schetzen (1965), Harris and 
Lapidus (1967), French and Butz (1973), Frentch and Butz (1974), Bedrosian and Rice 
(1975), (Yum (1975), Clifford and Tomlinson (1989), Bendat (1990, 1998), Bendat and 
Piersol (1993), for kernel and higher order FRF measurements. The frequency domain 
treatment of Volterra and Wiener series offers an easier computation and more intuitive 
interpretation, in comparison to time domain treatment. 



A frequency domain methodology has been developed, in the present study, for nonlinear 
rotor-bearing systems, whereby the Wiener kernels are estimated from measurements of 
the system response and applied white (broad band) excitation. Excitation is applied at 
the bearing stations and resultant vibration responses are picked up at the bearings and 
other required stations along the rotor. Volterra form of response representation is 
employed to synthesise mathematical expressions of the kernels in terms of linear and 
nonlinear parameters of the system. The relationships between Volterra and Wiener 
kernels of identical order, based on orthogonality conditions, are then employed to obtain 
estimates of the system parameters. 

fhe parameter estimation procedure is developed and illustrated for the following rotor 
configurations ~ 

(i) rigid rotor in bearings without cross-coupling 

(ii) rigid rotor in bearings with cross-coupling 

(iii) flexible rotors 

fhe case oi' a rigid rotor supported in bearings with negligible cross-coupling effects, as 
in the case of ball bearings, has been modeled as a single-degree-of-ffeedom system. The 
analytical development, though general in form, is done through an assumption of a cubic 
tvfx- of nonlinearity in the system. An engineering approach is suggested through a third 
order response representation of Volterra and Wiener series. Damping is taken as linear. 
Using frequency domain analysis, the first to third order kernels are extracted from 
measurements of the applied force and response. A third order kernel factor is 
svnthesised from the first order kernel transform and is processed along wth the 
mciisured third order kernel transform for nonlinear stiffness estimation. The procedure is 
illustrated through numerical simulation. 


The Volterra kernels are obtained from the Wiener kernels for a white Gaussian 
excitation. However, in practice, a broad band excitation can be provided and finite 
samples of data are used to obtain an estimate of a kernel. Though, the estimates are 
expected to differ from the true values, by amounts dependent on the variance of the 
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means, formed from the finite samples, errors can be made to diminish by averaging over 
more data. The analysis is developed in terms of nondimensional parameters. A 
numerical simulation of the response and subsequent parameter estimation is carried out 
and errors involved are illustrated for various sets of nondimensional parameters. 

A rigid rotor supported in bearings with cross-coupling effects, as in the case of fluid film 
bearings, is modeled as a two-degree-of-freedom system. This case is more involved, 
since in addition to the set of direct kernels, cross-kernels also need to be defined in this 
case, which can be convolved with the applied excitation for response representation in 
terms of Volterra and Wiener series. Starting with the equations of motion, Laplace 
transforms are employed to derive expressions for the firet and higher order direct and 
cross-kernels. Again, analytical modeling is carried out for a cubic type of nonlinearity. 
Damping is taken as linear and without cross-coupling, in order to keep the model simple. 
The first and third order kernels are extracted from the force and response measurements. 
Third order direct and cross-kernel factors are synthesised from the measured first order 
kernel transforms. These are subsequently processed with the measured third order kernel 
transforms for parameter estimation. The procedure is illustrated through numerical 
sinuiiati<m, for various sets of nondimensional parameters and the effect of 
approximations and enors involved are discussed. 

fhc procedure has been further developed to incorporate the influence of shaft flexibility. 
Influence coefficient ct>ncept.s are employed to obtain shaft stiffness for the analysis. In 
order to reduce algebraic complexities, the shaft stiffness is assumed to be linear and only 
the hearing nonlinearity has been included. The analysis has been carried out in matrix 
form, to the extent possible, in order to retain the generality of approach. The parameter 
estimation procedure is numerically illustrated for a set of nondimensional parameters. 

'I'he application aspects of the procedure are studied through experimental investigations 
on a laboratory rotor rig. The experimental work has been restricted to a simple 
configuration, a rotor supported in ball bearings, which is treated as a single-degree-of- 
freedom system. The excitation mechanism, instrumentation and the results obtained are 



discussed. The estimates of the parameters are compared with those obtained from a 
previous study based on Markov process approach. The results are also validated through 
an independent check, which employs Hertzian contact theory for obtaining stiffness 
estimates of bearings in isolation of the shaft. 

To summarize, linear and non-linear bearing stiffness estimation procedures have been 
developed for cases of rigid rotors in bearings with and without cross-coupling and 
flexible rotors. The procedures are based on the theoretical foundations of the kernel form 
of nonlinear response representation of the Volterra and Wiener series. The procedures, 
though analytically developed for cubic type nonlinearity, have a general form and can be 
employed for estimation of a general polynomial form of nonlinearity. The procedures 
arc also applicable to multi-degree freedom systems. The algorithms are illustrated 
through numerical simulation and experiments on a laboratory rotor-bearing test rig. The 
results and errors involved due to various engineering approximations are discussed. 
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CHAPTER 1 


INTRODUCTION 

Dynamic characterization of bearings constitutes a major area in rotordynamic studies. 
The response and stability of a rotor system are critically dependent on bearing 
characteristics. It has become evident in recent years, that an important class of rotor 
bearing phenomenon cannot be studied without adequately accounting for the nonlinear 
forces produced by the bearings. Bearing nonlinearity assumes a greater role for high 
speed and low weight rotor applications, where vibration amplitudes tend to be relatively 
large. 

While considerable amount of work has been done in recent years in estimation of linear 
elastic parameters of bearings, research in identification and estimation of the nonlinear 
parameters has been relatively scarce. This is true for fluid film bearings as well as rolling 
element bearings, which despite their mechanical simplicity, are known to display highly 
nonlinear behaviour and present some very complex rotor problems. The approximations, 
thus involved in the dynamic characterisation of bearings are responsible for some of the 
unreliability in prediction of the response and stability of a rotor system. Parameter 
estimation is an inverse problem of determination of the elements of the system from 
measures of the forcing functions induced upon it and the resultant response. Estimation 
of linear parameters of rotor-bearing systems has been mostly carried out within the 
framework of classical dynamics, where closed form solutions are generally available to 
model the system response. Work in nonlinear parameter estimation has been 
handicapped, due to absence of such general models to exactly represent the system 
response. 

The present study attempts to employ the kernel form of response representation of 
nonlinear systems, developed by Volterra euid Wiener, for estimation of nonlinear 
stiffness parameters of rotor-bearing systems. These theories have been in engineering 



literature for quite some time and used by researchers, mainly in electrical control 
systems area, for kernel identification and synthesis. However, they have been scarcely 
used in mechanical engineering applications. An engineering approach is adopted in the 
present work, whereby these theories have been utilised beyond first and higher order 
kernel identification to the problem of estimation of both linear and nonlinear stiffiiess 
parameters of rotor-bearing systems. 

The structure of Volterra and Wiener series offers a rigorous theoretical platform for 
analysis of the response of a nonlinear system for parameter estimation. The Volterra 
series models the relationship between the system response and input by a series of first 
and higher order convolution integrals. The convolution is carried out between the 
applied excitation and the kernels of the system. The kernels of the system are understood 
as multi-dimensional imit impulse response functions. The response can be represented as 
a series of linear, bilinear, trilinear terms and so on. Two basic difficulties associated with 
the practical application of the Volterra series are the convergence of the series and the 
measurement of the Volterra kernels of the given system. Their measurement is possible 
only if the contributions of each of the system’s Volterra operators can be separated from 
the total response. These difficulties are circumvented in the Wiener theory. Wiener 
functionals are formed by a set of orthogonal fimctionals from the Volterra functionals, 
for an input which is a white Gaussian time function. The Wiener kernels can be 
determined through cross-correlation techniques. This determination can be carried out in 
the time domain or frequency domain. 


A general procedure has been developed in the present study, whereby the Wiener kernels 
are estimated from measurements of the system response and applied white (broad band) 
excitation. Volterra form of response representation is employed to synthesise 
mathematical expressions of the kernels in terms of the linear and nonlinear parameters of 
the system. Relationships between Volterra and Wiener kernels of identical order, based 
on orthogonality conditions, are then employed to obtain estimates of the system 
parameters. A frequency domain approach is adopted for analysis. 



The problem is initially attempted for a rigid rotor supported in bearings with negligible 
cross-coupling effects, as in the case of ball bearings. Such a rotor-bearing configuration 
is modeled as a single-degree-of-ffeedom system. The analytical development, though 
general in form, is done through an assumption of a cubic type of stiffness nonlinearity in 
the system. An engineering approach is suggested using a third order response 
representation of Volterra and Wiener series. Damping is taken as linear. Using frequency 
domain analysis, the first and third order kernels are extracted firom measurements of the 
applied force and response. A third order kernel factor is synthesised from the first order 
measured kernel transforms and is processed along with the measured third order kernel 
transform for nonlinear stiffness estimation. The procedure is illustrated through 
numerical simulation. Volterra kernels are developed from the Wiener kernels for a white 
Gaussian excitation. However, in practice, a broad band excitation can be provided and 
finite samples of data can be used to obtain an estimate of a kernel. Though, the estimates 
are expected to differ from the true values, by amounts dependent on the variance of 
means formed from the finite samples, errors can be made to diminish by averaging over 
more number of samples. The analysis is done in terms of nondimensional parameters 
and numerical simulation of the response is carried out. Subsequent parameter estimation 
and the errors involved are illustrated for various sets of nondimensional parameters. 


The procedure is developed further for treatment of a rigid rotor supported in bearings 
with cross-coupling effects, as in the case of fluid film bearings. This case is more 
involved and is modeled as a two-degree-of-fireedom system. In addition to the set of 
direct kernels, such a system needs cross-kernels to be defined, which can be convolved 
with the applied excitation for response representation in terms of Volterra and Wiener 
series. Starting with the equations of motion. Laplace transforms are employed to derive 
expressions for the first and higher order direct and cross-kernel transforms. Analytical 
modeling is carried out for a cubic type of nonlinearity. Damping is taken as linear and 
without cross-coupling, for simplicity. The first and third order kernels are extracted from 
force and response measurements. Third order direct and cross-kernel factors are 
synthesised from the measured first order kernel transforms. These are subsequently 



processed with the measured third order kernel transforms for parameter estimation. The 
procedure is illustrated through numerical simulation, for various sets of nondimensional 
terms and the approximations and errors involved are discussed. 

The case of parameter estimation in flexible rotors is considered next. In the cases 
mentioned above, the rotor shaft is taken as rigid. The treatment has been further 
extended to incorporate the influence of shaft flexibility. Influence coefficient concepts 
are employed to incorporate shaft flexibility. In order to reduce algebraic complexities, 
shaft stiffness is assumed to be linear and only the bearing nonlinearity has been 
incorporated. Analysis has been carried out in matrix form to the extent possible, in order 
to retain the generality of the approach. Parameter estimation procedure is developed 
along the route described in the previous cases. Numerical illustration is carried out for a 
set of nondimensional parameters. 

Experimental investigations have been carried out on a laboratory rotor rig, to further 
study the application aspects of the procedures developed. The experimental work has 
been restricted to a simple configuration of a rotor supported in ball bearings, which is 
treated as a single-degree-of-freedom system. The excitation mechanism, instrumentation 
and results obtained are discussed. Estimates of the parameters are compared v^th those 
obtained from a previous study based on the Markov process approach. The results are 
also validated through an independent check, which employs Hertzian contact theory for 
obtaining stiffness estimates of bearings in isolation of the shaft. 



CHAPTER 2 


LITERATURE REVIEW 

Analysis of the response of non-linear systems subjected to random excitation have 
attracted considerable attention in the recent past, mainly due to growing awareness of 
the significance of the random nature of forces produced by a number of physical 
phenomena like wind, sea waves, earthquakes etc., which must be considered in limit 
design of aircrafts, buildings, off-shore structures, ships etc. The present study, deals with 
the inverse problem of the estimation of system parameters from a knowledge of the 
applied force and the resultant system response. Inverse problems in nonlinear analysis 
require techniques with rigorous theoretical base, which can provide valid routes to 
parameter estimation. Volterra series (Volterra; 1889, 1930, 1959), provides a basis for 
these requirements. It has, by now, been well established as a powerful tool in the 
analysis of nonlinear systems (Billings, 1980; Rugh, 1981; Korenburg and Hunter, 1990). 
The theory of Volterra series involves modeling the relationship between system 
response and input force in terms of a series of first and higher order convolution 
integrals. Wiener series (Wiener, 1958) circumvents the problems associated with 
Volterra series, i.e, convergence of the series and measurement of individual kernels. The 
Wiener series is formed by an orthogonal set of functionals from the Volterra series for a 
white Gaussian input. Such white noise analysis is considered to be an effective method 
for gaining a maximum of information with a minimum number of assumptions about the 
system. 

Volterra and Wiener approaches are adopted, in the present work, for linear and nonlinear 
parameter estimation of rotor-bearing systems. Problems involving nonlinear vibrations 
arise quite frequently in rotordynamics, especially those caused through random 
vibrations of bearings - both, rolling element and fluid film, excited by random loads 
(Ragulskis et ah, 1974; Harris, 1984; Dimentberg, 1988; Bendat, 1990; Stolarski, 1990; 
Choy et al., 1992; Childs, 1993; Soong and Grigoriu, 1993; Lin and Cai, 1995; Zhou and 
Hashimoto, 1995). In some cases, deterministic models prove to be inadequate or at least. 



extremely complex and the phenomenon can be adequately described only within the 
framework of statistical models. Statistical dynamics, concerned with the study of various 
random phenomena in dynamic systems, enriches the classical basic theory of oscillations 
and extends the possibilities for its applications. 

A review, of the various approaches to analysis of nonlinear random vibrations, followed 
by a discussion on the ongoing research on dynamic characterization of rotor bearings, is 
presented here. 

2.1 VOLTERRA AND WIENER SERIES 

Volteixa series has been in engineering literature for quite some time now. However, 
there have been few attempts towards its application to mechanical systems. Volterra 
series represents the response of a system in a functional form, through a series of first 
and higher order convolution integrals, involving explicit operations on the input to the 
system. The alternative form of Volterra series is the operator form. In the functional 
form, the emphasis is on the output at a specific time instant, for a given input function, 
while in operator form the focus of attention is the complete output function. The 
response is represented as a series of linear, bilinear, trilinear terms and so on. A linear 
operator is one for which the response to a linear combination of signals is same as the 
linear combination of the response to each individual signal. The functional form of the 
linear operator is a first order integral involving convolution between the first order 
kernel and input. The first order kernel is a unit impulse response function. The second 
order operator is one for which the response to a linear combination of signals is a 
bilinear operation on individual input signals. The functional representation of the 
bilinear operator is a two dimensional integral involving convolution of a second order 
kernel with the input. The second order kernel can be interpreted as a two dimensional 
impulse response. The trilinear and higher operators are similar extensions to higher 
dimensions. There have been many articles devoted to theoretical issues such as existence 
of Volterra series (Leon. 1978; Lesiak. 1978; Sandberg, 1982), computation of Volterra 
kernels of special systems, composition, feedback configurations, nonlinear circuits etc. 
(Brilliant, 1958; Barret, 1963; Flake, 1963; Lubbock, 1969; Chua, 1979, 1979a; Boyd et 



al., 1983) and the formal frame work for Volterra series (Halme, 1971; DeFigueiredo, 
1980, Sandberg, 1982). Volterra series can be treated in both — the time domain and the 
frequency domain. The frequency domain treatment provides a logical and appropriate 
way of extending the linear system theory to higher order Frequency Response Functions 
(FRFs). The higher order FRFs are calculated using techniques similar to those applied to 
linear systems. The obvious advantage of the frequency domain treatment over time 
domain treatment is that frequency domain kernels are easier to compute and have an 
easier interpretation. Two basic difficulties are associated with practical application of the 
Volterra series. The first concerns measurement of the Volterra kernels of the given 
system. Measurement is possible only if contributions of each of the system’s Volterra 
operators can be separated from the total response. No exact method of isolating an 
individual Volterra operator exists, for a system not predefined to be restricted to a finite 
order. The second difficulty is that Volterra series representation of a physical system 
may converge for only a limited range of the system input amplitude. Volterra series is a 
power series with memory, similar to Taylor series, and suffers from the same limitations 
as the Taylor series. These difficulties vvath the Volterra series are circumvented in the 
Wiener theory. 

Wiener functionals involve forming a set of orthogonal functionals from the Volterra 
functionals. These functionals are orthogonal when the input is a white Gaussian time 
function. The convergence of an orthogonal series is a convergence in the mean and 
Wiener functionals describe a larger class of nonlinear systems (Schetzen, 1980). The 
Wiener kernels render themselves to individual determination through means of cross- 
correlation. This determination can be carried out in the time domain or the frequency 
domain. White Gaussian excitation provides a means of exciting all possible nonlinear 
behaviour of a system, within a certain amplitude limit. However, in practice, a broad 
band excitation can be provided and finite samples of data can be used to obtain an 
estimate of a kernel. Though, the estimates differ from the expected values, by amounts 
dependent on the variance of the means formed from finite samples (Fakhouri et al., 
1981), errors can be made to diminish by averaging over more data. 



Considerable amount of research has been carried out on various aspects of Volterra zind 
Wiener series analysis of nonlinear systems. Some of the salient works are discussed 
below. 

Early studies: 

Fretchet (1910) developed a functional series, called as Fretchet-Volterra series, for 
representation of a continuous nonlinear system. Wiener (1958) considered the 
identification of nonlinear systems by representing each functional term by a Fourier- 
Hermite series in conjunction with Laguerre functions. Brilliant (1958) has used the 
frequency domain approach in the analysis of various systems composed of cascade, 
feedback, inversion, addition, and subtraction of the elements. A number of studies 
stemmed from this work, for synthesis of a system to be equated to an unknown system. 
Amorocho and Orlob (1961) used the basis functional approach in finding the kernels of 
rainfall-runoff process. Jacoby (1966) introduced a decomposition model for hydrologic 
systems, A number of calculation methods were developed by later investigators using 
direct solution, basis functions, time domain, and frequency domain approaches, a review 
of which has been made by Hung and Stark (1977). 

Time domain approach: 

Schetzen (1965) presented a method for measuring the Volterra kernels of a finite order 
nonlinear system. The kernels are obtained individually as a multi-dimensional impulse 
response. Multi-dimensional kernel transforms can also, be obtained by the method 
described. As an extension, a technique of obtaining Volterra kernels from a multi- 
dimensional step response was presented. This technique is useful for nonlinear systems 
which can be considered to be of a given finite order, for only a limited range of input 
amplitudes. In a subsequent paper (1965(a)), he has presented a method for synthesis of 
Volterra kernels. Later, Lee and Schetzen (1965) developed a practical and relatively 
simple method for measuring Wiener kernels of a nonlinear system. The method is based 
upon cross correlation techniques and does not involve orthogonal expansion of basis 
functions as in the Wiener method. 



Katzenelson and Gould (1962) developed an iterative method for obtaining kernels that 
can represent a system. They posed the system equivalently in terms of the solution 
needed for an optimum filter. Korenberg (1973) considered identification of differential 
systems having a Volterra series expansion, using a slowly exponentially decaying sum 
of sinusoids as input. The terms of the differential expansion are determined orthogonally 
using linear regression and simple averaging procedures. A method of directly identifying 
the Volterra kernels using an exponentially decaying function multiplied by a bounded 
zero-mean independent process, was also developed by Korenberg (1973). Identification 
is orthogonal, and the kernels are obtained by averaging procedure. Marmarelis and Naka 
(1973, a-c) and Watanabe (1975) discussed two input approach to extract kernel 
information. Krausz (1975) developed an identification method based on a random 
impulse train (Poisson process) rather than white Gaussian input. This method appears to 
be especially applicable to study of neuronal interactions, where the input is stimulation 
of a nerve axon. Klein and Yasui (1976) proposed the concept of ‘dual-space’ basis 
functions, which are then expanded in terms of Volterra kernels and further, related to 
Wiener kernels. Fakhouri (1978) developed an algorithm for identification of discrete 
Volterra kernels in terms of multidimensional z-transforms using higher order correlation 
functions and coloured Gaussian inputs. Fakhouri (1980) and Billings and Fakhouri 
(1980) analysed the nonlinear feedback system and developed an identification algorithm 
for open and closed loop nonlinear systems based upon pseudorandom excitation, which 
pro\ ides estimates of the individual component subsystems from measurement of the 
input and noise corrupted output. 

Freouenev domain approach: 

After the initial works of Brilliant (1958) and George (1959), a number of studies were 
carried out employing the firequency domain approach. Tick (1961) developed a measure 
of quadratic coherence, which is used to indicate the amount of quadratic effect on the 
output due to input. Barret (1963) introduced the use of functional expansions for 
calculation of output when the input is sinusoidal or random. Korenberg (1973) 
de\ eloped a general method for explicit solution of the transfer functions of a random- 
stimulus-input cascade system containing a chain of alternating linear and static nonlinear 



systems. Brillinger (1970) developed a procedure for an asymptotically unbiased and 
asymptotically normal estimation of the jth degree frequency transfer function for 
identification of a general unknown system. 

French and Butz (1973) developed a general frequency-domain method for calculation of 
higher degree Wiener kernels, analogous to the cross correlation method of Lee and 
Schetzen (1965) for the time-domain case. They have used complex exponential 
functions as a set of orthogonal functions for expanding the kernels. In another work 
Frentch and Butz (1974) have developed an algorithm based on expansion of Wiener 
kernels in terms of Walsh functions. The nonlinear system is described in terms of a set 
of kernels, which contain dyadic convolution operation, and identification is performed 
using the fast Walsh-Fourier transform. Crum (1975) presented a simple procedure for 
the simultaneous reduction and expansion of multidimensional Laplace-transform 
kernels. 

Bedrosian and Rice (1971, 1975) illustrated frequency domain methods for solution of 
some simple nonlinear combination of elements by means of harmonic probing, which 
involves simple harmonic functions as inputs and employs direct-expansion method. 
Gifford and Tomlinson (1989) illustrated the technique of calculation of higher order 
FRF’s, for systems with more than one degree of freedom. The technique involving 
correlation measurements using random excitation is illustrated by analysis of a nonlinear 
beam. They also introduced higher order FRF models of a discrete structural system and 
showed how this model can be used as a basis for curve fitting and extracting a 
parametric model of multi-degree of freedom nonlinear system. 

Nam et al. (1990) have presented a frequency domain approach for digitally estimating 
the system parameters, based on input/output measurements, of a nonlinear system which 
can be approximately characterised in terms of the first three transfer functions of 
Volterra series. As an extension of the work of previous investigators (Tick, 1961; Koh, 
1985; Nikias. 1987; Powers, 1987; Kim, 1988), they presented a unified approach for 
handling cubic systems with stationary random inputs, by providing a quantitative 



measure of power transfer between the system input and output. Odiari and Ewins (1992) 
have presented a procedure for identifying the vibration parameters of nonlinear vibratory 
systems, with particular emphasis on rotor-stator systems. They have derived the 
equations for calculating the Wiener and Volterra kernels in the frequency domain. 

Modal techniques (Ewins, 1984) form the basis of most of structural analysis of general 
degree of freedom systems. Bendat and Piersol (1982, 1986), have derived equations 
based on bispectral analysis techniques for finite memory square law systems. Later, 
Bendat (1990, 1998); Bendat and Piersol (1993) derived a range of formulae for finding 
the first order, second order and third order FRP’s through spectral analysis. Application 
of bispectral analysis methods to quadratic type system has also been reported by 
Yamamouchi (1974), for describing ocean wave characteristics and by Choi et al. (1985) 
and Vandiver and Jong (1986) for flow induced problems. Rice and Fit 2 patric (1988) 
have given a generalized method for identification of a large range of nonlinear systems 
subjected to random excitations using spectral method in frequency domain. Rice and 
Fritzpatic (1991) have also discussed a situation where the modal technique becomes 
unsatisfactory and an inversion approach has been adopted. Recently, Worden et al. 
(1997) have extended the conventional harmonic probing algorithm of Bedrosian and 
Rice (1975) to deal with multi-input multi-output form of Volterra series. 

System Characterisation: 

There have been two approaches to the system characterization problems. The first 
approach is the differential equation method (Eykhoff et al. 1966; Eykhoff, 1974), where 
the topology of a system is assumed to be known, so that a set of differential equations 
can be used to represent the system. Identification, therefore, reduces to determination of 
various parameters in the equations. The second approach (Hung, et al., 1977) is the 
integral equation method, where little or no a-priori assumptions are made about the 
topology of the system. Instead, the form of integral equations is fixed, and the 
identification method reduces to determination of values within the integrals, called 
kernels. When the internal structure of a system is known, the differential equation 
method is usually preferred to the integral method. Only a few parameters need to be 



identified in the set of differential equations, since often only a finite number of kernels 
need to be estimated. Another advantage of differential over integral equation approach is 
that the former is usually easier to interpret. On the other hand the power of the integral 
equation approach lies in its canonical nature (Wiener, 1958). The method can be applied 
to a wide class of systems, which are essentially time-invariant and have finite memory. 
Ho and Stark (1973) and Watanabe and Stark (1975) carried out interpretation of their 
estimated kernels. They have also emphasised the examination of the response dynamics 
of the kernel models to transient inputs, such as pulses, steps, and ramps etc. 

As mentioned in the introductory part of this chapter, statistical methods have also served 
as powerful tools for analysis of stochastic dynamics of nonlinear systems. A brief 
discussion of such methods is presented here. 

2.2 STATISTICAL METHODS 

Statistical Dynamics, concerned with study of various random phenomena in dynamic 
systems, extends the possibilities for application of the classical basic theory of 
oscillations to situations where deterministic models are inadequate. Recognition of white 
noise excitation as an efficient tool for extracting maximum information about the 
system, with least number of assumptions has served to underline the significance of such 
procedures for parameter estimation. 

Markov process approach: 

'Hie significance of the Markov process approach in parameter estimation lies in the fact 
that exact solutions of a limited class of vibratory systems can be obtained. A process, 
whose present probability distribution depends on only one previous time instant, is 
called a Markov process. The structure of a Markov process is completely determined, 
for all future times, by the distribution at some initial time and by a transition probability 
density function, which satisfies a linear partial differential equation known as the 
Fokker-Planck-Kolmogorov (Fokker, 1914; Planck, 1917; Kolmogorov, 1931) equation. 
The Fokker-Planck equation was derived and exact stationary response was obtained for 
certain cases of two^degree-of-freedom non-linear dynamic systems, by Ariaiamam 



(1960). The theory was generalized later by Caughey (1963) for multi-degree of freedom 
cases. The drift and diffusion coefficients in the FPK equation can be derived from the 
nonlinear equations of motion of the d)mamic system (Crandall, 1966). A general, closed 
form solution to FPK equation is yet to be found. However, the first order probability 
distribution for the stationary response distribution is readily obtained for a limited class 
of vibratory systems with nonlinear restoring forces and special forms of nonlinear 
damping (Fuller, 1969; Caughey, 1971; Caughey and Ma, 1983; Pradlwarter et al., 1991). 
Exact stationary solution in terms of probability density function, for a class of non-linear 
systems driven by a non-normal-delta-correlated process has been obtained by 
Vasta(1995). 

For the response to be approximated by a Markov process, it is necessary that the 
excitation be approximated by ideal white noise. This restriction can be, in principle, 
removed at the price of increasing the complexity of the system, by introducing linear 
filters between the ideal white-noise excitations and the system. It can be shown by using 
stocha.stic averaging principles that, under certain conditions, the response of a nonlinear 
dynamic system to non-white excitation can be approximated by a Markov process. 
Relevant works in this field are by Stratonovitch (1967), Khasminskii (1966), 
Papanicolaou and Kohler (1974), Dimentberg (1980), Roberts (1983), Spanos (1983) and 
Zhu (1983). These works involve extensive applications to nonlinear random vibration of 
mechanical and structural systems. 

Perturbation methods: 

Perturbation methods are generally employed, when the amount of nonlinearity in a 
system is controlled by a small scaling parameter. The solution is sought in terms of a 
power series in a small scaling parameter and successive terms are evaluated as linear 
responses to nonlinear functions of the preceding ternis. This classical approach for 
deterministic nonlinear problems (Stoker, 1950) was extended to random vibration 
problems by Crandall (1963). In practice, calculations are seldom carried beyond the first 
perturbation. Functional series methods offer an alternative approach to developing an 
expansion, based on the linear solution. An example of application of this method is the 



work of Orabi and Ahmadi (1987). They used a Wiener-Heimite expansion and presented 
a formal procedure for deriving the deterministic equations governing kernel functions, 
arising in the expansion (Roy and Spanos, 1990). 

Method of moments: 

A set of differential equations for various statistical moments, or related quantities known 
as cumulants (or semi-invariants) and quasi-moments (Stratonovitch, 1967) of the 
response, as function of time, can be obtained by multiplying the FPK equation by 
suitable functions and integrating over the probability space. Equivalent sets of equations 
can be derived directly from the dynamic equations of motions or the equivalent Ito 
equations. Approximate solutions have been proposed, based on ad-hoc closure 
assumption (Bolotin, 1979; Ibrahim and Roberts, 1978; Roberts, 1981, Crandall, 1985). 

Method of equivalent nonlinear equations: 

An alternative generalization of statistical linearization has been proposed by Caughey 
(1986). The idea is to replace the original set of nonlinear differential equations by an 
equivalent nonlinear set, where the latter belong to a class of problems, which can be 
soh ed exactly. This class is, at present, very limited, and thus the range of applicability 
of the technique is correspondingly restricted. Results have been obtained for oscillators 
with nonlinearity in damping and stiffness (Caughey, 1986; Cai and Lin, 1988; Zhu and 
Yu. 1989), Equivalent statistical quadratization methods have been applied to non-linear 
multi -degree-of-freedom systems subjected to random excitation by Spanos and Donley 
(1991, 1992). It has been demonstrated that the method is very effective as a means of 
predicting the probability distribution of response, with reasonable accuracy. 

Method of computer simulation: 

Numerical simulation or the Monte Carlo method (Shinozuka, 1972; Bolotin, 1979; 
Spanos, 1981; Rubinstein, 1981; Spanos and Mignolet, 1989) consists of generating a 
large number of sample excitations, computing the corresponding response samples and 
processing them to obtain the desired response statistics. The backbone, of any digital 



simulation study, is an algorithm, which provides a set of pseudo-random numbers, 
belonging to a population with specified probability density function. Proper processing 
of this set of numbers can yield values of sample functions of random process excitations, 
with pre-selected frequency content and temporal variation of intensity, at successive 
discrete equi-spaced times. Upon generating a single sample of excitation, commonly 
available subroutines for numerical integration of differential equations can be employed 
to obtain the system response. Another sample of excitation can then be generated and 
the computed values of system response can be used to update its statistics. The 
procedure, in principle, is very general and applicable to stationary or non-stationary 
response of systems of any degree of complication. 

2.3 ROTOR-BEARING CHARACTERISATION 

Parameter estimation of the elastic parameters of bearings involves establishing a 
relationship between the incident load on the bearing and its resultant deformation. Some 
of the relevant studies available in literature in the area of dynamic characterisation of 
rolling element and fluid-film bearings are discussed in this section. 

Rolling element bearings: 

Hertzian contact theory (Hertz, 1896) for the solution of local stress and deformation of 
two elastic bodies apparently contacting at a single point provides a platform for 
determination of the elastic parameters of rolling element bearings. Hertz's analysis is 
applied to surface stresses caused by a concentrated force, applied perpendicular to the 
surface. In determination of the contact deformation versus load curve, concentrated load 
applied normal to the surface alone, is considered, for most rolling element bearing 
applications. Methods of calculation of surface and subsurface stresses under a 
combination of normal and tangential (traction) stresses are complex, (Zwirlein and 
Schlicht, 1980). Owing to infinitesimally small irregularities in the basic surface 
geometries of the rolling contact bodies, neither uniform normal stress field nor a tmifoim 
shear field are likely to occur in practice (Sayler et al., 1981; Kalker, 1982). Rigorous 
mathematical/numerical methods have been developed to calculate the distribution and 



magnitude of surface stresses in any line contact situation, that is. including the effects of 
crowning of rollers, raceways, and combinations thereof (Kunert, 1961; Reusner, 1977). 
Additionally, finite element methods (FEM) have been employed (Fredriksson, 1980) to 
perform the same analysis. 

It is possible to determine how bearing load is distributed among the balls or rollers, after 
having determined how each ball or roller in a bearing carries load. To do this, it is 
necessary to develop load-deflection relationships for rolling elements contacting 
raceways. Most rolling bearing applications involve steady-state rotation of either the 
inner or outer raceways or both. Rolling element centrifugal forces, gyroscopic moments 
and frictional forces and moments do not significantly influence this load distribution in 
most applications. Theoretical models (Palmgren, 1959; Ragulskis et al., 1974; Harris, 
1984; Eschamann et al. 1985; Stolarski, 1990) are available for estimation of bearing 
stiffnesses under static loading conditions. Recently, Chen and Lee (1997) presented a 
method for estimating linearised coefficients of rolling element bearings, which uses 
relations of unbalance responses and known system parameters in the construction of an 
estimator. 

Bearing vibrations caused due to geometric imperfections of contact surfaces were first 
analyzed by Lohman (1953) and Gustavsson (1962). Comprehensive investigations have 
been carried out on the high frequency response of bearings (McFadden and Smith, 1984) 
and its relation to surface irregularities ( Sunnesjo, 1985; McFadden and Smith, 1985; Su 
et al., 1993). Lim and Singh have analyzed the vibration transmission through rolling 
element bearings in series of publications (1990a, 1990b, 1991, 1992, 1994). 


A method for determination of the nonlinear characteristics of bearings using the 
procedure of Krylov-Bogoliubov-Mitropolsky has been suggested by Kononenko and 
Plakhtienko (1970). Honrath (1960) and Elsermans et al. (1975) have experimentally 
examined the stiffness and damping of rolling element bearings. Walford and Stone 
(1980) designed and fabricated a test rig for direct measurement of the relative 
displacement of the shaft and bearing housing for the oscillating force applied to the 



bearing housing, which is used to obtain the stiffness parameters. They found that the 
interfacial play between races and housing and shaft play a significant role in the 
detemiination of bearing stiffness and damping. Kraus et al. (1987) presented a method 
for extraction of rolling element bearing stiffness and damping under operating 
conditions. The method is based on experimental modal analysis combined with a 
mathematical model of the rotor-bearing-support system. Effects of speed, preload and 
free outer race bearings on stiffness and damping have been investigated. Muszynska 
(1990) has developed a perturbation technique for estimation of these parameters. The 
technique involves a controlled input excitation to be given to the bearings. Goodwin 
(1991) reviewed the experimental approaches to rotor support impedance measurement, 
with particular emphasis on fluid film bearing impedance measurement. A general 
procedure for identification of restoring force nonlinearity from system’s response to a 
white-noise excitation has been discussed by Dimentberg and Sokolov (1991). Nonlinear 
stochastic contact vibrations and friction at a Hertzian contact have been studied by Hess 
et at. (1992). Analytical and experimental studies are carried out by them, using the 
Fokker-Planck equation and simulating the vibrations to the contact region, either 
externally by a white Gaussian random normal load or internally by a rough surface 
input. 

Tiwari and Vyas (1995, 1997, 1997a, 1998) in a series of papers described procedures for 
e.stimation of non-linear elastic parameters of bearings based on the analysis of random 
response signals picked up from the bearing caps. The procedures did not require an a- 
priori knowledge of the random excitation force induced on the system. The dynamics of 
the rotor-bearing system was modeled as a Markov process and Fokker-Planck equations 
arc fonnulatcd. The Fokker-Planck equations were solved and the response was 
processed, for the inverse problem of parameter estimation. The procedure was 
developed for cases of rigid rotors; single disc flexible rotors and multi-disc flexible 
rotors. The study was extended further, to include harmonic excitation to the non-lmear 
system along with random excitation and the case of an unbalanced rigid rotor was 
discussed. The algorithms were also verified experimentally, for a laborator>' rotor- 

bearing test rig. 



Fluid-film bearings: 


Vast amount of literature is available on the linear aspects of rotors supported in fluid 
film bearings. Reference can be made to the texts by Fuller (1984), Childs (1993) and 
Rao ( 1 998) for a detailed study. Some of the relevant works are discussed below. 


The idea of representing the dynamic response characteristics of a journal bearing by 
means of stiffness and damping coefficients originates with Stodola (1925) and Hummel 
(1926). Their aim was to improve the calculation of the critical speed of a rotor by 
including the flexibility of the bearing oil film. Concurrently, Newkirk (1924, 1925) 
described the phenomenon of bearing induced instability, which he called oil whip. 
Robertson (1946) analysed the whirling of a journal in a sleeve bearing. The influence of 
fluid film journal bearing on the stability of rotors was reported by Hagg (1946). 

The bearings have been approximated, in the initial stages, as short bearings (Oevirk, 
1952) or a.s infinitely long bearings (Gross, 1962; Booker, 1965). Both types of 
approximation.s have been usefully employed by researchers for bearing analysis. The 
analytic finite length bearing model is given by Warner (1963). Later works include those 
of Kirk and Gunter (1970, 1975, 1975a) who employed the short bearing model for rotor 
dynamic anaiysi.s; Simandri and Hahn (1975) and Tonneson (1975) used the short bearing 
model in the analysis of squeeze film dampers. Vance and Kearton (1975) examined the 
appropriatenos.'; of the long bearing model for long squeeze film dampers with end seals. 

Several iin e.-Jtigator.s have investigated the combined transient rotor bearing dynamic 
problem b\ numerically solving the Reynold's equation for the bearing reaction force, 
while simulianeously integrating the rotor equations of motion. Lund and Stemlicht 
(1962) initially calculated the linear stiffness and damping coefficients through numerical 
differentiation of a finite difference solution to the Reynold's equation. Blok (1965) 
soh'ed the Re>-nold's equation directly via finite differences to obtain mobility 
descriptions for finite length bearings. The numerical approach is also followed in the 
works of Kirk and Gunter, mentioned earlier, for the short bearing model and by Myrick 
and Ry lander (1975) for finite length bearing. Orcutt and Arwas (1967) used a similar 



analysis approach to obtain stiffness and damping coefficients for both laminar and 
turbulent conditions. 

Lund (1968) developed a perturbation solution to the Reynold's equation, which 
eliminates the requirement of numerical differentiation. A direct analytic method for 
deriving a complete set of analytic stiffness and damping coefficients from impedance 
descriptions is developed and demonstrated for the cavitating finite length bearing 
impedances by Childs et al. (1977). The development of the concept of spring and 
damping coefficients for journal bearing is briefly reviewed by Lund (1987). Baheti and 
Kirk (1994, 1994a) used finite element method to solve the nonlinear coupled 
hydrodynamic and thermal equations to obtain the stiffness and damping coefficients of 
seals. 

The direct experimental determination of stiffriess and dam ping coefficients of journal 
bearing was initiated by Glienicke (1966) and Morton (1971). They used the technique of 
excitation of bush sinusoidally in two mutually perpendicular directions, and in each case 
mea.suring the amplitude and phase of the resulting motions. One way of improving the 
accuracy of measurement of stiffness and damping coefficients is to collect experimental 
data over a wide range of forcing frequencies, and to use all of this data to determine the 
coefficient values. In publications by Stanway et al. (1979), by Burrows et al. (1981), and 
by Burrows and Sahinkya (1982) the system is forced at all frequencies within the range 
of interest, simultaneously in horizontal and vertical directions. A method of treating the 
unbalance response data from a test rig to determine bearing stiffness and damping 
coefficients was suggested by Sahinkaya and Burrows (1984). A similar approach 
through the transient response to an impulse or step force was given by Nordman and 
Schollhom (1980). Goodwin (1991) has given a concise account of the experimental 
techniques available for use in the measurement of bearing impedances and discussed 
the potential advantages and disadvantages of each approach. The frequency swept 
rotating input perturbation techniques, for identification of fluid film rotor bearing system 
dynamic characteristics are described by Muszynska and Bentley (1990). The method 
developed by Rouvas and Childs (1993) is designed to account for an unmeasurable input 



to the system, namely the force generated by turbulence and cavitation in fluid film, by 
assuming that the applied excitation and the fluid film induced forces are statistically 
independent. 

Krodkiewski and Ding (1993) developed an on site identification algorithm for multi- 
bearing rotor systems based on monitoring the trajectories of the relative motion of the 
journals with respect to the bearings. Arumugam et al. (1997) have carried out 
experimental investigations to identify static and dynamic characteristics of journal 
bearings under the influence of various cases of twisting misalignment. Goodwin et al. 
(1997) described a combined theoretical and experimental investigation of linear fluid 
film stiffness and damping coefficients. A multi-frequency force signal was used to 
excite the bearing and measurements of the relative movement between the shaft and 
journal have been made. 

Investigations into the nonlinear aspects of fluid-film bearing dynamics have been 
relatively few and recent. Choy et al. (1992) have examined the nonlinear characteristics 
and their effects on the dynamic performance of a hydrodynamic journal bearing at 
various operating conditions. Krodkiewski and Ding (1993a) considered the problem of 
static indetermination in the mathematical model for nonlinear dynamic analysis of multi- 
bearing systems. Recently, Chu et al. (1998) have described a nonlinear dynamic model 
for hydrodynamic bearings which incorporates the nonlinear stiffness and damping 
parameters as functions of static bearing stiffriess. A finite difference approach is used to 
solve for these coefficients. For a detailed account on available experimental data on 
hydrodynamic bearings, reference can be made to the extensive survey by Swanson and 
Kirk (1997), which catalogues about a 100 published experimental works. 



CHAPTER 3 


PARAMETER ESTIMATION IN RIGID-ROTORS 
SUPPORTED IN BEARINGS WITH NO CROSS-COUPLING 

A rigid-rotor supported in bearings with no cross-coupling effects is considered as a first 
exercise in linear and nonlinear parameter estimation. This may be the case for a rotor 
supported in rolling element bearings, where the cross-coupled stiffhess terms are negligible 
in comparison to the direct stiffness terms. Similarly, damping is taken to be comprising of 
the direct terms only. The effect of the rotor shaft flexibility is also not accounted for, in this 
chapter. Rolling element bearings are generally characterized by a cubic nonlinearity. For 
algebraic ease, only stiffness nonlinearity is considered. Damping is taken to be linear. This 
reduces the rotor-bearing set-up to a single-degree-of-freedom system with cubic stiffhess 
nonlinearity. An engineering approach for parameter estimation is developed through a third 
order Volterra kernel representation of the system response. Using frequency domain 
analysis, the first to third order kernels are extracted from measurements of the applied force 
and response. A third order kernel factor is synthesised from the first order kernel transform 
and is processed along with the third order measured kernel transform for estimation of the 
nonlinear parameter. The procedure is illustrated through numerical simulation. The 
assumptions in\olved and the approximations are discussed. The influence of excitation 
force, linear damping parameter and probable measurement noise on the estimates is 
illustrated through nondimensional simulation. 

3.1 Governing Equation And Volterra Series Representation of Response 

For the rotor-bearing system shown in Fig. 3.1. modeled as a single-degree-of-freedom 
system, the governing equation of motion can be written as 


mx + cx + kx + k^ = f{t) 


(3.1) 




where m is the mass, is the unknown non-linear stiffiiess, while c and k are the unknown 
linear damping and stiffness terms. /(/) in the above equation represents the random 
excitation given to the system. 


The equation (3.1) is rewritten in nondimensional form as 

Ti\t) + 2^T]\r) + tj{t) + Xri^ir) = /(r) 

where ( ' ) denotes differentiation with respect to t , and 
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Employing Volterra theory the system response is taken to be of the following kernel form 
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where Hq is a constant and the n* order Volterra kernel 

h„(Ti r„) = 0 for tj{0, y=l,2 n 

The response can be, alternatively expressed, in operator form as 

m = H[7(r)] 
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3.2 Synthesis of Higher Order Volterra Kernel Factors 

The approach suggested by Schetzen (1980) is employed to define the Laplace transform of 
the first-order kernel, Hi (s) , in terms of the linear parameters of the system, namely G>„ and 

The expressions for higher order kernel transforms H 2 is), H^^s) etc. are consequently 
synthesised from the first-order transform, Hy ( 5 ) , and the nonlinear parameter A . 


Replacing the applied force /(r) by cf{t), the system response, fi-om equation (3 .4) is 

Sc"//„[7(r)]= £c";;„(r) 


where, c being a constant and for convenience 

Substituting the new force and the response of equation (3.7) in equation (3.2) gives 

+ 2^|ic”;7„(t)| +|xc";7„(r)| + A|lc";7„(r)| =c7(t) (3.9) 

I'he above power series representation of the governing equation is solved by equating the 
coefficients of like powers of c. 




Equating the coefficients of the first power of c, one obtains 

(3 10) 

Noting from equation (3.4) that q, (r) , the first term in the Volteira representation of > 7 (r) , 
is the solution of the linear part of the differential equation (3.2), i.e. 

( 3 . 11 ) 

[7(f)] ='7.(0 

the Laplace transform of H, , is 


//,( 5 ) = 1 /( 5 ^ + 2 ^ + 1 ). 

Similarly, equating the coefficients of gives 


(3.12) 



rj'i {t) + 2^7^ (r) + 73 (r) = 0 


(3.13) 


The above requires the second order kernel /12 to be identically to be zero, i.e. 

A2(ri,T2) = 0 ^ 2 . 14 ) 

Equating the coefficients of gives 

73 (r) + 2 ^ 7 ^ (r) + 73 (r) = (r) ( 3 . 15 ) 

Similar to that in the case of equation (3.10), the above requires 

73(r) = -l/:r,[7f(r)3 (3.16) 

Now, since 

%(r) = H,[fir)] (3.17) 

one gets 

= (3.18) 

and in terms of Laplace transforms (refer the block diagram representation of equation (3.18) 
in Fig. 3 . 2 ) 

= (3.19) 

where 'F3(5|,52»'^3)* christened as the synthesised third order kernel factor and is 

mathematically expressed as 

%is„s„s,) = -//,( j ,) 7 /,(^,)//,( 53 )//.(^. +^2 +^3)- (3-20) 


3.3 Wiener Kernel Representation of Response 

The practical application of Volterra series analysis is known to present two major problems, 
namely, measurement of individual kernels and convergence of the series. Wiener (1958) 




Fig. 3.2 


Block diagram representation for the operator, . 






suggested a set of orthogonal functionals in order to circumvent the above problems. If /(r) 
is a stationarj' Gaussian white noise excitation with variance A (which is same as the power 
spectrum, Sjj- (&>) , of input force for a Gaussian process), the nth order Wiener kernel written 
as 

>>^o{«);/(^)] (3.21) 

is an orthogonal functional formed from Volterra functionals such that it is orthogonal to any 
homogeneous Volterra functional of degree less than n, i.e. 

,^o(„)lf~(r)] = 0 (3.22) 


The system response can be readily written in terms of the Wiener functionals, as below (Lee 
et. al., 1964) 
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(3.23) 


whereby owing to the orthogonality property the Wiener and Volterra kernels can be shown 
to bear the following mutual relations (for a third order system response representation; 
Schetzen (1980)) 

^2(rj,r2) = M-2(r,,r2) 

Aj(r,) = M-i(ri) + Wi(3,(r,) 

Ay = n'l, + w„(2) 

with 
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3.4 Extraction of Wiener Kernels from Measured Response 

In view of the stated difficulties in the measurement of Volterra kernels, the measured system 
response is employed to extract the Wiener kernels. These Wiener kernels are then employed 
to generate the Volterra kernels, using the relationships of equations (3.24). 

Extraction of Wiener kernels from measured response involves enormous amount of data 
processing, since the kernels are multi-dimensional. The Laguerre filters proposed by Wiener 
(1958) or the alternative approach of using cross-correlation techniques and time delay filters 
suggested by Lee and Schetzen (1965), present formidable amount of data processing. Use of 
a complex filter, in the frequency domain (French and Butz, 1973), reduces the 
computational effort and is also suitable for such analysis since the Wiener kernel theory 
involves multidimensional convolutions. The scheme, employing a complex exponential 
filter, is graphically shown in Fig. 3.3. The Fourier transform representation for the variables 
in equation (3.23) may be e.xpressed as. 
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Substituting equation (3.25) into equation (3.23) and using the relationship that 
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one obtains 




{■'ig. 3.3 


Scheme for evaluation of the first order Wiener kernel transform, fV^ (cd) . 
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The output z(r) . from the exponential filter, can be written as 
z(r)= |c^'‘''^'/(r- r,)f/r, 

“ OC ' 


The ensemble average of the output of the circuit (Fig. 3.3) is obtained as 

{rj{T)z{T))^(F\co^e-j‘^^WQ 
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(3.27) 


(3.28) 


(3.29) 


2 


Since /(r) is stationary Gaussian white noise with zero mean and variance A, the Fourier 
transform. F(d)), is also a stationary Gaussian white noise process and the ensemble 
averages of the products transformed functions are (Raemer, 1969) 



(F(fy)) = 0 

(F(<y,)F(G) 2 ))= + 6 ) 2 ) 

(F(6>,)F(02)F(6)3)) = 0 

(F(<a,)f(®2)n‘»3)f(®4))=^^[5(ffl, +<»2)#(ffl3+ffl4)] + ^^[^(ffl, +m,)S{ci>2 +< 04 )] 

+^^[5(0), 4-(i),)] (3.30) 

The relations of equation (3.30) reduce equation (3.29) to 

{tj(t)z(t)} = AH^,(a)) (3 31 ) 

However, due to the equivalence of time and ensemble averages, the ensemble average 
can also be written as 


1 r /2 

(;7(r)z(r))= lim - | T]{T)z{r)dr 

T-*x T -TI2 

giving 

(T}(T)ziT)) = F*(co)T]ia}) 


(3.32) 

(3.33) 


Equations (3.31) and (3.33) give 
AW^ia)) = F*i£o)rj{co) 

from which the exprcs.sion for the Fourier transform of the first order Wiener kernel is be 
obtained as 

W^{a)) = F*ieo)TjiQ})/A (3.34) 

For measurement of third order kernel transform, a circuit invohing three exponential delay 
filters as shown in Fig. 3.4 is considered. The output, z(t) , from the exponential filters is 

z(r)= (3.35) 

-»oo ^ 

The above, after some algebraic manipulations, reduces to 
z(T) = Fi-m^)Fi-a)2)Fi-co^)e ^ 2 3 




Fig. 3.4 Scheme for evaluation of the third order Wiener kernel transform, (o), co, w) 






The ensemble average of the output of the circuit can be expressed as 
{7(r)7(r)> = A^U\(a>0S(-co2 - 6 ) 3 >•'■^"‘“ 2 -^ 3 )- 

+ n\ico2)S(-o>y 

+ W',(6>3)<?(-6>, ]+6A%(6)y,0>2,CO3) 

The equivalence of time and ensemble averages gives 
(7(r)7(r)) = lim - irj(T)z(T)dT 

J / ^772 


(3.36) 


(3.37) 


giving 

{r]{T)ziT))= Fi~0)^)F(-C02)F(-(i)2)rj{Q)^ +6>2 +fi>3) (3.38) 

Equations (3.36) and (3.38) give the expression for the measurement of the Fourier transform 
of the third order W'iener Kernel as 

,ru2.o>3) = (1 / (>A^)^F* {coy)F* {( 02 )F* {ci) 2 )ri{co^ +6)2 + m3)j 

- (1 / 6 v 4 )[W^, {cox)d{a >2 + CD^) + l^\{ 6 } 2 )S{o)i +6)3) + fF, (^>3 ) J(£u, + 0)2)^ 

(3.39) 

,^ 2 .ry 3 ) forms a multi-dimensional surface on the (tu, ,< 2 ) 2 »<y 3 ) axes. 
Measurements are made for special trispectral kernels with coi = CO 2 = cdj = co . These 

kernels are functions of only one variable co and are much easier to compute and interpret. 
(Bendai (1990) has termed such single function transforms for Volterra series as Special 
Trispectra! Kernel Transforms). The Special Trispectral Wiener Kernel Transform can be 
readily written from (3.39) as 

»7o;,r^,^.ri = (l/6^^)[{F*(fij)}^77(36))]-(l/24r,(cy)J(iu)] (3.40) 


3.5 Parameter Estimation 

The Special I rispectral Wiener kernel Transforms can be extracted from the measurement of 
the applied randont force and system response and employing equations (3.34) and (3.40). 


Subsequently, for a third order representation of the system response, noting the equivalence 
between the Volterra and Wiener kernels (Eqn. 3.24), the third order Special Trispectral 
Volterra kernel transform can be computed from 

H^{Q},0},C0) = W^((O,C0,Q>) 

= (I/6a4^)[{F*(6?)}3 7(3©)]-(l/2a4)[r,(G>)fy(e>)] 

Similarly, from the relations of equations (3.24), the expression for the first order Volterra 
kernel transform, in terms of the measured Wiener kernel transforms, becomes 

//,(&>) = JTj (6?) + (3 42) 

where lf,(<»)is as given in equation (3.34) and lF,( 3 )(<o)is the first order derived kernel 
transform, 

00 

^^i(3)(<y) = -3^ lW^{co,0)2-a)2)da)j^ (3.43) 

-00 

The linear parameters, cOfi and ^ can be readily obtained by equating the measured first 
order Volterra kernel (of equation (3.42)) to its analytical expression given in equation 
(3.12). Standard curve fitting techniques can be used. These, estimated linear parameters, 
m fj and ^ , are employed further in the estimation of the nonlinear parameter A . 


The estimate for A is obtained by equating the synthesized expression (Eqn. 3.20) and the 
measured value (Eqn. 3.41) of the third order Volterra kernels of the system. Thus 


(1/6.4^); 




■ (1/24)^, Wwlyt'P, (<»,<»,(»)] 


(3.44) 


3.6 Computer Simulation 

The procedure is illustrated through numerical simulation of the response of the 
nondimensional equation with cubic nonlinearity, Eqn.(3.2). Owing to the statistical nature of 
the estimation procedure, illustration is carried out for various values of the nonlinearity 
parameter, X . The procedure is repeated for various values of the damping ratio, ^ . 



The forcing function in equation (3.2) is a normalised random force, /(r), with zero mean 
value. The excitation force is simulated through random number generating subroutine and is 
normalised with respect to its maximum value . A typical sample of this excitation 

force is shown in Fig. 3.5(a), while Fig. 3.5(b) shows a typical nondimensional response, 
Tjir ) , for nonlinear parameter /I =0.1 and a damping ratio, ^ = 0.01 . The response has been 
numerically generated for 4096 number of instances in the nondimensional time ( r ) range 0- 
2048, using standard fourth-order Runge-Kutta subroutine. The response is computed for 
2000 number of samples of the simulated random force (The influence of the number of 
samples in the ensemble is discussed later). 


3.6.1 Estimates of Linear Parameters 

Fig. 3.6(a) depicts the power spectrum of the random force averaged over the ensemble of 
2000 samples. The corresponding ensemble average of the power spectrum of the response 
can be seen in Fig. 3.6(b). The first order Volterra kernel, H^ico), Fig. 3.7(a), is then 
computed, using the expression of equation (3.42), over the ensembles of the force and 
re.sponsc. The linear parameters o)„ and ^ are computed from //,(ru), through routine 
modal analysis procedures (Ewins, 1984). The linear parameters, thus estimated are 

0.1 606 (cycles / r) = 1.009 (rad / r) 
and ^ ~ 0.012. 


I'he curve of Fig. 3.7(b) shows the error incurred in the estimate of /f, (ui) , due to the 

slatistical nature of the Fast Fourier Transform computational procedure and the finite length 
of an individual sample (4096 in the present case). The normalised random error, as known 
(Bendat and Piersol. 1986). can be seen to be the maximum in the vicinity of natural 

frequency of the linear part of the system (Irad / r i.e. 0.1592 cycles / r). The error in th 

frequency range 0.0 - 0.10. is less than 4%. It can be readily inferred that the normalised 
error for the higher order kernels would show a similar trend and the error in the estimate of 
the nonlinear parameter 2 can be expected to be less in the frequency zone 0.0-0.10 cycles / 
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Fig. 3.6(a) Power spectrum of the input force (averaged over 2000 samples). 
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Fig. 3.6(b) Power spectrum of the response (averaged over 2000 samples) 

(for simulation values of A = 0.10; 4 - • 
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3, 7(a) Estimate of first order Volterra kernel transform, Hi (cd) 
(for simulation values of A = 0.10; ^ = 0.01) . 






3.6.2 Estimates of the Nonlinear Parameter 

Nonlinear estimation is carried out for a range of values of the nonlinear parameter Z , and 
damping ratios The response of the nondimensional equation (3.2) is numerically 
simulated for A = 1.00, 0.10 and 0.01, while keeping the damping ratio fixed at 0.01. It 

may be noted from equation (3.3) that A includes both, the nonlinear stiffness term and 
^max ■ maximum value of the applied force. A low value of ;i , for a fixed value of 

implies a low value of , while a high value of A , for the same implies a high Fj,^^ 

and vice-\ ersa. The results for such a nondimensional parameter can be readily employed to 
design experiments and decide the excitation level for ah expected nonlinearity of a given 
system. 

7'he estimated results for the third order kernel are depicted in Fig. 3.8. The third order kernel 
factor synsthesised from the measured first order kernel transform /f, (ru), and 

the measured third order kernel transform H^{co,a),(o) are shown in Fig. 3.8 (a), 3.9(a) and 
3.10(a) for ?. = 1.00, 0.1 and 0.01 respectively (^= 0.01). It is to be noted here, that while 
the first order kernel is estimated in the entire available frequency range 0.0 - 1.0, the third 
order kernels, involving a 3a> factor, have to be restricted to one-third of this frequency zone 
(i.e. 0.0 - 0.33). It can be observed, from the figures, that while the measured third order 
kernel traii^lonn //^(ry, <w,£y) is reasonably accurate in showing the harmonic at /3 (at 

nondimen-sional frequency = 0.053), the identification of the harmonic at 0)„ (at 
nondimensional frequency — 0.159) is weak, the best approximation being in the case of A — 
1 .0. The estimation of A , from these kernels, is therefore restricted to the frequency zone of 
0.0 - 0.10. The estimates of the nonlinear parameter A, obtained in accordance with the 
relationship (3.44), are shown in Figs. 3.8(b), 3.9(b), 3.10(b). A fourth order polynomial 
curve regressed through the estimates of A, over the frequency range, is also shown in 
Figures 3.8tb). 3.9(d) and 3.10(b). The mean values of the estimates of ;i are found to be 
1.08, 0.14 and 0.05 respectively. The order of the magnitude can be seen to be estimated 
correctly in all three cases, while a good accuracy can be seen to be obtained for A =1 .0 (Fig. 
3.8(b)). It is to be noted here that the response representation of equation (3.23) has been 
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restricted to include kernels only up to the third order, in order to keep the computations to a 
manageable level. Inclusion of higher order kernels (5**’., 7*^ order) in the response 
representation can be expected to improve the accuracy of the estimates at increased 
computational effort. Another source of inaccuracy in the estimates is the finite length of 
samples and the ensemble size. The nondimensional time interval for sampling has been 
taken as 0.5 and 4096 response instants are collected for a sample data set which gives a 
frequency bandwidth of ± I.O cycles / r and a frequency resolution of 0.488 x lO’' cycles / 
T . Additionally it may also be noted that the numerical fourth order Runge-Kutta procedure 
of response simulation is also a source of error. 

3.6.3 Influence of the Number of Samples 

It was observed that while increasing the sample size beyond 500 has insignificant influence 
on the first order estimates, the influence on the third order estimates can be seen from the 
curves of Figs. 3.1 1(a) - (d). The figures show the power spectrum of the input for ensemble 
sizes 500, 1000, 1500 and 2000 along with the corresponding third order kernel factor 
and the measured third order Wiener kernel transform, H 2 {o), 0 ), 0 )) , which are 

both observed to get refined with increasing number of samples in the ensemble. In the 
{iresent study the ensemble size has been limited to 2000. 

3.6.4 Influence of Measurement Noise 

I-rrors can be expected during an experiment, in the measurement of the excitation force and 
the response. The influence of measurement noise is studied by contaminating the simulated 
force and response signals, individually and simultaneously, with 5% simulated random 
noise. The frequency range is split into two and Fig. 3.12(a) shows the third order kernel 
transform H and the third order kernel factor 4 ^ 3 ( 61 , 61 , 61 ) (for 5% noise in both 

input and output), in the frequency zone of interest 0.00-0.03 cycles / r , on a magnified 
scale, while the remaining portions of the curves are shov,n in Fig. 3.12(b). Other estimates 
are also observed to be similarly robust to measurement noise influence. 
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Fig. 3.12 Effect of measurement noise on third order estimates 
(for simulation values of A = 0.10; = 0.01) . 






3.6.5 Influence of Damping on the Accuracy of Nonlinear Estimates 

Apart from the nonlinear parameter X , the other nondimensional parameter contained by the 
governing equation (3.2) is the damping ratio ^.Numerical simulation is carried out to check 
the accuracy of the estimation procedure and the response is simulated for tv^o different 
damping values ^ = 0.01 and 0.001 (for X- 1.0). Figs. 3.13(a,b)-3.14(a,b) show the 
estimates in the two cases. The estimates of ( 1 .08 and 1 .72, respectively) can be seen to be 
sensitive to damping being accurate for higher damping. 

3.6.6 Identification of the Sign of Nonlinear Parameter 

In the results depicted so far, nonlinear estimates have been arrived at by dividing (refer 
equation 3.44) the absolute value of the measured third order kernel transform, 
H-^{a),co,a }) , by the absolute value of the third order kernel factor, Both these 

terms arc, however, complex quantities and their absolute values have been considered while 
estimating the nonlinear parameter X , solely for achieving better smoothening of data. This 
involves the loss of information on the sign (whether positive or negative ) of A . The 
nonlinear parameter X being a real constant, information on its sign along with its 
magnitude, can be readily obtained if the real part of H^{(0,(0,(d) were to be divided by the 
real part of %((-u,fy,6t) (or the imaginary part of H^{o),co,co) were to be divided by the 
imaginary part of 'P,(<y,e;,fy) ). Consideration of equations (3.19) and (3.20) reveals that 
T,(6>,ru,r(;), a term synthesised from the first order kernel transform is a function solely of 
the linear parameters r^„and ^ and will be the same whether the nonlinearity in system is 
(■+ X ) or (- A ). The information on the sign of A is contained in Hi{co,co,co ) , which shows 
a reversal of sign, with change in the sign of A . This is illustrated in Figs. 3.15-3.17. The 
real and imaginary part of the synthesised third order kernel factor %{a>,a>,co) are shown 
Fig. 3.15(a), (b). The real and imaginary components of the measured third order Volterra 
kernel transform for a positive nonlinearilty (i.e. +A) are shown in Fig. 

3.16(a), (b), while those for a negative nonlinearity (i.e. - A ) are given in Fig. 3.17(a), (b). 
The peaks, occurring at a frequency of 0.056 cycles / r , can be seen to bear similar signs 
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(Figs. 3.15, 3.16) for and H^(^co,co,co) , for positive nonlinearily. However, 

they bear opposite signs for negative nonlinearity. Figs. 3.15 and 3.17. 

TTie results shown in Figs.3.15-3.17 pertain to nonlinearily X equal to +0.01 and -0.01 and a 
value of damping parameter ^ equal to 0.01. Similar results are obtained for other values of 
the nonlinear parameter and damping ratios. 

3.7 Remarks 

The procedure developed gives good engineering estimates of the nonlinear parameter. The 
estimates are satisfactory for a range of system damping. It also appears to be robust to 
measurement noise. The analysis presented is in nondimensional form and can be suitably 
employed to design experiments. The accuracy of the estimates, shows improvement with the 
increase in the number of force and response samples, over which averaging is carried out. 
Extension of the present procedure to systems with more than one degree of freedom and 
involving coupling is explored in subsequent chapters. 



CHAPTER 4 


PARAMETER ESTIMATION IN RIGID ROTORS 
SUPPORTED IN BEARINGS WITH CROSS-COUPLING 

The parameter estimation procedure, becomes more involved, if the cross-coupling effects in 
the bearings are taken into account. This may be the case for a rotor supported in fluid film 
bearings. In contrast to the no cross-coupling case, which could be treated as a single-degree- 
of-freedom system, bearings with cross-coupling pose a nonlinear two-degree-of-ffeedom 
problem. It is a system where an excitation, say, in the x-direction causes motion in the both 
the X and jp-directions. Such a case, therefore would require cross-kernels to be defined, 
which can be convolved with the excitation for the response representation, in terms of 
Volterra or Wiener series. 

Starting with the equations of motion, Laplace transforms are employed in this chapter, to 
derive expressions for the first and higher order direct and cross-kernel transforms. The 
algebra is kept simple by taking the damping to be linear and without cross-coupling. The 
nonlinearity is taken to be cubic and analysis is carried out in the frequency domain. First and 
the third order direct and cross-kernels are extracted from the measurements of the excitation 
force and response. Third order direct and cross-kernel factors are synthesised from the 
measured fiirst order kernel transforms. These synthesised third order factors are then 
processed along with the measured third order kernel transforms for estimation of the 
nonlinear parameters. The procedure is illustrated through numerical simulation. The 
estimation is carried out for various values of the nondimensional nonlinearity parameter and 
damping ratio. Statistical considerations and effect of probable measurement noise are 

illustrated and discussed. 

4.1 Governing Equations and Volterra Series Response Representation 

The equations of motion for a twodegme-freedom idealisation of the rotor-bearing system 
shown in Fig. 4.1 are 




Rgure 4.1 Rigid rotor in bearings with cross- coupling. 


mx + c„x + k^x + k^y + k^x^ = f^{t) 

my + .V + kyyy + ky^x + k^y'^ + = f^ (t) 

In the above equations k^^, kyy are the direct linear stiffness terms, k^,k^ are the cross- 

coupled stiffness terms, k^,k^ are the direct nonlinear stiffness terms, while k^.^k^ 

represent the nonlinear cross-coupling terms. The damping is linear and comprising of direct 
terms c . c . The excitation forces in the x and 3/ directions are /, (t), f 2 (0 respectively. 


Defining 


T = p.t, p = ^ky^ / m 

Xr) = x/;f,, , yri(j)=^ylX,p, (4.2) 

^ xt ~ ~ f ^ ’ / — 1,2 

and substituting in (4.1) the equations of motion can be written in nondimensional form as 


r;'(r ) + 2^ ,, I7’(r ) + x'- ‘T](r) + Ai ,. + x'L’7]^{z) = 


where 


p - J} 

2mp 




kn 


^XX 


k^ 

_ "■y •* max 


x,y; j^x,y 


(4.3) 


(4.4) 


^XX 


The solution of equations (4.3) is represented in terms of Volterra operators as 

'';7(r)='//[/,(r),/2(r)] 
where, (with /c denoting x or y) 


*7/[/,(r),/2(r)]=^/:ro+ 

/= 1,2 /= 1 , 2 ;;= 1,2 

+ y; [/, (r), /; (r), A (r)] + 

/«I,2;y=l,2;^=t,2 


(4.5) 


(4.6) 



The individual operators of equation (4.6) are given in kernel form by 

oc 

for/ = 1,2 




2 (^1 )fM - ^2 )dr^dr2 


(4.7) 


for i = 1,2; j = 1,2 

-00 

X dTydT 2 dT-^ for i = l,2; j = \,2\ A: = 1,2 


00 00 00 


-00 — 00 —00 


with nth order Volterra kernel 

K{Ty rj= 0 


for r_<0, p= 1,2,3,. 


.n. 


(4.8) 


For convenience, writing the Volterra operators as 




the response of equation (4.5) can be written as 

‘m- I'-?""’ + 

,.1.2 /= I .2;;=1,2 ,=1,2;;= I ,2;/:=1,2 

Z>1"+ 

r ..).2 l>^ 1 . 2 ;y*I ,2 /=1.2;y=l,2;*=1.2 


(4.9) 


(4.10) 


4.2 Synthesis of Higher Order Volterra Kernel Factors 

The Volterra operators are now determined as follows. The excitation forces /i(r), / 2 (r) 
are replaced by c/,(r) and c/^dr) respectively, e being a constant Noting equations (4.7), 
the resulting response of the system, becomes 

,= 1.2 /=1,2;;=1,2 

•>(r)= Zc-'V,”- 



Z,3 

1-1,. 

2;;=1,2;1:=1,2 

,<'•» H- 

Z,3,,(..M) 






/=: U 2 ;/= 1 . 2;^=12 


(4.11) 



Substituting equations (4.1 1) and the derivatives in equations (4.3) one gets 
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Summing up the responses of equal order, one can write 


I or der response 


/=1,2 




'7, = E"’7! 

z=l,2 


(4.13) 


II order response 
/ = 1.2;y=1.2 

^'^2= 

j=.I,2;y^l,2 


(4.14) 


III order response 
''73= E'#"*’ 

i=K2;y=l,2:A‘«I»2 

■'■'73= E-''7;'-"‘’ 

;=^1.2;y=l,2;A' = 1.2 


(4.15) 


Noting the following symmetry of kernels (Schetzen, 1980) 


x„^Uj) _x (,hi) . 
V2 “ V2 ’ 

”2 “ ”2 ’ 


X „U.J..i)_x _x UJ>0 

V3 - V3 - V3 

V3 “73 - 73 


; etc. 


(4.16) 


and u.sing equations {4.13)-(4.16), equations (4.12) are written in condensed form as. 


00 

Ic' 

7=1 


l + ;i^ 

•r 

00 

Sc" 

X ^ 

In 

3 

+4 

00 

Sc" ^77„ 


_« = i 



.7 = 1 


-i3 




00 


'''7;+4 

/2 = 1 


L X 
yx 


ln\ 


] + /l^ 

00 

Ic" 

1 


00 

Sc r]„ 

_n=\ 

j 




= cf 2{.r) 
(4.17) 


Equations (4.17) are power series in c with coefficients of c” being r]„ or ^rj^. The 
responses, ^ and ^ are determined by equating the like powers of c as follows 


c’ terms : 

7i ^rjy = /,(r) 

^'^r+ 2 ^,,, ^'; 7 i' + 4 ' '> 1+4 ""^1 = / 2 (^-) 

c^ terms: 

''72+2^xx ''72+''72+4 ^^72=0 

^^2 + ^'ri2 + 4 ^ni + 4 "^^2 = 0 


3 

c terms: 

^73+^73 +4 "^3 +4 V+4" 
+7; + 2^,, >'73 + 4 "^3 + 4 '^3+4 ^7? + 


7 ? =0 

4 


0 


( 4 . 20 ) 


Equations (4. 1 8)-(4.20) can be solved sequentially. Taking Laplace transforms of equations 
(4.18), for zero initial conditions, one obtains 


2 ^ 7, ( 5 ) + ' 7, (5)+^ 7, (s) + 4 ^ ^1 i^)] = 0 

■’' 7i (^) + 2^ ,, 5 7i (^) + >' 7, ( 5 ) + 4 ^ 7, ( 5 )] = F 2 (s) 


Solving the above two simultaneous equations, the solutions for ^ 7, ( 5 ) and ^rji (s) are 


, (^'+2f„,s + 4)F,(.)-4F,(^) 

7i (^) = "T”' "" — 0 1 T~T~ 

+2^„j + 1Xj' +2i„s + A’-^)-(X%x‘-^) 

>rj (5) = SL_ f. 

(;y ^ + IX^ ^ + l^yyS + 4 ) " (44 ) 


(4.22) 


Referring to the notations (4.9) and (4.13), the individual Volterra operators, from the above, 

are 





(s + 


(S^ + + l)(i 2 + 2{^,s + A^,) - ) 


^Hl^Hs) 


%■ 


( 5 ^ + 2<^^,S + 1 )( 5 ^ + 2^y^,S + ^^yy.)-iZ^^2L) 


^'Hl^\s) 


-A 


yx 


(5 ' + 2^,,. + l)(s^ + 2^ ^s + ) - (2^4 ) 


(4.23) 


yii{2)( s_ -^2^ ^s + \) 

‘ ( 5 ' + 2 ^^. + l){s^ + 2^^,s + 4) - 

For the second order kernels, the Laplace transforms of equations (4.19), similarly give 

"r7,(5) + 2^„5 "t]2{s)+"7]2(s) + A^ >';7^ (5)1 = 0 

^ T (4.24) 

[5^ ■ ;7;(5) + 2^,,.5 ''tJ2(s) + A^'^ ^’tJ2{s) + A^^^ ^;;2(-S)] = 0 


which yield 
'TJ2is)-- 0 
>'77,(.v)=0 


(4.25) 


that is. the second order kernel is identically zero - 

) " 0 (4 2 

■'772(r,,r2) = 0 

In order to synthesise expressions for third order kernels, equations (4.20) can be written as 


'7;' + 2.;„ ' ri-f + X^y’’ rj^ =9, 

■' rj'i + 2f „ ■' 7 ; + X% ’'!?, + X';^ ^ m = <!2 

where the following abbreviations have been used 

<7 I ~ A 77 ) A Jiy Tl\ 

qj = -A „ 7j Ay^ 77] 


(4.28) 



Equations (4.27) are linear in 7 / 3 , similar to equations (4.18) and therefore the solution in 
terms of Volterra operators is 

I?, 


that is 

However, since and ^ are linear operators, one obtains 

(4.31) 


Noting the abbreviations (eqn. 4.31), the terms on the right hand sides of the above 
equations, are individually expanded as 

(4.32) 

' -i;, ’//<"( '77?]-^"; '7?1 


In the abo\ e equation.s 

>i = 


’W|"'[/l(0]4 W,‘"’|/2(Ol 


- 7i + 7i 


Taking Laplace transforms of equations (4.32) gives 


+ +3 + 3 1,2,2 ^^vj/j-2,2.2^ 


XX VJ,^ 


1-1,2,2 , xxu/1-2.2,2' 


(4.33) 





In equation (4.34). the terms and have been 

called in this work as third order kernel factors and are defined as 


”4']' +42+43) '<’(4,) Vf'Usj) 

'> T]'-'" ‘'(s,,42,ij) = - '//{'•(s, +42 +43) ^«,">(4,) ''Hf4>(42) (4.35) 

/=U;3-1,2; *-1,2; I -1,2 

The third order kernel factors above, can be readily constructed, using equations from the 
first order kernels 

Similarly, the Laplace transforms of other terms on the right hand side of equations (4.31) 
can be worked out to be 


)] = (■’^■Vi/2-l.U ^3 ;0'V|,2-].1.2 ^3 x>vj/2-l,2,2^xyxj,2-2,2,2^ 

+ +3 _^3 xxvj,^2-l,2,2_^xxvj/2-2,2,2^ 

•'7/;'>(Pi{5i. 5,,S3)| = +3 ^VJ/I-1.1.2 ^3 >arvj,j-l,2,2_^>t3(vj/j-2,2.2^ 

+3 +3 }y^^<2,2_^yy^]-2,2,2^ 

■"i7;^>l^2(5,,52,53)l = (2^v{/2-UJ +3 W'q/2-U2^3 >yvj/2-l,2,2^XVvj/2-2,2,2 ^ 

+ A^ (yxvj/2-1,1.1 _|,3 >-Jvj/2-1.1.2 ^3 >!Xxj/2-1,2,2_j^>«v}|2-2,2,2^ 


(4.36) 


(4.37) 


(4.38) 


In equations (4.36)-(4.38), the following third order kernel factors have been used in addition 
to tho.se defined in equation (4.35). 


■'’•'fl' ^■‘■''(S,.42.4 j) = 
^'T''^'-''‘’(4„ 42,43) = 


- >'«!"(4, +42+4,) ^//f*»(42) >'i/f'>(4,) 

- +42+43) ''//P’( 42 ) 'Wf'*( 43 ) 

i = l2;J=l,2; k=l,2; 1=1,2 


(4.39) 


The Laplace transforms of the third order kernels (equations (4.30)), can now be expressed as 



’J^3(S,.S..S3) = ;i'^ +3 

+ /I"r4'3'-'''-' +3 


xx\j/j-U,2 «vp^l-l,2,2_^»:vj/J-2,2.2^ 
Jrvvj/ W,I,2 ^ wv{/j-l,2,2_^jrxxpl-2,2,2^ 

^'VJ/2-I,1,2 jr>'Vj/2-l,2,2_j_j^^2-2,2,2^ 




1,2,2 _^jcc 


Vj/2-2,2,2 ' 


>7/3(5,,.V2.53) = A;^, 


+ 3 ^3 w\j/^-l,2,2_|^>«xj,j-2,2,2^ 

+ 3 +3 y}' yj/ 1-12,2 1-2,2, 2^ 

+ 3 XVv^S-l.U _|_3 w\j/2-l,2,2_^>yvj/2-2,2,2^ 
+ 3 >''^'J/2-I,1,2 ^3 >^\j/2-1,2,2^>«vj/2-2,2,2^ 


(4.40) 


4.3 Measurement of Wiener Kernels 

As stated in Chapter 3, measurement of individual Volterra kernels is not possible, while 
equivalent ’^'iener kernels can be extracted from the measured response if the excitation to 
the system is w'hite and Gaussian. These Wiener kernels can then be used to generate the 
Volterra kernels. 


In the present case, the Wiener kernels of the nonlinear system are extracted by application of 
white Gaussian forces /i(r) and one at a time, i.e. first a white Gaussian force 

/, (r) with variance , is applied in the x-direction, while keeping the >'-direction force, 

resulting response (in both the x and _y directions), is employed to extract the 
direct, x-direction, Wiener kernels and the cross (x^y) kernels. In the next instance, a white 
Gaussian force /2 (r) with variance A 2 , is applied in the ;;-direction, while keeping the x- 
direction force, /i(t")= 0- The system response, in this instance, is employed to extract the 
direct >»-direction Wiener kernels and the cross (y-x) kernels. 


The system response, in terms of Wiener kernels, is now expressed, in the two individual 


cases as - 



(4.41) 


/, (r) Gaussian white, with variance and /^(r) = 0 : 

';?(r)='»'[/,(r). 0]='»|/,(,-)] 

>' 7 (r)=' »|/,(r). 0]=‘'»|/|(r)] 

/2 (r) Gaussian white, with variance and (r) = 0 : 

';7(r)='r[o,/3(r)]='r[/2(r)] 

'7(r)=’»'[0,/,(r)p(f[/2(r)] 


In the above equations 

(with K denoting x or j') for / = 1 or 2 


(4.43) 


I he individual ojx’rators of equation (4.43) being given in kernel form by 
(01- J"'.;"(r,)7;.(r-r,Mr, 

ur y 00 

-x ‘f I* ~oo 

y / / 

I J j'''^f(rj,r2,T3)/i(r-Tj)/,(T-T2)/,(T-T3)dridT2dr3 

' / - - y 

y » 

-3A, fj*Hf(rj,T2,T2)/,(r-Tj)dTidr2 (4.44) 

The relationship between the Volterra kernels of equations (4.7) and the Wiener kernels 
above, as in the pre\ious chapter can be shown to be (for a third order response 
representation ) 



(4.45) 


hj^^l ’ ^2 )“*^^2 ’ ^ 2 ) 

"A!'kr,)="H-l‘)(r,)+"wi;])(^i) 

^■/zo='^Wo 4 ^h-(>( 2 , 

with 

H-;;J) (r, ) = -34 ]"■ wf (r, , r2 , r2 )dT2 

-OC' 

""^0(2) ="^< J*^H4'^(r,,r,)^/r, 


Employing the following Fourier transforms, 


/,.(r)= J/;(r^)e^“'Vru 

»«or 


^;7(r) = 


I *■ ri((t))e '“’^' d(0 


M'j" ( r ) =-" I ^ (<y , )eJ^ » ^do) 1 

-.jr 

'^M4'kr|.r2)= J 

•hx -“or. 


^ M'3^’(fl.f2-^3)" 


--€if . »-OC. -*00 


(4.46) 


and using the relationship 

= <5t<w) , the Dirac Delta function, 

-OC 


(4.47) 


the responses (4.41 ) and (4.42) can be expressed as 



— OO 

qO OO 

-OO — OC) ^ 

— OO 

Xj OO OO 

-OC ~€f^ -0C» 

Of) 00 

-3Ai J J^JF3'^(<5),,<y2^-^2W^l<^<^2 


+ , 


(4.48) 


+ 


'r)(r)=’'fVc + 

-00 

00 00 00 

J I (a, . ©3 )f;. (CO, )F, K )e *“ 0 'dco,dm^ - A, \ '-r W (fflj -o, )d ®3 

-00 - OC ^ -00 

X.) 00 00 

4 J J \yWl'\o),,C 02 . 0 }^)Fi{( 0 {)F.X( 02 )F,i(D^^^^^ 

-OC* '“00-00 

00 00 

- 3 A, J |'*'l^^ 3 '^{<y,,iy 2 


(4.49) 


A complex exponential filter, similar to the one described in Chapter 3, has been used for 
measurement of the individual Wiener kernel transforms. Referring to Fig.4.2, the ou^ut of 
the filter is 


2 3 * ■ 

z^'^(r)= fc-'‘"’'‘/,(r-rj)r/rj 

-«00 

and the ensemble averages of the outputs of the circuit are, 


67 




Fig. 4.2 


Scheme for evaluating the first order direct and cross Wiener kernel 
transforms, If/' ^(6)). 





or ex' _ 

i i )(^- (®i )r; K )r;’ j 

- A,{f- (< o))e-i’^' '‘rw<‘\m„-ci,,)da,, 

—00 

or oc oc ^ 

1 J J''^3'^(ryi,ry2,ry3)(^(^y,)^(6)2)^.(6j3)^*(fi,)\ 

-a - oc -oci ' / 

34 / -®2)(f;(o,);;*wy-><®i 


(4.51) 


1 )>IVl'\cOt.a>yF,(<a0FA‘^2)P'(‘^)y'"^^^^^ 


- A,{F;(m))e-^‘^^ j>'Wj‘\m„-o,2)d02 


) "j )>Wi‘H0„0^,0,){FMme>2)F,(‘»,)F;(0)) 


X 


d0,d02do>2 


34 "j J’'W',‘'*(<B|,®2.-®2)(^(®l)^*Wy 

«ar ^ ‘ 


(4.52) 


Since /,(r) is slalionar>- Gaussian white noise with zero mean and variance 4- , the Fourier 
transform. F^io)), is also a stationary Gaussian white noise process and employing 
equations (3.30) from the previous chapter for ensemble averages of the products 
transformed functions, equations (4.51) and (4.52) can be reduced, after some algebra to 



(4.53) 

(4.54) 


(^^ri{T)z'-'\T))= A, ^W}'\co) 

{>'rj{r)z^‘\T))= A,ywl'\co) 

However, due to the equivalence of time and ensemble averages, the ensemble average 
r 7 (r)z''*(r)^ can also be written as, 


= lim ~ J ^Tj{T)z^‘\T)dr 
\ I T-*^'T-ri2 

= F*{co)'‘rj{Q)) (4 55 ) 

Equations (4.53) and (4.55) give 

A^’nf\o»^F'(a) Xffl) (4.56) 

from which the expression for the first order Wiener kernel transform is obtained as 
"W'^‘\(o)^^F*{(i)) 'vCru)/^^ (4.57) 

Similarly, 

ywl^\a>) = F*i(0) yjmiA. (4.58) 


Since / take.s values I and 2 the direct kernels yw^^\a>) and the cross kernels 

■‘ W'^'’(f«i) can be extracted from the measured responses ^//(^y) and 

the applied force ( Fi((o) and its variance A^ or F 2 (iy) and its variance A^ ) , through 
equations (4.57) and (4.58). 

For measurement of the third order kernel transform, a circuit involving three exponential 
delay filters, as shown in Fig. 4.3, is considered. The output, z^'\t), from the exponential 


filters is 









z«(r) = Jr«'3/,(r-r3)*3 

—00 ~oc 

= (4.59) 

and the ensemble averages of the outputs of the circuit are 

7(r)z^'^(r)^ = Af\^Wl'\(o^)5i-co2 -6>3)e-'(-'"2-‘“3)^ 

{coj )5{-(j) 1 - 6)3 


+^W^‘\co^)S{-cOi 


CA^fV^'^cOi, 0)2,0)^} 


However, the equivalence of time and ensemble averages gives 


^ 7]{T)z^'\r)) = lim - I ^ rj{r)z^‘\T)dr f 

= Fi-0)^)F{-0)2)F{-0)2>) "" 77(6)1 +6)2 +6)3) 


(4.60) 


(4.61) 


Equations (4.60) and (4.61) give the expression for the measurement of the third order 
Wiener kernel transform as 


6 A. 

I 


6 A. 


similarly we can have. 




6 A? 


6 A- 


F. (6),)f (0)2) f. (6)3) ^77(6), +6)2 +6)3) 

^nf\o)^)S((Oi + (i)^)+'‘wl'\o)2)S(a}i +6)3)+''lF/'^(£y3)J(6)i +6)2) 

(4.62) 

F* {a)i)F* i 02 )F. (6)3) ^77(6)1 +6)2 +6)3) 

>'r/'’((yj)<J(a), +6)3)+^Pr/'\6)2)<!)(6), +6)3)+^lf/'^(6>3)5(6)i +6)2) 


(4.63) 



The third order kernel transforms, , 6 ) 2 ,c) 3 ), ywl‘\a>^,co 2 ,co^) form multi- 

dimensional surfaces on the (fit)] ,<i >2 ><^ 3 ) axes. Measurements are made for special 
trispectral kemeks with iy , = ^2 = <^3 = ^ • As stated in the earlier chapter, these kernels 
are functions of onl j' one variable o) and are much easier to compute and interpret. For such 
trispectral kernel transforms the expressions (4.62) and (4.63), for their measurement, reduce 
to 









{a))S(o))+yw}''> (<u)<J(ru)+^lF/') (o))S(o)) 


(4.64) 


4.4 Parameter Estimation 


The third order Special Trispectral direct Wiener kernel transforms ^W^^\o},ci),a)) , 
(co,(i),(t}) and the cro.ss-kemcl transforms ^W^^\ci},a),(i>) , ^}V^^\o),cd,ci}) are 


extracted from the mea.surements of the responses ^r}{3o)), ^tj{3q)) and the applied force 
( /'l(fy) and its variance or Fjioj) and its variance ). Subsequently, for a third order 

representation of the system rcsp<mse, noting the equivalence between the Volterra and 
Wiener kernels (equations 4.45), the third order Special Trispectral Volterra kernel 
transforms can be computed as - 
direc t-kernels: 






^ ^ (ru, ry, ry )= •'’ rf * (<y, ry, ar) 

3 3 


tj{3o)) 


1 

6Ai 




1 






TjQeo) 


J. 

6A2I 


y {a))Sio})+ ^ ’ (u>)^ (fi>)+ ^ (<y)<^ (<y) 



cross-ker nels: 


(<y, CO, (o)=^ (fo, ( 0 , (o) 



1 

“ IaI 

{^^*{6>)p7(3<y) 

1 

6 A 2 

(a>)S (£y )+^ {co)S(o))+^W^^^'> ico)Sico) 

(a), a),co )= '' H ’ (co, co, co) 



1 

~ 6^,^ 

n (<y)f'''^(3iy) 

1 

6 Ai 

{co)Si0)+>w}^^ i(0)S(0})+ywl^^ (o})S(o)) 


(4.65) 

Similarly from ecjuations (4.45), the expressions for the first order Volterra kernel 

transforms, in terms of the extracted Wiener kernel transforms become 

direct-kernels: 

cross-kernels: 

(4.66) 

where and >'r« (re) , with i aking a value I or 2. are kernels derived from the 

third order kernels through the following expressions. 


00 


*00 

00 


2 2 ^ 




-00 


2 ' 2 ' 


(4.67) 


The linear parameters can be obtained from the above estimates of the first order Volterra 
kernel transforms, ^ hI^\q}), >‘H^^\co) . Noting • the algebraic 



expressions of these kernel transforms in equations (4.23), the ratios of the estimated kernel 
transforms can be expressed in terms of the linear parameters as follows 

The parameters 2'„,. A';,. 

^xx estimated through application of a curv'e 

fitting routine. 

Nonlinear parameters are computed from the estimates of the Special Trispectral Volterra 

Kernels, . Noting that these 

kernel transforms are estimated by application of a single white Gaussian force at a time, the 
expressions (4.40), for the third order Volterra kernel transforms, synthesised in Section 4.2, 
also get reduced, in the two individual cases, to 


/](r) Gaussian white, with variance ^4^ and f^(r) = 0 ; 






fj ( r) Gau-ssian while, with variance and /j (r) - 0 : 


(4.69) 


Substitution of the estimated third order kernel transforms on the left-hand side of the above, 
gives four equations, 



6A{ 




*7(3ty) 
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)j 'r}{3M) 
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W,^^\C0)S(C0)^ = 4 (-vj/I-U.l ) ^ (^xpl-l.l,! ^ 




- • ) -^ 

F^ (<y)| *7(3^) 


!^L 




6^f 


* I’ 

i (^i>)j "rj{3co) 


2A, 


+ 4 + 

= 4 ) + 4 ) 

+ 4 r'F2-2-2’2) + 4(«vp2-2.2,2^ 

+4 (-’''}^ 3 ^'‘’’’‘ + 4 (^%^~'’‘’') 


(4.70) 

which can be soh ed .simultaneously for the four unknowns, namely the nonlinear parameters, 
A'^' A'"' />'' 

^xx ^ ^ n ' n ' vx • 


4.5 Computer Simulation 

The procedure is illustrated through numerical simulation of the response for the 
nondinicnsional coupled equations (4.3). The forcing functions in the equations are 

nonnalised, zero mean random forces, /)(r) and / 2 (r). The excitation forces are 
simulated through random number generating subroutines and are normalised with respect to 
the maximum value of fi(t). A typical sample of such excitation is shown in Fig. 4.4(a). 
Fig. 4.4(b) shows the corresponding power-spectrum of the excitation force averaged over 
2000 samples. 

Owing to the statistical nature of the problem, the procedure is illustrated for various sets of 
linear and nonlinear stiffness parameters and damping factors. Various case studies have 



Auto-Power ^xsctrum {dB> ' Normalised force 



— ' ^ ^ ^ ^ ■ ! 1 ^ ^ U 

400 800 1200 1600 2000 

Dimensionless Time 


.4(a) Typical sample of input force. 



Fig. 4.4(b) 


Power-spectrum of the input force (averaged over 2000 samples). 



been designed to study the influence of these parameters and the errors involved. The 
follow'ing scheme is adopted. 

. Keeping the linear parameters, ^ (with X^y,=X%), fixed, the relative values 
of the nonlinear parameters and X^yy are varied. 

• Keeping the nonlinear parameters, and , fixed, the relative values of the linear 

parameters , X^^y , 4^, (with X^^^ = ) are varied. 

• For computational ease the nonlinear parameters X^^ and X^^ are taken as zero. 

yx 

• Damping is fixed and taken to be identical in the two directions, i.e. = 0.01. 

^ XX ^ yy 

(lire variation in the values of damping is taken up later) 

The response is computationally simulated by solving the governing equations through a 
standard fourth-order Runge-Kutta subroutine. These responses are fed as inputs to the 
parameter estimation algorithm. The estimated parameters are compared with those 
originally used for the simulation of response. 

4.5.1 Ca.se Study 1 (a) 

In the first instance, the nonlinear parameters are taken to be numerically equal, for response 
simulation, 'fhe direct linear parameters are also taken to be identical. The linear cross- 
coupled paramcter.s are taken as half of the direct linear terms. The numerical values are 
given below 

/i ( r) = normalised Gaussian white, as shown in Fig.4.4(a) ; /2 (G = 0 i 

A* =4=0.10 

4 ' 4 = O'® 

= 1.00 
= 0.50 
= 0.01 



With the above values, the simulated nondimensional responses, and ^ tj { t ) are 

typically shown in Figs. 4.5 (a) and 4.5 (b). These responses have been numerically 
generated for 4096 number of instances in the nondimensional time (r) range 0-2048. 2000 
number of such samples of response are obtained from 2000 different samples of the 
simulated random force. ( I he influence of the number of samples in the ensemble is 
discu.s.sed later) linsemble average of the power spectrum of the responses can be seen in Fig. 
4.6 (a) and 4.6 (b). 

These spectra are fed as inputs to the parameter estimation algorithm. 

E.stimat es of Lin ear Parameters 

The finst order direct Volterra kernel transform, and the first order cross Volterra 

kernel transform, ^ , Figs. 4.7 (a) and 4.7 (b) respectively, are computed, using 

equation.s (4.66). These kernel transforms exhibit peak responses at frequencies equal to 
0.112 and 0.195, which correspond to two critical frequencies (refer the denominators of 
equation.s (4.23)). 

=(1/2((4, +1)±((4+1)2+4(A^,-4,))'^2)jI/2 (4 71^ 

ITe linear paranteters 4 ’ » 4 ’ computed from these direct and cross 

kernel transforms, as explained in equation (4.68), through routine modal analysis procedures 
(Ew'ins, 1984). The linear parameters, thus estimated are 

4 “ 4 ==0.1 598 (cycles / r) 

= 1.0046 (rad / r) 

4 = 4 = 0.0799 (cycles / r) 

= 0.5025 (rad / r) 

=0.015. 



NondimeiuicMi>l response Nondimensional response 





Auto-Power Spectrum (dB> Auto-Power Spectrum (dB) 



Power-Spectrum of response in x-direction (averaged over 2000 samples): 
Case 1(a). 



Power-spectrum of response inj^-direction (averaged over 2000 samples): 
Case 1(a). 


Fig. 4.6(b) 





Nondlmensional Frequency 


Fig 4.7(a) Estimate of the first order direct Volterra kernel transform, {co) 
Case 1(a), 



Fig 4.7(b) Estimate of the first order cross Volterra kernel transform, ^ {co ) : 
Case ] 



The exact values of the above parameters are those used as input for numerical simulation of 
the response. While, linear stiffness parameters can be seen to be estimated with a good 
degree of accuracy, the error in the damping estimates is higher. Damping has been estimated 
by applying the standard half-power method on the first order kernels. The curves of Figs. 

4.8 (a) and 4.8 (b) show the errors incurred in the estimate of and due 

to the statistical nature of the Fast Fourier Transform computational procedure and the finite 
length of samples (4096 in the present case). (The exact values of these kernels are those 
obtained from the expressions (4.23), after direct substitution of the numerical values of the 
linear parameters employed for the simulation of the responses). The normalised random 
error, as known, can be seen to be the maximum in the vicinity of the peak responses at 
critical frequencies ^ ^ . Since the numerical error is higher in the vicinity of the peaks, the 

error in damping estimates, obtained by application of half-power method about the peaks, 
will be higher than those in the stiffness estimates. 


The first order kernel transforms and estimated by applying white 

Gaussian excitation in the y-direction, while keeping that inx-direction as zero, are shown in 
Figs. 4.9(a), (b) respectively. 


Estimates of Nonlinear Parameters 

As observed from Figs. 4.8 (a) and 4.8 (b), the statistical errors in the estimate of 

and (6>) , the error is less than 6%, in the frequency range 0.0 - 0.10 , for the ensemble 
size of 2(K)0. It can be readily inferred that the normalised error for the higher order kernels 
would show a similar trend and the error, in the estimate of the nonlinear parameters and 

, can be expected to be less in the frequency zone 0.0 - 0.10. 


The third order kernel 


factors, «).<»), 

(( 0 , 0 , 0 )), synthesised fix)m first order kernel transforms in 


accordance with equations (4.35) and (4.39) are shown in Figs. 4.10 (a)-(d). 
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I'ig. 4.8(a) Error in the estimate offirst order direct kernel transfono, . 



Fig. 4.8(b) Error in the estimate of first order cross-kernel transform, ^ (©) 
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I'ig. 4.9(a) Estimate of the first order cross Volterra kernel transform, {(o ) : 

Case 1(a). 



Fig. 4 . 9(b) Estimate of the first order direct V olterra kernel transform, {co) . 
Case 1(a). 




TMrd order keroel factor (dB| 


Nondimcnsional Frequency 


(a) Estimate of 6), 6y). 



0.01 0.1 
Nondimcnsional Frequency 


(b) Estimate of 


fiy. 4.10 


Third order kernel factors: Case 1(a) (Contd.) 



Tim4 kamel {Mmm (dB) Thkd order kernel factor (dB) 





The system response is then processed to extract the set of third order measured kernel 
transfornis "//f (^th /,(r)9t0 and / 2 (r) = 0) and 

y H^^\(o,co,co) (with / 2 (r)?i 0 and /,(r)=0). 

Howe\cr, since the nonlinear cross-coupling terms, are taken to be zero, only one 

of the above sets, oi third order kernel transforms, is sufficient for parameter estimation 
(refer equations 4.70). Presently, the first set ^(r»,6),ry), , as shown in 

Figs. 4.1 1 (a) and 4.1 1 (b), is extracted for use in parameter estimation. While the first order 
kernels are estimated in the entire available frequency range 0.0 - 1.0, the third order kernels, 
involving a 3rt> factor, have to be restricted to one-third of this frequency zone (i.e. 0.0 - 
0.33), as mentioned in the previous chapter. It can be observed, from the figures, that while 

the measured third order kernel transform ^ H 2 \o),(i),a)) is reasonably accurate in showing 
the harmonic at <i)„ ^ ^ / 3 (at nondimensional frequencies = 0.037 and 0.067), the 

identification of the harmonic at ^ ^ nondimensional frequencies = 0.1 12 and 0.195) 
is weak, due to higher statistical errors, mentioned earlier. The estimation of nonlinear 
parameters and , from these kernel transforms, is therefore restricted to the 
frequency zone of 0.0 - 0. 10. 

The estimate.s of the nonlinear parameter and X^yy , obtained in accordance with the 

relationship (4.70), are shown in Figs. 4.12 (a) and 4.12 (b). A fourth order polynomial curve 
regressed through the estimates of these nonlinear parameters , over the frequency range, is 
also shown in these figures. The mean values of the estimates are found to be 

= 0.098 

X^yy, =0.112 

The exact values of the above nonlinear parameters are those chosen for response simulation, 
that is, = /I'l = 0.10 . The nonlinear parameters can be seen to be estimated within the 




0.01 


0.1 

Nondimcnsional Frequency 


Fig. 4.1 Ha) Third order measured kernel transform, : Case 1(a) . 



Fig. 4. 11 (b) Third order measured kernel transform {(O, co, (d) : Case 1 (a). 





Fig. 4. 1 2(a) Estimate of the nonlinear parameter, : Case 1(a). 



Fig. 4.1 2(b) Estimate of the nonlinear parameter, 2 . Case 1(a). 



correct decimal orders, ttith the cTOrs being 2 % and 12%. The ptimarj- source of error is the 
restriction of respon.se representation to third order kernels only. This has been done to keep 
the computations to a manageable level. Inclusion of higher order kernels (5*., order) in 
the response representation can be expected to improve the accuracy of the estimates at 
increased computational effort. Another source of inaccuracy in the estimates is the finite 
length of samples and the ensemble size. The nondimensional time interval for sampling has 
been taken as 0.5 and 4096 samples are collected for an ensemble, for this numerical 
simulation which give a frequency bandwidth of ± I.O cycles / t and a frequency resolution 

of 0.488 X 10 cycles / r. Increased sample size and/or increased ensemble size can be 
expected to yield more accurate estimates. 

4.5.2 Case Studie.s 1 (b) and 1 (c) 

In the case study 1 (a), the values of the nonlinear parameters, were kept identical 

(= O.I), in response simulation. As the next cases, dissimilar values of these parameters are 
chosen for response simulation and their subsequent estimation from the simulated response. 

Case 1(b): =0.1; /l^^.= 1.0 

Ca.se 1(c): ^«^0.1; /l^j, = 0.01 

The remaining parameters remain the same as in Case 1(a). 

Ibc results obtained are shown in the figures listed in Table 4.1 


Table 4.1 Figures of Cases 1 (b) and 1 (c) 


Parameters 

Case 1 (b) 

Case 1 (c) 

^//f'>(£y), 

Figs. 4.13 (a)-(d) 

Figs. 4.17 (a)-(d) 


Figs. 4.14 (a)-(d) 

Figs. 4.18 (a)-(d) 


Figs. 4.15 (a), (b) 

Figs. 4.19 (a), (b) 

lA' 2^ 

‘^xxf '^yy 

Figs. 4.16 (a), (b) 

Figs. 4.20 (a),(b) 










First order kernel transfcMm (dB) First order kernel transform (dB) 
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0.01 0.1 

Nondimenslonal Frequency 

(a) First order direct kernel transform, {(d) . 
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Estimates of the first order direct and cross-kernel transforms: Case 1(b) 
(Contd.) 




Fif« order kerne) transform (dB) Frst order kernel transfonn (dB) 



0.01 0.1 1 
Nondimensional Frequency 


(c) 


First order cross- 


kernel transform. 



(d) First order direct kernel transform, . 


kernel transforms! Case 1(b) 


Fig. 4.13 


Estimates of the first order direct and cross- 



Ttlird ofdcr kernel fictw (4B) Third order kernel factor (dB) 
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(a) Estimate of ^ ^3 ^'^’^{(D,0),o}). 



Nondimcnsioaal Frequency 

(b) Estimate of a ) . 


Fig. 4.14 


Third order kernel factors: Case 1(b) (Contd.) 



Third order kernel factor (dB) (^,3^ 
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Nondimensional Frequency 


Fig. 4.1 6(a) Estimate of the nonlinear parameter, : Case 1(b). 



Fig. 4. 1 6(b) Estimate of the nonlinear parameter, : Case 1 (b). 
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Fkni CKcto* kmm^ trarmfCOTi (dB) First order kernel transform (dB) 








Hiirci wder kernel factor (dB) Third order kernel factor (dB) 
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0.01 0.1 
Nondimensional Frequency 



(a) Estimate of ^ (co, o), co) . 



(b) Estimate of 'Fj {(O, (o, co) . 


Fig. 4.18 Third order kernel factors: Case 1 (c) (Contd.) 
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0.01 

Nondimensional Frequency 


(c) Estimate of 'F 3 (o),Ci),a)) . 
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-25 


0.01 


Noadimcnsional Frequency 


(d) Estimate of ■^'^3 (o),<i},(o) . 


Fig. 4.18 Third order kernel factors: Case 1 (c). 
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Third order kernel transform (dB) Third order kernel transform (dB) 



(a) Estimate of ^ Hf\(D,CD,a)) . 



(b) Estimate of ^ (<», o), co) . 


Fig. 4.19 Third order measured kernel transforms; Case 1 (c). 
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Nondimensional Frequency 


Fig. 4.20(a) Estimate of the nonlinear parameter, : Case 1(c). 



Nonditnensional Frequency 


Fig. 4.20(b) Estimate of the nonlinear parameter, : Case 1(c). 
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The parameter estimates are given in Table 4.2 


Table 4.2 Parameter estimates in Cases 1 (b) and 1 (c) 


Parameters 

Case 1 (b) 

Case 1 (c) 

II 

1.109 

1.005 

4 ' (='‘ 5 .) 

0.499 

0.503 


0.016 

0.011 


0.089 

0.088 


1.18 

0.012 


In these cases too, the linear stiffness parameters are estimated accurately. Damping 
estimates involve higher errors, due to reasons cited earlier. The estimates of the nonlinear 
parameters show an error between 11-20%. 

4.5.3 Case Study 2 

For the next set of illustrations, the nonlinear parameters and A are kept fixed at 0.1, 

as in Case Study 1 (a), and the linear parameters are varied. In all the Cases 1, the direct 
linear stiffness parameters were chosen to be identical (= 1 .0) and the cross-coupled stiffness 
parameters were taken as half of the direct ones (= 0.5). 

In Case 2 (a) the following sets of values are now chosen for response simulation and their 
subsequent estimation from the response. 

At, = 1-0; At,. =2.0; = 

In Case 2 (b), the negative cross-coupling effect has been incorporated, with 
At, = = 1.0; At = At =-0.50 

The figures for these two cases are listed in Table 4.3. 
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Table 43 Figures of Cases 2 (a) and 2 (b) 


Parameters 

Case 2 (a) 

Case 2 (b) 


Figs. 4.21 (a)-(d) 

Figs. 4.25 (a)-(d) 





Figs. 4.22 (a)-(d) 

Figs. 4.26 (a)-(d) 


Figs. 4.23 (a), (b) 

Figs. 4.27 (a),(b) 

nN 

Figs. 4.24 (a), (b) 

Figs. 4.28(a), (b) 


The estimates are given in Table 4.4. 

Table 4.4 Parameter estimates in Cases 2 (a) and 2 (b) 


Parameters 

Case 2 (a) 

Case 2 (b) 

X^ 

1.01 

1.03 

x^ 

^}y 

2.03 

1.07 


0.992 

- 0.495 


0.013 

0.012 

x^ 

0.078 

0.108 

Xiyy 

0.107 

0.111 


The estimate of the direct linear term has been made by replacing term 1 on the r.h.s of 

equations (4.68) by itself It can be seen from Table 4.4 that, in Cases 2 (a) and (b), 

where the direct linear stiffness terms are dissimilar (unlike Case 1 (a)), the pattern and 
accuracy of estimates is similar to the previous cases. It can also be noted, in Case 2 (b), that 
negative cross- coupling effect has been correctly identified. 
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First order kemef transform (dB) First order kernel transform (dB) 



0.01 0.1 1 
Nondimenslonal Frequency 


(a) First order direct kernel transform, . 



(b) First order cross-kernel transform, . 

Fig. 4.2 1 Estimates of the first order direct and cross-kernel transforms; Case 2(a) 

(Contd.) 
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Frst tmmfcrm (dB) Plrst cH'der kem^ trarrsfomn (dB) 


16 



(c) First order cross-kernel transform, {co) . 



(d) First order direct kernel transform, . 


Fig, 4.21 Estimates of the first order direct and cross-kernel transforms; Case 2(a). 
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Tfmd OTder kernel factor (dB) Third order kernel factor (dB) 



(a) Estimate of ^3 (o), 0), co) . 



(b) Estimate of . 


Fig. 4.22 Third order kernel factors: Case 2(a) (Contd.) 
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(d) Estimate of 
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Fig. 4.24(a) Estimate of the nonlinear parameter, : Case 2(a). 
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Fig. 4.24(b) Estimate of the nonlinear parameter, /l^ : Case 2(a). 




Rrst order kernel transform (dB) PIrst Q^der kernel transform (dB) 
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0.01 0.1 1 
Nondimensional Frequency 

(a) First order direct kernel transform, (co) . 



(b) First order cross-kernel transform, 


Fig. 4.25 Estimates of the first order direct and cross-kernel transforms: Case 2(b) 
(Contd.) 
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TTiird order kernel factor (dB) Tliird order kernel factor (dB) 




Third order kerne! factor (dB) Third ^rder kernel factor (dB) 





Third order kernel transform (dB) Third order kernel transform (dB) 
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(b) Estimate of ^ Hf\a),a),co) . 


Fig. 4.27 Third order measured kernel transforms: Case 2(b). 
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lamlxte (dB) Estimated lambda (dB) 



Fig. 4.28(b) 




Apart from the linear and nonlinear stiffness parameters, other nondimensional parameters 
contained by the governing equation (4.3) are the damping ratios. All the previous case 
studies have been carried out for 1% damping 0.01). Numerical simulation is 

carried out to further, investigate the estimation procedure for a lower damping, with 
^xx~'yyy~ 0 00^ (keeping the values of the remaining parameters the same as those in 

Case 1 (a) ). The parameters and the corresponding figures are given below 



Figs. 4.29 (a)-(d) 

^ {(O, CO, CO), {CO, CO, CO) 

yx yyl-W >y XJ/ 2-U.l 

Figs. 4.30 (a)-(d) 

''H^^\ 60 ,co,o)), yH^^^{co,co,co) 

Figs. 4.3 1(a), (b) 


Figs. 4.32 (a), (b) 


The estimates are: 


It 

1.07 

fl 

0.488 

11 

0.002 


0.150 


0.170 


The estimates can be seen to be sensitive to damping and comparison with the results of Case 
1 (a) shows that, the results are less accurate in the present case, ■with lower damping. 
Averaging over increased sample and ensemble sizes can be expected to reduce the statisucal 
errors in the extracted kernels and improve the accuracy of the estimates. 
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(a) First order direct kernel transform, . 



Fig. 4.29 Estimates of the first order direct and cross-kernel transforms; Case 3 
(Contd.) 
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First order kernel traftslonn (dB) pj^gt kernel transform (dB) 


15 



(c) First order cross-kernel transform, {co) . 



(d) First order direct kernel transform, 


Fig. 4.29 Estimates of the first order direct and cross-kernel transforms; Case 3. 
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Third order kernel factor (dB) Third order kernel factor (dB) 




Third order kernel factor (dB) Third order kernel factor (dB) 



(c) Estimate of '^3 ^’^’^{ci), 0 ),{o) . 



(d) Estimate of ^ T 3 (o), (O, O)) . 


Fig. 4.30 Third order kernel factors: Case 3. 
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Third order kernel transform (dB) Third order kernel transform (dB) 



(a) Estimate of ^ H^\(D,0),co) . 



Fig. 4.3 1 Third order measured kernel transforms: Case 3. 
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Esttnated lan^da (dB) 



Fig. 4.32(a) Estimate of the nonlinear parameter, . Case 3. 





4.5.5 Influence of Ensemble Size 


It W3S observ ed thst incressing the ensemble size beyond 500 has insignificant influence on 
the first order estimates, the influence on the third order estimates is investigated for 
ensemble sizes 500, 1000, 1500 and 2000. The parameters are kept the same as in Case 1 (a) 
and the results are arranged, for clarity, as given below 

Ensemble Size 

500 1000 1500 2000 

Fig. 4.33(a) Fig. 4.34(a) Fig. 4.35(a) Fig. 4.36(a) 

Fig. 4.33(b) Fig. 4.34(b) Fig. 4.35(b) Fig. 4.36(b) 

^^3 " {(0,00. CD), 

^ H\^\co,(o,co) 


It can be seen from the above figures that the third order Wiener kernel transforms, 
^ H^\(0,co.(o) and ^ H^\co,a),0)) get refined with increasing number of samples in the 

ensemble. However, the refinement in the kernel factors, ^ {co.co.oo) etc., which are 

synthesised from first order kernel transforms, is marginal. The ensemble size was not 
increased beyond 2000, in the present study, due to data storage limitations of the computer. 

4.5.6 Influence of Measurement Noise 

The influence of measurement noise, which can be expected during an experiment, in the 
measurement of the excitation force and the response signals, is studied by contaminating the 
simulated force and response signals with 5% simulated random noise. The parameter 
estimation algorithm was found to be robust in the presence of this noise and the change in 
the estimates is marginal. For illustration, only the third order kernel factors, 

>^^'l~^'^'\co,(0,(o), yy^l~^'^’^{(j0,O},0)) and the third order Wiener kernel transform, 
yH^^\(O,0),a)), are shown in Figs. 4.37 (a) and 4.37 (b). The data corresponds to 
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0.01 0.1 


Nondimensional Frequency 

(a) Estimates of 6), fi?), ^'^l~^'^'^{co,0},co) and {a),co,co). 



0D1 0'' 


Nondimenslonal Frequency 


Estimates of ^ {(0,(0, (d),^ (o), O), a) and (o},0),co). 


Fig. 4.33 Influence of sample size on third order estunates: Case 1(a) 
(sarcple size = 500). 
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NondlmensionaJ Frequency 



0.01 0.1 
Nondimensional Frequency 


(b) Estimates of ^ (cD,co,a)),^ (o), o), (o) and (cd, o), co) . 

Fig. 4.34 InjQuence of sample size on third order estimates: Case 1(a) 

(sample size = 1000). 
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0.01 0.1 
Nondimensiona] Frequency 


(a) Estimates of ^ 'F 3 ^’^’^{a),CD,a)) and ^H 2 ^\o},co,(d). 



Estimates of ^ (o), o), o) , ^3 • 


Fig. 4.35 Influence of sample size on third order estimates; Case 1(a) 
(sample size = 1500). 
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Estimates of ^ ^3 ^^’^{co,co,Ci}) and ^H 2 ^\(D,co,ct)) . 



0.01 O"' 

Nondlmenslonal Frequency 


Estimates of and ^H^^\co,co,Ct)). 

Influence of sample size on third order estimates: Case 1(a) 

(sample size = 2000). 
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- Psi3yx with no noise 
-Psi3yx with 5% noise 
• Psi3yy with no noise 
-PsiSyy with 5% noise 
--H3y with no noise 
■H3y with 5% noise 


Nondimensional Frequency 

in the jfrequency range (0.0-0.03). 



Nondlmenslonal Frequency 


(b) {CO, (o,(o),yy {CO, CO, CO) , {co, co, co) 

in the frequency range (0.03-0.33). 


Fig. 4.37 Effect of measurement noise on third order estimates: Case 1(a). 
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Case 1 (a). The figures show the results for two different situations - (a) no measurement 
noise and (b) 5% random noise in both i. e. input as well as response. (The noise signals are 
uncorrelated random signals with rms values equal to 5% of the nns values of the excitation 
force and the response) . In Figures 4.37 (a) and 4.37 (b), the frequency range has been split 
into two. Fig. 4.37 (a) corresponds to the frequency zone, of interest for parameter 
estimation, 0.00-0.03 cycles / r , on a magnified scale. The remaining portions of the curves 
are shown in Fig. 4.37 (b ). 

4.6 Remarks 

The salient features of the parameter estimation procedure for a rotor-bearing system, with 
cross-coupling effects in the bearings, are the definitions of cross-kernels, which upon 
convolution with the excitation represent the response in terms of Volterra or Wiener series. 
For the nonlinear two-degree-of-freedom problem, Laplace transforms are employed to 
derive expressions for the first and higher order direct and cross-kernel transforms. The 
numerical illustration reveals satisfactory performance of the procedure for various 
nondimensional parameters of the governing equations. Such a study can be usefully 
employed for design of experiments. The rotor has been treated as a rigid body in this 
chapter. The case of flexible rotors is discussed next. 



CHAPTER 5 


PARAMETER ESTIMATION IN FLEXIBLE ROTORS 

The case of parameter estimation in flexible rotors is considered in this chapter. In the 
analysis so far, the rotor is treated as a rigid body and the shaft flexibility is not accounted 
for. The treatment in this chapter incorporates the influence of shaft flexibility. An 
attempt has been made to reduce the algebraic complexities by adopting the same rotor- 
bearing configuration as in Chapters 3 and 4. However, in order that there is no loss of 
the generality of approach, equations have been written in matrix form. The shaft 
flexibility has been included in the analysis through influence coefficients concepts. The 
shaft stiffness is taken to be purely linear with nonlinearity present in the bearing stiflBness 
terms. Bearing cross-coupling effects which have already been discussed in the Chapter 4 
has not been included in the analysis. Similarly, damping is treated to be linear with no 
cross-effects. The parameter estimation procedure follows previously described route of 
deriving expressions for the first and higher order direct and cross-kernels. These are to 
be extracted from the excitation force and response measurements. Third order kernel 
factors synthesised from measured first order kernel transforms are processed along with 
the measured third order kernel transforms for nonlinear parameter estimation. The 
procedure is illustrated through numerical simulation. 

5.1 Governing Equations And Response 


The equations of motion for a balanced rotor with a centrally located disc on a massless 
flexible shaft supported in bearings, shown in Fig. 5.1 are written as 
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(5.1) 




Figure 5.1 Flexible rotor in bearings. 
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In equation (5.1) Wj, are the system masses at stations 1, 2, and 3 respectively. 

The external white noise excitation provided at the stations 1 and 2 are /,(/), /aCO- 

^4 ’ unknown linear stiffness terms of the two bearings and and are 

their unknowTi nonlinear stiffness terms. The shaft stiffness parameters are represented by 
ky terms, which are defined as the resistive force at the ith station corresponding to a 

unit deflection at the y'th station, with all other deflections held to zero. The values of the 
ky terms can be obtained from the Strength of Materials formulae. 

linear damping terms .' 


Defining 

t- pt. 

iij =Cijl^m,p, 

and 


P = ylk33lm^; 

^ St ~ -^maxi ! ^iP ’ 

Xy = ky ! rriip'^ ■, i,y=l,2, 3 



^b, ~ ^max 


.N rrl 


3 „6. 


Im^p 


(5.2) 


for /-I, 2 

and considering only the direct damping terms, in order to keep the algebra simple 
equation (5.1) reduces to 
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The solution of the equation (5.3) is represented in terms of Volterra operators as 


''>7(r)='«f/.(r),/,(r)] = i 'W.[/,(r). A(r)] 


(5.4) 


with K denoting Xj , X 2 or X 3 . 



Further, 

'= 1.2 i = 1 , 2 ; 7 = 1,2 

+ 

;=l,2;y=l,2;t = l,2 

and the indiv idual operators are given in kernel form by 

OCi 

H\‘\f,{T)]= \'''h\'\r)l{r-v^)dt^ for / = 1,2 

- oC ' 

oc oc 

<«*OC “-OC 

for / = 1,2; j = 1,2 

00 00 00 

'^//^'■■'’*>[/,(r),/,(r),/ar)]= J j J''/jf^’*^(r,,r2,r3)/;(r- r^Z/r- r2)/i(r- ^3^ 


(5.6) 


-OC -00-00 


X dTy dtjdr^ for /=1,2; ; = 1,2; k-l,2 
For convenience, writing the Volterra operators as 


IC „U) 


r), /,( r). /*( r)] 

the response of equation (5.3) can be written as 

Xr)= 

'-1.2.;*=U '-1.2;y.l.2;*.l,2 


(5.7) 


(5.8) 


5.2 Synthesis of Higher Order Volterra Kernel Factors 


The Volterra operators are now determined on the same lines as in the previous chapters. 
The excitation forces /,(t), Mr) are replaced by c/,(r) and cf,(r) respectively, c 
being a constant. Noting equations (5.6), the resulting response of the system, becomes 
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i»l,2 


/«1.2;;=>1.2 


3 K AiJ,k), 

+ 1)1 + 
i=l,2;y=l,2;A:=l,2 


(5.9) 



Substituting equations (5.9) and derivatives in equations (5.3), one gets 


Jc 2;<;2 '■'n(W)- 


/= 1,2 


i=l,2;y=l,2 


+ ~'Vi 

/=1,2;>l.2;*=l,2 


f-l.2;j-l.2 i-1.2;2.1,2;*.l,2 


Sc’ 


Sc + 


1=1,2 


Sc’ + 

;=1,2;;=1,2 


3 x.{ij,kr 


i=l,2;y=l,2;/t=l,2 


+ 2 


0 

0 ^22 
0 0 


0 

0 

^33. 


Z^‘’’7r+ S';’ S<’’'‘’7p‘’’ 

_'=^’2 i=l,2;;=l,2 i=l,2;y=l,2;it=l,2 

Sc'=7,‘'’'+ S'^’ S'’’ 

.'=1.2 /=l,2;y=l,2 i=l,2;y=l,2;A=l,2 


3 Xi„iiJ,k)' 


+ 




^21 


Ml 


"^12 ^13 

^22 

A, 


23 


M2 


M3 


Sc'’>7f"+ S^’>F + Sc- -'72 

,i=1.2 /=l,2;j=l,2 /=l,2;>I,2;i-1.2 


Sc 7,‘''> + S'’’ '' 7?''’ + Sc’ '■ 73''’’*’ 

i=l2;y=l,2 /=l,2;y=l,2;Jt=1.2 

Sc'’ 7,"’+ Sc’ '-’75'-”+ Sc’”7f’'‘' 

/=1.2;7=1,2 i=l,2;j=l,2;k=l,2 

Sc'’7,''>+ Sc’'’7l'”+ Sc'’'’7^-‘’ 

1=1.2 /=1,2;7=1,2 (=I,2;y=l,2;t=l,2 




3^ 



r 

Sc '' 7 ,'"+ Sc’''‘ 75'-”+ Sc’''7f”‘> 




1*1,2 l•»l,2;y>*l,2 /=1.2;y=l,2;*=l,2 



0 

0 

1 

— 

3 


/i(^') 

0 4 

Sc '’ 7 ,“’+ Sc’'’7^'+ Sc’ '’7l'”'‘> 



/2(^) 

0 0 0 

,»1,2 /=1,2;7-1,2 i=1,2;;=1,2;*=I,2 

3 


0 

V 


Sc '’ 7 ,"'’+ Sc’'’75'-”+ Sc’ ''’7F’*’ 




i«i,2 i*I,2;7»1,2 /=1,2:7=1,2;/c=1,2 
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Summing up the responses of equal order, as follows 


I order response 

7i = 

E'7!'’ 



i=l,2 

II order response 

II 




i=I,2;y=l,2 

III order response 

II 



(5.11) 


i=\2j=U\k=\,l 

Noting the following symmetry of kernels (Schetzen, 1980) 


V2 - V2 » 




(5.12) 


and using equations (5.1 1)-(5.12), equations (5.10) are written in condensed form as, 




]c 7j +c" 772 + C 


.2 X2 

\C TJy +C" 772 + C 


3 X2 ” 

m\ 


^3„"^^2 ^3„" ^^3 ^3^^ 


\ + 2 


■^11 0 0 

0 ^22 0 

.0 0 ^33J 


[c ^'77! +C^ ^^112+0^ 
[c ^^Vi+c^ 


\\[c ■‘^71+^ ^Vi+c 
[c +c^ '‘^Vi+c 


^3]] 

% 


X^\+X^ 


^21 


X, 

X22 x^ x^ 


'23 


A31 .^32 


^33 


f[. 


^Ir, J.^2 Xt„ _l^3 ^1«J1 


[c '77, +C ’772 +C ’73 

|][c +c^ "^r]2+c^ ^'773 jy 

[c ""^77, +c^ ^^772 +c^ "^^773 


1.V 


0 0 
0 Aj; 0 

0 0 0 


[c"’77j +<7^ '^'t72+^" "‘^3 

3 f 


J Xj. 


fl 


j.^2 Jt2„ 4.^-' 


\C "^TJ^ -k-C" "^772 
[c ^377, +c^ ^'772 +C^ ^'773 


) = c 


fM 

/2W 

0 


(5.13) 


Equations (5.13) are power series in c with coefficients of c being '77^, ^77„ or 
^3 ri„ in = 1,2,3). The responses, 77„, ^2 n„ and ^3 r]„ are determined by equating the 
like powers of c as follows 


c terms: 
x\. 


^2 


^3 


1\ 

V\ 

l"\ 



^11 

0 

0 ' 

V ' 


+ 

1 

^^12 

/?.,3 



7i(r)' 

> + 2 

0 

^22 

0 ■ 



A 21 

^22 

^^3 


> = ^ 



0 

0 

^33 _ 

^ 377 ; 


A 31 

Xy2 

^^33 



0 


(5.14) 



terms: 












M ^ 

M 


^11 

0 

0 ■ 

""'72 


+ 

1 

All 




O' 

% 

^ + 2 

0 

^22 

0 

■^^ 7 ; 


X 21 

^22 + 

^3 ■ 


► = ^ 

0^ 

'72 


0 

0 

^33. 

V J 


A 31 

^32 

"^33 



0 

^ J 


(5.15) 


terms: 


r ^ 

'73 


■^11 

0 

0 ■ 

r 1 '' 

^'73 


'^11 +^6, 

Ai 2 

^i 3 

r 

""’73 

73 

> + 2 

0 

^22 ■ 

0 

-^^7; 

> + 


^22 + ^\i 

A22 ■ 

^273 > 

v J 


0 

0 

^33 _ 

''^73 


A31 

A22 

^33 




1 

0 

0 

1 

f ;ci „3 " 

7i 


'o' 

+ 

0 

0 

X2 3 

7i 

► = ^ 

0 > 


0 0 0 

J? 

V 


0 

< J 


Equations (5.14)-(5.16) can be solved sequentially. Taking Laplace transforms of 
equations (5.14), for zero initial conditions, one obtains 

[M]{7,W} = {F(s)) 
with 


[M]. 


S +2^jj5 4-/lu -f 


H2 


^21 


^ ^2 ^t>2 


1 

^3] 

1 ^2 


{7i(‘y)} = ^ 

"'rji (s) 

(s) 

>■ and {F(s)} = < 

Fiis) 

F2(s)> 


"^7, (s) 


0 

< ./ 


^\3 

^3 

+ 2^335 + /I 33 


The solution for rj^ (s ) , rjy (s) and (s) , from the above is, 

{rj,(s)}=lH]{ns)} 

where, the kernel transform matrix 

[H]=N-' 

whose individual elements can be worked out to be 


(5.17) 


(5.18) 


(5.19) 




^^Hl^\s) 
^3// (’>( 5 ) ^^Hl^\s) 


' ^ C *^) - ( '^ 23^32 +('^' + 2 s ^22 +^^22 + 25^33 +'^ 33 ))^-^ 

= (^13^32 +2 s ^23 +'^ 33 ))/ 

^'NI^^(s) = (-^13 (-5^ +25^22 +^22 +^^>+^ 12 ^ 23 )/^ 

^/fj \s) = (^23^31 ~^l(‘^^ +25^33 + A 33 ))/ D 

^ H\ ^( 5 ) = (—-^13^31 +(-S^ +25i^jI + Ajj + A^ ){s^ +25^33 ■^X^t,)')! D 

^ = (Aj 3A2 i —{s^ + 2s^ii +'^11 )'^3)/ 

^ ( j) = (“(■? + 25 (^22 "^22 )'^3 1 ■*■ ^1-^32 ) 

^( 5 ) = ( AJ2A3J — (5 + 25 ^ jI + Ajj + A ^ ) A 32 )/ Z ) 

= (- A , 2 A 2 , +( 5 ^ + 25 ^,, + A,J + A ^ + 25^22 +^2 +4 


D = (5^ + 25^1, + A„ + Ai, + 25^22 + >^2 + 4 X-"' + 25^33 + ^33 ) 

— ( 5 ^ + 25 ^jj + Ajj + A ^ )^ 23^32 ~ 2 - S '^22 ■*■ ^^2 ■*" )'^31 ( 5 - 21 ) 

~ ^I2‘^2l(^^ +25^33 ^33)'^ ^I2‘^23^3l ■^'^13'^1^32 


Taking Laplace transform of the equation (5.15), similarly give 


5^ 4- 25(^jj + Aji + A^ Aj 2 Aj 3 



0 

A 21 + 25^22 + "^2 + ^4 ^^23 


>. = ■< 

0. 

A 3 , A 32 5^ +25^33 +A33^ 

‘='72W 

V 


0 

J 


(5.22) 


which yield 
^2is) = 0 

’‘2 tj 2 (s) = 0 (5-23) 

^3 772 ( 5 ) = 0 

which shows that the second order kernel is identically zero i.e., 


/22 (r, , T 2 )='^'/j 2 (^ 1 » ^2 )='''^2 (^ 1 . ^2 ) = 0 


(5.24) 



In order to synthesise expressions for third order kernels, equations (5.16) can be written 

as 




■^11 

0 

0 ■ 

^ V ' 

'73 


^^ri '3 

> + 2 

0 

^11 

0 

< ^^T]'3 

> 4- 



0 

0 

^33 _ 




^21 
A3 1 


where the following abbreviations have been used 
gi(r) = -A^^'7j^(r) . 

q2ir) = -Al^^Vx{r) 


^\1 

2 i 3 



7i(r)' 

A 22 + 2 ^ 

^23 ■ 



72 (^) ^ 

^32 

3^33 

^^73 


0 


(5.25) 


(5.26) 


Equations (5.25) are linear in 73 terms, similar to equations (5.14) and therefore the 
solution in terms of Volterra operators as 

"2;73(r)='^//,[?„?2] (5,27) 

that is 

'^■f^3[/i(>-),/2(''))='^W,[«i,92] (5.28) 


However, since and are linear operators, one obtains 

r,,(t)=^'H,[q„q^\ 

'^'73(7')='^^i[9i-92] 

'3 ,3(r)='^W, [?„?,] 


(5.29) 



Noting the abbreviations (5.26), the terms on the right hand side of the above equations, 
are individually expanded as 


In the above equations 

">7, = 775*4^' ?j\^^ 

='=//<'>[/, (r)] 

'3,1 = '3,(l) + 3rj^(2) 

= '3//, «[/,(,)] + 


Laplace transforms of equations (5.30) gives 


+ 3 ^ 1 - 1 , 2, 2 ^ \ j / 1 - 2 , 2,2 ^ 


(5.32) 

where, , ^2 ,-S' 3 ) is the third order kernel factor which has been defined 


+53)^>//p>(5,) ">i/f*>(52) 
i = l,2;y=l,2; A:=l,2; /=1,2 (5.33) 

The third order kernel factors above, can be readily constructed, using equations from the 
first order kernels 

Similarly, the Laplace transforms of other terms on the right hand side of equations (5.29) 
can be worked out to be 

''/?!-’[ft(3tl.3t2.J3)] = +3 V3 >p 7-'’'.2 +3 ','!q32-1.2.2+V3tp7-W-2) 

'3 f/J ‘ ’ la ( J, , J, )] = SKJ-'-'’’ + 3 '3'' 'J'j-'-'-' + 3 *3'3 >P +'='! ) 

'3H; = ’ia(s,,S,,S3)] = Aj; ('3'.>P,^->.>’U3 +3 



^3//f')[(2i(5,,52,^3)] = ^3 orjiupM.U ^3*n|/l-l,2.2^^3:t, vj,l-2,2,2. 

3 3 ' 

''3//p^[^2(>yi>-y2>‘S’3)] = ^^ (""3^2 ^2- 1,1,1 ^3 x^x2^2~\X2 ^3^2i{/2-l,2,2_^jt3X2vj,2-Z2,2v 

3 3 

(5.34) 

In equations (5.34), the following third order kernel factors have been used in addition to 
those defined in equation (5.33). 

.2-., +.3 ''^<*>(.3) 

VI '¥<‘-J-*J\s„l2,S3) = - +S3 + J3)"'//<^'>(i,) '^<‘>(^2) ''.fi’<'>(i3) 

V2 'pV-‘-y-‘\s„S2,S3) = - +S2 +S3)‘^ffl'>(s,) ‘2//<*>(S3) 

/=l,2;y=l,2; A:=l,2; / = I,2 (5.35) 


The Laplace transforms of the third order kernels (equations (5.28)), can now be 
expressed as 

Bj (S, ,S2,S3) = A^ (*>^' 'Fj-*’*’* + 3 + 3 +^'^i v}/j-2.2.2 ) 

'■h ^ 2 ^1-^2 .Xj^2 vj/2 --1,2,2 ^jcjj:2 vj/2 “2,2,2 ^ 

//j (5, , 52 , 53 ) = 'Tj"’’’’’ + 3 + 3 ^2^' vpj-f 2,2 ^X2X, vj3j-2.2,2 ^ 

/’^2^2 ^2^2Vf;2-l,l,2 ^ 2^2 v}/2”’l,2,2_j^jr2X2 vp2~2,2,2x 

/>2 ^ 3 3 3 3 ^ 

7/3 (5i , 52 , 53 ) = (^'^' 'Fj"*’’’’ + 3 + 3 v}/j-2,2,2 ) 

^2 ^3^21|/2“M2 -^3^2Vj/2-l,2,2^jr3jr2 v|/2-2,2,2^ 

1^2 ^ 3 3 3 3 


( 5 . 36 ) 


5.3 Measuremeot of Wiener Kernels 

As stated in Chapters 3 and 4, measurement of individual V olterra kernels is not possible, 
while equivalent Wiener kernels can be extracted from the measured response if the 
excitation to the system is white and Gaussian. These Wiener kernels can then be used to 
generate the Volterra kernels. 



In the present case, the Wiener kernels of the nonlinear system are extracted by 
application of white Gaussian forces /,(r) and one at a time, i.e. first a white 

Gaussian force /, (r) with variance , is applied at station 1, while keeping the force 

at station 2, /2(r)= 0. The resulting responses ^>7, and are employed to 
extract the direct, Xj — x j -(at station 1 ), Wiener kernels and the cross ( x j - X2 and Xj - X3 ) 
kernels. In the next instance, a white Gaussian force ^2(^) with variance A 2 , is applied 
at station 2, while keeping the force, f\{T) = 0. The system responses, in this instance, 
are employed to extract the direct X2 -coordinate Wiener kernels and the cross (X2- 
Xj,X 2 -X 3) kernels. 

The system response, in terms of Wiener kernels, is now expressed, in the two individual 
cases as - 

/i (t) Gaussian white, with variance A^ and f^{T) = Q : 

(5.37) 

/2 (r) Gaussian white, \vith variance ^2 /jCf) = 0 : 

'5 7( r)= W[0, A ( r)]= ( r)] 

In equations (5.37)-(5.38), the Wiener operators are 



(with K denoting Xj, X2 or X3 for i = 1 or 2 ) 
and the individual operators are given in kernel form by 


(5.38) 


(5.39) 





! 

|-O0 

Cxi) 00 

"■«2'^[7;(0]= j J''>^2^(fi,r2)i;(r-r,)7;.(r-r2)cfri^/t-2 -4- J'"M'J’kr2,r2)rfr2 

-do -00 ^ 

( —00 

OC' 00 00 

'^f^3'^[./i(^)]= I J J'^^*^(^i.^2>^3)^(^-^‘i)yi(i’-r2)^(r-r3) c/r, rfr2i/r3 

—00 —00 — oo 
00 00 

-3A/ J \'''w^\rx,T2,T2)fi{r-Ti)dT^dt2 

-0X0 


(5.40) 

The relationship between the Volterra kernels of equations (5.6) and the Wiener kernels 
above, is the same as written in the previous chapter (for a third order response 
representation). These relationships are reproduced below. 

'^;2]'^(r,,r2,r3)='^wJ^(rj,r2,T3) 

*’/I0='^Wo+'^H'o(2) 


with 


00 

^ 1 ( 3) (^1 ) = ~ 34 j J *" (^1 > ^2 ’ ^ 2 ^2 

woo 

I 

00 

>^0(2) (^1 ) = -^i ' J'" 

-00 


(5.42) 


Employing the Fourier transform relations and the property of the Dirac-Delta function, 
5{o )) , given in the previous chapter by equations (4.46) and (4.47), the responses (5.37) 

or (5.38) can be expressed as 



’^?7(r)=‘fVo + 


{0)2 -0)2 )do )2 


oc 00 00 




— OC —00 —00 


3 v 4 , J 


+ K being x, , X2 or X3 and i = 1,2 

(5.43) 

A complex exponential filter, similar to the one described in previous chapters, has been 
used for measurement of the individual Wiener kernel transforms. Referring to Fig. 5.2, 
the output of the filter is 


—00 


(5.44) 


and the ensemble averages of the outputs of the circuit, for k equal to Xj ,X 2 and X 3 are. 




jcor . cK ivO'h 






-aJf; r (fflj ,-®2 )d0. 


+ J \''Wj‘\0^,a)2,CD^){Fii(o^)F,ico2)Fii(o^)F’{co) 


— 00 — CC —X 


^ dco^dco^dco^ 


-3A,'^\’^\’‘W^'\cd^,co2-co2)iFi{(o^)F*{(^^ ^^dco^dco^ 


(5.45) 



Mr)^F,i(D) 


Fig. 5.2 



Scheme for evaluation of the first order direct and cross Wiener kernel 
transforms, 


146 






Since / (r) is stationary Gaussian white noise with zero mean and variance A, , the 
Fourier transform, is also a stationary Gaussian white noise process and 

employing the properties (ref. equation 3.30) of the ensemble averages of these products 
transformed functions, equations (5.45) get reduced, after some algebra, to 

Due to the equivalence of time and ensemble averages, the ensemble average 

r]{T)z^'\T'^ can also be written as, 

('^7(r)z^'^(r))= lim - J Tj{T)z^'\T)dT 

\ / Y-^00 'J' 12^ 

= F*{co)’'ti{o}) (5.47) 

Equation (5.46) and (5.47) give 

= F*(o)) (5.48) 

from which the expression for the first order Wiener kernel transform is obtained as 

= (5.49) 


Since k takes values Xi,X 2 ,x^ and i takes values 1 and 2 the direct kernel transforms 
^'lF^*^(u)), and the cross kernel transforms 

* 3 pfO)(^) extracted from the measured responses ^’;7(cy), 

fj(ct)), T](o}) and the applied force (F,(ty) with variance A^ or F 2 (co) with 
variance A ^ ), through equation (5.49). 


For measurement of the third order kernel transform, a circuit involving three exponential 
delay filters as shown in Fig. 5.3 is considered. The output, z^ \t), of the exponential 


filters is 


rW(r) = 


-OO 






(5.50) 














The ensemble averages of the outputs of the circuit (Fig. 5.3) are 






+ ( 0)3 )S(-Ci}i - 0)2 )e ] 


+ 


6Af '»'«(ffl,,<i.2,a)3) 
(5.51) 


The equivalence of time and ensemble averages gives 

1 th 


'' T]{t)z^‘\t)) = lim- I ’'T}{T)z^'\T)dTi 
/ T-^ooT-tI2 


— F,(^ ^(^1 +< 2^2 ”^^ 3 ) 


(5.52) 


Equations (5.51) and (5.52) give the expression for the measurement of the third order 
Wiener kernel transform as 


( 0 )j , 0 ) 2 ,£i> 3 ) = — ^ k * (<^1 )^* (^ 2 )F*( 0 >Z ) 7(<yi + 6 ) 2 + 0 ) 3 ) 

■J y jJ III J 


6A 


6 A. 


^H/(0 


W^^•>{o}^)Sio)^ +o)^)+''wl‘\o)2)S{a)^ + 03 ) 


i \_+’^W^'\o}2)S(co^ +CO 2 ) 


(5.53) 


Measurements are made for special trispectral kernel transforms with 
6>] = 0)2 - ( 1 ) 3 = CO . As stated in earlier chapters, these transforms being functions of 
only one variable <2), are easier to compute and interpret. For such trispectral kernel 
transforms the expressions (5.53) reduce to 


■*1 w(0 


W^‘>(o},o>,a)) 


1 


6A: 


^ 2 «/(') 


Ify \( 0 ,( 0 ^co) 


6Af 


f=;*(ffl)f'''7(3<B) 
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6A. 

i 

lF/'\tu)+^' 1F/''(6))+^‘ wI'\co) • 


1 

6 A. 

i 

wl‘\o))+'‘^ w}'\co)+"^- w}‘\co) 
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■"^3 11^(0 


\(o,co,(o) 


6Af 

1 




*77(36)) 


6A. 


(6))+"'' wl‘\ico)+^^ (CO) 


S(co) 

S(Q)) 


(5.54) 



5.4 Parameter Estimation 


The first order direct and cross Wiener kernel transforms and the third order 

direct and cross Special Trispectral Wiener kernel transforms (co,a,Q>) are 
estimated using equations (5.49) and (5.54) respectively, from the measurements of 
spectral component F{q}) of the excitation force, its variance A, and the spectral 

components of the corresjxinding response. Subsequently, for a third 

order representation of the system response, noting the equivalence between the Volterra 
and Wiener kernels (equations 5.41 and 5.42), the first order Volterra kernel transforms 
and the third order Special Trispectral Volterra kernel transforms can be computed as 

(6), ry, co)=’^WP {co, (o, co) 

where 

1 -a> )da (5.56) 

^ ' ~00 

The linear parameters can be obtained from the above estimates of the first order Volterra 
kernels, //{’’(6i), Noting 

the algebraic expressions (eqns.5.20, 5.21) of these kernel transforms, the linear stiffiiess 
parameters of the bearings ,X^ are estimated through a complex curve fitting routine 
(Levy, 1959). In addition to the bearing stiffnesses, the damping ratios ^,,, ^ 22 » ^33 
and the mass ratios //j , Hi are also obtained from the curve fit routine. 

The nonlinear parameters are computed from the estimates of the Special Trispectral 
Volterra kernels , 

. Noting that these kernel transforms are 
estimated by application of a single white Gaussian force at a time, the synthesised 



expressions (5.36), for the third order Volterra kernel transforms, also get reduced, in the 
two individual cases, to 

/, (r) Gaussian white, with variance and /^(r) = 0 : 


'''H^{co,a},6)) = A^^ (^1^14/j-U 

’’) + >lJ^ (^>""2 4/2-1.1,1) 


^27/3 (di,n;,<y) = /l^'(^2^>'Fj"’’ 

X + X>2 (^'^''F3^~'’'’‘) 

(5.57) 

^37/3 ia),co,Q}) = A^^ (^3^1 4/ wj 

+ (^3^2 4/2-1,1.1) 



/2 (r) Gaussian white, with variance A^ and (r) = 0 : 

(6),<a,4y) = vj/j-2.2,2 ^ ^ vj/2-2,2.2^ 

"^2 7/3 (cD, 0 ), 0 )) = A^ vj/j-2,2,2 ^ ^ ^X2X2 Xf, 2-2.12 ^ 

^2 7/3 (fij, 6), 6>) = ) + A^ Xf> 2-12.2 ^ 

It is to be noted here, that application of two forces (individually) is required for 
estimation of the linear cross-coupling terms. However, since no cross-coupling has been 
included in the nonlinear terms in the present analysis, only two equations out of the set 
of six expressions (5.57) and (5.58) are sufficient for the estimation of the two nonlinear 

unknown parameters A^,A^ (cross-coupled nonlinear parameters being 

taken zero). Estimation of A^,A^ has been carried out here, using the first two equations 
from the set (5.57). For these two equations, using equations (5.55), we have 

^ [ fc * (n;)f tj(3co)] ( o))Sico)] A^ (^■"> ) 

6^,H J 2^iL -1 ^ 

+ A^ (V2 4/ 2-1,1,!^ 

(5.59) 

-Xr {F*(£0)f "''7(36)) =/lJ; 

J 2^iL J 

+ < (" 2 - 24 / 2 -i.U) 

The above equations are solved simultaneously for the nonlinear parameters A^ and A^^ . 



5.5 Computer Simulation 

The procedure is illustrated through numerical simulation of the response for the 
nondimensional equations (5.3). The excitation forces are simulated through random 
number generating subroutines and are normalised with respect to the maximum value of 
f\ (0 • ^ typical sample of the excitation is shown in Fig. 5.4 (a). The power spectrum of 
the input averaged over an ensemble of 2000 force samples is shown in Fig. 5.4 (b) 

The shaft stiffness matrix for the simply supported system, carrying a centrally located 
disc as (refer Fig. 5.1, for station numbering) is computed using Strength of Materials 
Formulae (Childs, 1990) 


( 5 . 60 ) 


^11 

^12 

^13 

_ 12£:/ 

■ 1 

0 

-i 

^21 

^22 

^23 

0 

1 

-1 

_^31 

^32 

^33 _ 

-1 

-1 

2 _ 


which gives the nondimensionalised matrix as 


^12 ^13 



0 

-( 0 . 5 ///,)' 

^21 ^ 22 ^23 

= 

0 

( 0.5 + 4)///, 

-(O.5///2) 

/I31 A32 ^33 


- 0.5 

- 0.5 

1 
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where /i, ,j^2 ^ ^2 ^^3 respectively. In the rotor 

configuration shown in Fig.5.1, m, and m 2 are the masses effectively seen by the 
sensors at the bearing ends and will be small in magnitude in comparison to m 3 . 
However, for simplicity in numerical simulation the mass ratios mj / and m 2 / m 3 are 
chosen to be each equal to 1.0. The values, of the individual elements of the above 
lambda matrix, therefore become 


111 _ ^2 - 0-5 

Aj3 = A 23 = 


/I33 =1.0 


b2 


A21 


^31 


A32 = - 0.5 


( 5 . 62 ) 


In order to illustrate the numerical results when shaft and bearing stiffness are both 
equally significant numerically, the values of the linear bearing stiffness parameters are 


taken as 
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Fig. 5.4(a) 
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Fig. 5.4(b) 



For the above set of values of the linear parameters the estimation procedure is illustrated 
for the following sets of the nonlinear parameters - 

Casel: ^^=0.1 =0.1 

Case 2; =0.01 

Damping is chosen as =^22 =^33 = 0-01 in both the cases. The same set of 

excitation as shown in Figs 5.4 (a), (b) is employed for the two cases. The governing 
equations are then numerically solved through a standard fourth-order Runge-Kutta 

subroutine, to obtain the responses 17 , ^ 17 and "3 . These responses are fed as inputs 
to the parameter estimation algorithm. The various first and higher order kernels are 
extracted from the responses and consequently parameter estimation is carried out. The 

output consists of the linear stiffriess parameters the damping ratios 

'all ’ “3 22 ’ * 333 ’ mass ratios //j,/i 2 and the nonlinear stiffness parameters 

The estimated parameters are compared with those originally used for the simulation of 
response. 

Casel: 

For the first set of values A^ =0.10, A^ = 0. 1 0 , the nondimensional responses rj , 

■*^2 Tj and ^3 fj numerically resulting from the force of Figs. 5.4 (a), (b), applied at 
station 1 , are shown in Figs. 5.5 (a), (b) and (c). (No force is applied at station 2). Their 
corresponding power-spectra are shown in Figs.5.6 (a), (b) and (c). The critical 
frequencies of the system can be noticed from these figures to exist at 0.097, 0.175 and 
0.219 (cycles/ r). These frequencies correspond to those obtained from the eigenvalue 
solution of the linear nondimensional stiffriess matrix described in equations 5.61-5.63 
and can be called It can also be noticed from the power-spectra of Figs. 5.6 

(a),(b),(c) that the system nonlinearity is not apparently, equally displayed by responses 
J 7 , and The plot of "*2 77 (Fig.5.6b) displays additional peaks at frequencies 
0.263, 0.291, 0.341, 0.403, 0.491, 0.525 and 0.569 cycles/r, which can be identified as 
{2(0^ -CO 2 ), (3a>j), ( 2 ft )3 ( 20 ) 1 + 0 ) 3 ), ( 0 ), + 0)2 + 0 ) 3 ), ( 30 ) 2 ) and ( 20 ) 2 + 0 ) 3 ) 












harmonics respectively. These harmonics are not distinctly, visible in the plots of tj 

(Fig. 5.6a) and ^3 tj (Fig. 5.6c). This is due to the fact that stations 1 and 3 being closer 
the station of force application (station 1), the linear response levels, are higher at these 
stations in comparison to that at station 2 and overlap the nonlinear contributions. Similar 
response plots can be obtained by application of the force at station 2, while keeping the 
force at station 1 equal to zero. 

The applied force and -resultant response at the three stations are employed in equations 

(5.49) and (5.55) to extract the first order Volterra kernel transforms, H^\a>) shown in 

Figs.5.7(a-f). These kernels show the three critical frequencies of the system, mentioned 
above. The errors in these estimated kernels are given in Figs. 5.8 (a)-(f). The errors are 
computed by comparing the estimated kernel transforms with those obtained from the 
exact analytical expressions (5.20). The errors are high in the vicinity of the frmdamental 
frequencies, as in the cases of rotor configurations, discussed in the earlier chapters. 
However, the errors are within 10% zone upto a nondimensional frequency of 0.07 
cycles/ r . Parameter estimation is therefore carried out in this frequency range. 

The third order kernel factors , 

and 

^ 3-*2 estimated from the first order kernel transforms, are shown in 

Figs.5.9 (a)-(f). However, in the present case, the last two factors, *’’’*(u>,tu,to) 

and ^<3^2 are not required since the nonlinear parameter estimation can 

be carried out from the first two equations only, from the set of equations (5.57). 

The third order measured kernel transforms ^'H^\co,0),co), and 

(ct),CL>,(0^ , extracted from the measurements of the force and response, in 

accordance with equations (5.54), are shown in Figs. 5.10 (a)-(c). As before, while the 
first order kernel transfonns are estimated in the entire available frequency range 0.0 - 
1 .0, the third order kernel transforms, involving a 3o) factor, have to be restricted to one- 
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Normalised error in the first order estunates: Case 1 (Contd.) 
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third of this frequency zone (i.e. 0.0 - 0.33). It can be observed, from the figures, that 
while the measured third order kernel transforms ate reasonahly accurate in showing the 
harmonic at hti 2,3 /3 (at nondimensional frequencies = 0.032, 0.058 and 0.073 

cycles/ r ), the identification of the harmonic at r», „ ( at nondimensional frequencies = 
0.097, 0.175, 0.219 cycles/r) is weak, due to higher statistical errors, mentioned earlier. 
The estimation of nonlinear parameters and , from these kernel transforms, is 
therefore restricted to the frequency zone of 0.0 - 0.07. 

The estimates of the nonlinear parameter and , obtained in accordance with the 

relationships (5.59), are shown in Figs. 5.11 (a)-(b). A fourth order polynomial curve 
regressed through the estimates of these nonlinear parameters , over the frequency range, 
is also shown in these figures. The mean values of the estimates are found to be 

A^ =0.09 and A^ =0.12 

6 , l>2 

(The exact values of the above nonlinear parameters are those chosen for response 
simulation, that is, A^ = A^ = 0 . 10 ). 

The linear parameters have been estimated from the first order kernel transforms of 
Figs.5.7 (a)-(f), using a complex curve fitting routine (Levy, 1959). The estimates are 
/^j =0.956 // 2 = 0.917 

^^=0.487 =0.918 

^11 = ^22 ~^33 - 0 01 1 


Case 2 

In the case study 1, the values of the nonlinear parameters, A ^ , A^ were kept identical 

(=0.10), in response simulation. As the next case, dissimilar values A^ =0.1, A^ =0.01, 

of these parameters, are chosen for response simulation and their subsequent estimation 
from the simulated response. The linear parameters are the same as in Case 1 and the 
forcing function is the same as shown in Fig. 5.4 (a),(b). The results obtained are shown 
in the figures listed below. 



Bstimated I 






nondimensional responses ^' 7 , 
power spectra of ■** 77 , jj 

first order Volterra kernel transforms, 
third order kernel factors 

.,.x, 

third order measured kernel transforms 

H^^\a),0,o)), Hf\a),(ji),(D) 

nonlinear parameters 

The estimates computed in this case are 
;ii!' =0.106 a!![ =0.04 

&l t>2 

and 

//j = 0.99 //j = 0.95 

;i^ = 0.49 = 0.95 

^ 11 =^ 22 =^ 33 = 0-013 

Reasonably good estimates can be seen to be obtained, in both the numerical cases 
illustrated above. The accuracy lies within the 20% zone for the nonlinear parameters, 

except , in the second case. However, the order of the magmtude of the estimate is the 

same as that of its exact value. It is to be noted that the response representation has been 
restricted to third order kernels only, in this study. Also, the ensemble size has been 
restricted to 2000. Increased sample size and/or increased ensemble size can be expected 
to yield more accurate estimates. 
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Power-spectra of the response: Case 2 
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Estimates of the first order kernel transforms: Case 2 (Contd.) 
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Estimates of the nonlinear parameters: Case 2. 




5.6 Remarks 


The pflrameter estimation procedure, developed in the previous chapters for parameter 
estimation in rigid rotors, has been extended to incorporate shaft flexibility, in this 
chapter. This has served to illustrate, the general nature of the procedures adopted for 
nonlinear parameter estimation. Relatively few numerical exercises have been carried 
out, in this chapter. The accuracies of the estimates with various other nondimensional 
parameters and measurement noise can be expected to follow trends similar to those in 
the earlier chapters. 



CHAPTER 6 


EXPERIMENTAL INVESTIGATIONS 

Experimental investigations have been earned out, in addition to the computer simulation 
described in earlier chapters, for further illustration of the parameter estimation 
procedures. The experimental studies have been carried out on an existing laboratory 
rotor-rig. These studies are restricted, due to laboratory constraints, to the case of a rigid 
rotor supported in ball bearings. For such a bearing, cross-coupling stiffiiess parameters 
are negligible in comparison to the direct stiffness coefficients. The experimental rotor- 
bearing system, if the shaft is rigid, corresponds to the single-degree-freedom system 
considered in Chapter 3. The rotor-bearing system is set into vibrations by providing 
white noise excitation at one of the bearing caps. Measurements for excitation force and 
the resultant vibration response are made simultaneously at the point of force application 
itself, i.e. at the bearing cap. The measured data is then processed in accordance with the 
procedure described in Chapter 3, for linear and nonlinear parameter estimation. The 
experimentation is carried out for various force excitation intensity levels. Variation in 
the excitation force levels provides with different values of the non-dimensional 
nonlinear parameter, i . The rig and the instrumentation are shown schematically in 
Figure 6.1. A detailed description of the rig, instrumentation and results follows. 

6.1 The Rig 

The laboratory rig consists of a disc centrally mounted on a shaft supported in two 
identical ball bearings. The shaft can be driven, if required, through a flexible coupling 
by a motor. The disc has a mass of 0.815 Kg. The shaft is 0.16 m long with a diameter of 
10 mm. The bearings are supported in pedestals, comprising of an upper cap and a lower 
case. The pedestal is mounted on a steel base, which in turn is bolted to a foundation. 

The bearings are SKF make with the following specifications - 




Fig. 6. 1 Schematic diagram of the rotor rig. 








Ball bearing Type SKF 6200 

Number of balls 6 

Ball diameter 6 mm 

Bore diameter 10 mm 

Outer diameter 30 mm 

Pitch diameter 20 mm 


Inner ring ball race radius 3.09 mm 
Outer ring ball race radius 3.09 mm 
Allowable pre-load 0 - 2 // m 

Rotor mass per bearing 0.41 kg 

Figure 6.2 shows the overall experimental arrangement. The rotor-bearing is excited 
through the bearing cap by a electro-dynamic shaker. The shaker is mounted on a 
specially designed rigid frame. Two-plane rigid-rotor balancing is carried out before 
proceeding with the experiments. 

6.2 Instrumentation 

White noise excitation is generated by a random signal generator (Bruel& Kjaer make, 
Type 1405). This noise generator has a frequency range of 20Hz to 100 kHz and can 

provide a maximum uniform spectral density of 10~^ V^/Hz upto 50 kHz frequency. The 
white noise signal is amplified through a power amplifier and fed to the shaker. The 
Power Amplifier is Briiel & Kjaer make, Type 2706 with a frequency range of 10 Hz to 
20 kHz, with a signal distortion limit of approximately 0.2%. The electrodynamic shaker 
is also of Briiel & Kjaer make, with a frequency range of 0 - 18 kHz and a peak-to-peak 
force rating ranging upto 10 N and a maximum displacement of 6 mm. The excitation 
force and response are measured by an impedance head attached between the shaker and 
the bearing cap. The close up of the rig including the shaker and the impedance head is 
shown in Figure 6.3. The impedance head is Briiel & Kjaer make. Type 8001 with 

accelerometer sensitivity of 3 pC/ms force gauge sensitivity of 370 pC/N and 1- 




Close-Up of the 


shaker and impedance head 




10,000 Hz frequency range. The signals from the impedance head is fed to a computer 
through standard Briiel & Kjaer, Type 2635, charge amplifiers. 

6.3 Data Acquisition 

The analog signals from impedance head through charge amplifier are digitized in the 
computer through a 32 bit Analog to Digital conversion card. The data is acquired into 
the computer through a virtual instrumentation (VI) program using LabVIEW software. 
The data acquisition c^d (National Instruments, Texas; Type-AT-MIO-16E-10) has a 12 
bit resolution and a maximum sampling rate of 100 kS/s. The front panel and back panel 
of the VI program developed in LabVIEW for the data acquisition are shown in Figures 
6.4(a) and 6.4(b). The data for the excitation force and the vibration response is acquired 
at 4096 number of equispaced time (/) instants, for a total period of 0.32 sec. Such a set 
of data is treated as a set of force and response samples. A total of 2000 such samples are 
collected for both the force excitation and the acceleration response and treated as an 
ensemble of force and response samples. Statistical averaging is carried out over this 
ensemble for further processing. 

6.4 Case Studies 

As mentioned earlier, experimental studies are carried out for various force excitation 
intensity levels. These force intensity levels are set by choosing the spectral density level 
on the white noise generator. The maximum force, , obtained in a force sample, is 

incorporated in the definition of the non-dimensional nonlinear parameter A (ref. 
equation 3.3). For the same rotor-bearing configuration, different values of provide 
different values of X and the results obtained for the estimates can be compared in a 
statistical context, as in Chapter 3. 

Case 1: 

As a first case the spectral density level on the white noise generator is set such that it 
gives = 3.27 N, as measured by the impedance head on the bearing cap. A typical 
sample of the wiiite noise excitation force is shown in Fig. 6.5 (a). The response in terms 
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Fig. 6.4(a) Front Panel of the LabVIEW VI for data acquisition 






Fig. 6.4(b) Back Panel of the LabVIEW VI for data acquisition. 
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Fig. 6.5(a) Typical sample of the input force: Casel 
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Fig. 6.5(b) Typical sample of the response: Casel 


of displacement, is obtained throngh integration of the acceleration signal. The typical 
sample of the displacement is sho™ in Fig. 6.5 (b). In order to check the noise level in 
the .nea.surements, all the instroments, including the eleetromagnetic shaker, were put on. 
The gain of the power amplifier between the signal generator and the shaker was adjusted 
at the lowest level so as to have no excitation conditions. Typical signals recorded, from 
the impedance head, under this condition, are shown in Figs. 6.6 (a) and (b). Comparison 
of these with the force and response signals of Figs. 6.5 (a) and (b) respectively, shows 
that the noise to signal ratio, in the extreme cases are 2.8% and 3.5% respectively. The 
influence of measurement noise on the parameter estimation procedure was discussed in 
Chapter 3 and the procedure was illustrated to be robust in the presence of 5% noise. The 
present noise levels are rather small and it can be safely assumed that their influence on 
the estimates will be insignificant. The power spectrum density of the excitation force is 
obtained by averaging over the ensemble formed from 2000 individual samples force 
signals. The power spectrum density for this case can be seen to be equal to 


0.36//'^ ///z (-4.436dB) from Fig. 6.7 (a). The power spectrum density of the 

displacement response, similarlj averaged over the ensemble is shown in Fig. 6.7 (b). 
The frequency in these two plots ranges from 0 to 6400 Hz. The two peaks are observed 
in Fig. 6.7 (b), a prominent one at a frequency equal to 1125 Hz and a minor one at a 
frequency of 3375Hz. The predominant frequency of 1125 Hz, represents the 
fundamental frequency of the system. (This was confirmed separately also, through 
repeated rap-tests, using an impact-hammer.) The minor peak at 3375 Hz suggests the 
presence of a nonlinearity in the system which can be approximated by a cubic function. 
l*hat the system nonlinearity is of cubic type was further confirmed by applying a 
harmonic force on the system, through the Signal Generator, instead of a random white 
noise. A sweep test with harmonic excitation was performed and the response was noted. 
It was observed that along with a prominent peak corresponding to the harmonic 
excitation frequency, the plot of the response contains a minor peak at the 3 multiple of 
the harmonic excitation frequency. This observation is displayed in the series of Figures 

6.8 (a)-(f)- 
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Fig. 6.7(a) Power spectrum of the input force: Casel 

(averaged over 2000 samples). 
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Fig. 6.7(b) Power spectrum of the response: Casel 
(averaged over 2000 samples). 
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Fig. 6.H(a) Harmonic excitation at ty = 1000 Hz. 
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I'ig. 6.8«c! Hamionic excitation at 0) = 2000 Hz. 
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Fie. 6.8(d) Resmmsc to 
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Fig. 6.8(c) Harmunic excitation at ty = 3000 Hz. 
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g. 6.8(f) RcsjHnise to the harmonic excitation of (e) above. 




The white neise excitation and the response measurements are employed to construet dte 
first order Volterra kernel transform thmugh use of equations (3.34) and (3.42> 

(3.4.3), Ihe first order Volten-a kernel transfonn thus obtained is shown in Fig. 

6.9. I he algebraic form of H^{co) from equation (3.12) is 

HiUo) , ^ 

4 Jc(oi k 

I/w ■ (6-1) 

{m; (o")^ j2^(0(o„ 

1 he linear par. imetei.s and ^ estimated from //|(ry). The natural frequency has 

been readily identified a.s 

a)„ 1125 11/ 706K,58 rad/.s 

I'he estimate of damping, using a least square fit of the curve of Fig. 6.9 with equation 
(6. 1 ), i.s found to be 
I 0.008. 

Fhc third order kernel factor is synthesised from the first order Volterra 

kerne! tr.ue.fnni of Fig. 6.9, in accordance with equation (3.20). The absolute 

value.s of 4'j(ft^^j,wr) and are shown in Fig. 6.10 (a), (b). As in the case nf 

the eonif'utei simulatinn of ('haptcr 3, it can be observed, from the Figs. 6.10 (a),(b) that 
while tfie measuied third order kernel transform is reasonably accurate in 

showing the harmonic at / 3 (=375 Hz.) , the identification of the harmonic at 
m„ 1 1 25 Hz) is weak. Improvement in these plots can be expected if the number of 
samples in the ensemble, over which averaging is earned out, is increased. In the present 
study, however, due to data storage limitations in the computer, the number of samples in 
the ensemble has been restricted to 2000. The estimation of A, from these kernel 
tramsfurms, is therefore restricted to the frequency zone of 0-1000 Hz. The nonlinear 

stiffness A is obtained from the plots of and hy employing 

equation (.3.44), The estimate of ^ ^ , thus obtained is plotted in Fig. 6.1 L A fourth order 
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polynomial curve regressed through the estimates of over the frequency range, is 
also shown in Fig. 6.1 1. The mean value of the estimates of is found to be 2.99 x 
10 N/m . Assuming that both the bearings are identical and act parallely, the 

nonlinear stiffness of each bearing can be written as 1.49 x 10^^ N/m^ . The 
corresponding value of the nondimensional nonlinear parameter X is calculated as 
0.366. 

The sign of can be readily identified as negative, from the discussion of Section 
3.6.6, and by observing the signs of /3 peaks of the real and imaginary components 

of 'F 3 (cy,<y,fu) and H^{ci),0),a)) shown in Figs. 6.12 and 6.13. 

Case 2: 

As the next case, the excitation force level is lowered down, such diat = 1.1 N. 

Typical samples of the excitation force and the displacement response are shown in 
Figs.6.14 (a),(b) respectively. The power spectrum density of the excitation force (Fig. 

6.15(a)) is 0.038 ! Hz (-14.2dB), in this case. The power spectrum of the 

displacement response, shown in Fig.6. 15(b), as in Case 1, shows a major peak at 1125 
Hz and a minor one at 3375 Hz. The first order Volterra kernel transform, Hi(co) is 

shown in Fig.6. 1 6. The natural frequency identified from the plot of {of) is 
co„ = 1125 Hz = 7068.58 rad/s 
while damping has been estimated as 
^ = 0.0078. 

The third order kernel factor, ^ 2 ,^( 0 , a),(o) and the third order Volterra kernel transform 
H 2 {o), 0 ),co) have been plotted in Fig. 6.17 (a) and 6.17 (b), respectively. The estimate 

of the nonlinear stiffness, , is shown in Fig. 6.18. The mean value of k for each 

bearing, in this case, has been found to be 3.58 x 10^^ N / . It can be observed, 
through a comparison between the plots of Cases 1 and 2, that while the identification of 
the /3 peak is sharper in H 2 {o), 0 ),(o) if the excitation force level is lower (Case 2), 
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Fig. 6.14(a) Typical sample of the input force: Case2 
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Fig. 6. 1 4(b) Typical sample of the response ; Case2 





tm 

I 

(/> 


Cu 

6 

< 



Frequency (Hz) 


Fig. 6. 1 5(a) Power spectrum of the input force: Case 2 



Fig. 6. 1 5(b) Power spectrum of the response: Case 2 
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Fig. 6.18 



the estimate of the nonlinear parameter, , shows better uniformity over the frequency 
range for a higher excitation force level (Case 1). The corresponding value of the 
nondimensional nonlinear parameter A, is calculated in this case as 0.1 14. 

Case 3: 

As a third case, the excitation force level is lowered down still further, such that = 

0.288 N. Figs.6. 19 (a),(b) respectively shows typical samples of the excitation force and 
the displacement response. The power spectrum density of the excitation force 

(Fig. 6.20(a)) is 0.003 / Hz (-25.68dB), in this case. The power spectrum of the 

displacement response is shown in Fig.6.20 (b), while the first order Volterra kernel 
transform, (fi?) is shown in Fig.6.21 . The natural frequency identified from the plot of 

Hi((d) is 

= 1 1 28 Hz = 7068.58 rad/s 
while damping has been estimated as 
^ = 0.0079. 

The third order kernel factor, 'y^ia>,a>,a)) and the third order Volterra kernel transform 
H^{co,co,co) have been plotted in Fig. 6.22 (a) and 6.22 (b), respectively. The estimate 

of the nonlinear stiffiiess, is shown in Fig. 6.23. The mean value of k for each 

bearing, in this case, has been found to be 9.50 x 10^° The corresponding value 

of the nondimensional nonlinear parameter A is calculated in this case as 0.0199. 

As in the first two cases, the trends observed in this case seem to be sumlar, that while 
the identification of the t»„/3 peak is sharper in H-iia},a>,(D) if the excitation force 

level is low, the estimate of the nonlinear parameter, , shows better uniformity over 
the frequency range for a higher excitation force level. Similar observations were made in 
Chapter 3, during the computer simulation of the procedure. The computer simulation 
was carried out for three values of the nondimensional parameter, A (1.0, 0.1 and 0.01), 
and the accuracy of the estimates was seen to be the best in the case oi A = 1.0. The 
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Fig. 6. 1 9(a) Typical sample of the input force: Case3 
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Fig. 6. 1 9(b) Typical sample of the response:Case3 
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Fig. 6.20(a) Power spectrum of the input force: Case 3 
(averaged over 2000 samples). 



Fig. 6.20(b) Power spectrum of the response: Case 3 
(averaged over 2000 samples). 
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values of the nondimensional nonlinearity parameter, in the three cases during 
experimental investigations have been 0.366, 0.114 and 0.0199 and the mean values of 

the estimates of the nonlinear stiffness in the three cases have been 2.99 x 
10^^ N / ,7.16 X 10^^ N / and 19.0 x 10^® N/m^ , respectively. The order of the 

nonlinearity of can be concluded to be as 10^^ N/m^ from the three results, while 
on the basis of the computer simulation results of Chapter 3 and the uniformity of the 
estimates over a frequency range in the experimental investigations, the correct 
magnitude appears to be one in the range suggested by the first two values ( 1 . 49 , 3.58) 
rather than the third (9.50). The major source of error, apart from experimental noise, as 
mentioned in Chapter 3, is that the estimation procedure is restricted to kernels up to the 
third order only, in order to keep the computations to a manageable level. Inclusion of 
higher order kernels (5'*’., order) in the response representation can be expected to 
improve the accuracy of the estimates at increased computational effort. Another source 
of inaccuracy in the estimates is the finite length of samples and the ensemble size. 

6.5 Validation 

An independent check is carried out to validate the estimates of the nonlinear stififiiess. 
This check is based on the analytical formulations of Harris (1984) and Ragulskis et al. 
(1974), which treat the bearing in isolation of the shaft and employ the Hertzian contact 
theory. 

Referring to Fig.6.24, the total elastic force at the points of contact of the ith ball with the 
inner and outer races is expressed as 

F; =A:„(g + JC cos ;?,+>’ sin 77 ,)^'^ (6-2) 

and its projection along the line of action of the applied force is 

F, =K„ig + x COST], + y sin Tj,y^^ COST], (6.3) 

where g is the radial preload or preclearance between the ball and the races and x and y 
are the displacements of the moving ring in the direction of the radial load and 
perpendicular to the direction of the radial load respectively, t], is the angle between the 




Figure 6.24 Schematic diagram of a loaded bearing 


lines of action of the radial load (direction of displacement of the moving ring) and the 
radius passing through the center of the rth ball. K„ is a coefficient of proportionality 
depending on the geometric and material properties of the bearing. The value of K„ , for 

the test bearing is estimated by the method suggested by Harris (1984) as 2.82 x 10^ 
N/mm'". 

The total elastic force in the direction of the applied force is 

^ = ( 6 . 4 ) 

/-I 

where n is the total number of balls in the bearing. 

Using the condition of zero elastic force in the direction perpendicular to the elastic load, 
the deformation, y, perpendicular to the radial force line is expressed as 

= Z [g + Jfcos(77, sin (f7,)/^[g + xcos ( t }, )]''^ sin^ (77, ) ( 6 . 5 ) 

;•=! i»l 

Equations (6.3) and (6.5) are used in equation (6.4) and the bearing stiffness is 
determined as a function of the deformation x as 

k{x) = ^/3c ( 6 . 6 ) 

It can be seen that the bearing stiffiiess is critically dependent on the preloading, g, of the 
balls. While the manufacturer, may, at times, provide the preload range, the exact value 
of the preloading of the bearing balls in the shaft-casing assembly, especially during 
operations which have involved wear and tear, would be difficult to determine. The 
stiffness of the test bearing is plotted in Fig. 6.25 as a function of the radial deformation, 
X, for various allowable preload values, g. 

The same experimental rotor-bearing set-up was employed by Tiwari (1995), for linear 
and nonlinear parameter estimation. He adopted a Markov process approach and modeled 
the governing equations of the system into a Fokker-Planck equation. No external 
excitation was given to the system and the random disturbances during the rotation of the 
shaft, caused due to imperfections and deterioration of the rolling surfaces as well as 




Fig. 6.25 Stiffiiess comparision 
7 - Present study 

1-5 Harris (1984) and Regulskis et al (1974) with preload 0.2, 0.3, 0.4, 0.5 
and 0.6 jj. m respectively 
6 - Tiwari (1995). 




from other random sources, like inaccuracies in alignment etc., were approximated in an 
engineering sense as white and Gaussian. An exact solution of the Fokker-Planck 
equation was obtained, in terms of the probability density functions of the response, a 
curve fitting algorithm was developed to process the solution for parameter estimation. 
The results obtained by Tiwari (1995) along with those from the present procedure are 
also plotted in Fig. 6.25. This comparison is also listed in Table 6.1. The expressions for 
the theoretical stiffiiess in Table 6.1 have been obtained by curve fitting the stiffiiess 
values obtained from equation (6.6), through a quadratic in x. 

Table 6.1 Comparison of bearing stiffness parameters 


j 

Hertzian Contact Theory 
(isolated bearing) 

Tiwari 

(1995) 

Present Study 


Preload (microns) 


Case 1 

Case 2 

Case 3 


0.2 

0.3 

0.4 

0.5 

0.6 





k (KtV/m) 

12.0 

14.7 

16.9 

18.9 

20.8 

22.3 

10.0 

10.0 

10.3 

i' (ESIm')' 

-40.1 

-21.8 

-14.2 

-10.2 

-6.1 

-85.0 

-149.0 

-358.0 

-950.0 


• EN ! rr? = exa-Newton! tr? = Newton In? 


6.6 Remarks 

A good resemblance can be observed between the results obtained in the present study 
' and those obtained by Tiwari and the theoretically possible values. While there is 
variation in the exact magnitude of the estimates, the order of the magmtudes is 
practically identical. This serves to illustrate the correctness of the experimental exercise 
and the practical application aspects of the estimation procedures. 



CHAPTER 7 


CONCLUSION 

The present study has concerned itself with the inverse problem of parameter estimation 
in nonlinear rotor bearing systems. Volterra and Wiener theories of nonlinear analysis 
have been employed as the theoretical platforms for analytical development of 
procedures. A frequency domain approach has been adopted and nonlinear response of 
the system is expressed through first and higher order Volterra and Wiener kernel 
transforms. The procedure involves extraction of Wiener kernels from measurements of 
the applied white noise (broadband in practice) excitation and the resultant response of 
the systems. The Wiener kernels are transformed to Volterra kernels, which are then 
processed for parameter estimation. The following rotor systems were considered 

(i) Rigid Rotor in bearings without cross-coupling 

(ii) Rigid Rotor in bearings with cross-coupling 

(iii) Flexible Rotor 

The procedure was developed in steps - the first rotor configuration mentioned above, 
was treated as a single-degree-freedom-system; the next two cases involved more than 
one degrees of freedom, where the cross-kernel concept was developed through Laplace 
transforms. The procedures have been extensively illustrated through numerical 
simulation. Reasonably good estimates have been obtained for both — linear and 
nonlinear stiffness parameters. 

One major approximation in the present study was the restriction of response 
representation, to include kernels upto the third order only. This has been done in order to 
keep the algebra at a manageable level. The accuracy of estimates is related to the 
accuracy of res}X)nse representation through kernels upto third order. The check, in the 
present study is carried out by performing numerical simulation for various values of the 
nondimensional nonlinear parameters. This nondimensional nonlinear parameter includes 
the linear stiffness and force magnitude terms, in addition to the nonlinear stiffness terms. 
The nonlinear parameters were found to be estimated within the correct order of the 
magnitudes of their exact values, though the accuracy of the values vary with the values 



of the nondimensional parameters chosen for computer simulation. The present study can 
be usefully employed to design experiments to choose appropriate excitation force levels 
for an expected set of values of the stiffiiess parameters. However, it would be 
worthwhile, as a future exercise, to the extend the procedure to include kernels of higher 
order. This can be done for a single-degree-freedom-system, without a great deal of 
additional algebra. However, for systems, with more than one degree of freedom, it 
would help to first develop some generic forms for response and kernel expression. 

Despite the fact that the expressions look formidable with Volterra and Wiener series 
(this would be the case, anyway, for a nonlinear analysis) the computer implementation 
of the algorithms is rather simple and quick. The ensemble size of the force and response 
samples, in the present study, is kept at 2000, due to data storage limitations of the 
computer (A large amount of data was to be stored and retained in the computer, during 
the entire course of the study). The accuracy of the estimated first order kernel transforms 
was checked and the parameter estimation was restricted to a frequency range, upto 
which the errors were low. An attempt can be made, as a future exercise to work with a 
larger ensemble size to improve the accuracy and thereby also increase the usable 
frequency range for parameter estimation. 

The analytical development, in the present study, has been carried out for cubic form of 
nonlinearities. However, the procedures can deal with any general polynomial form of 
nonlinearity. Damping was treated to be linear during this entire exercise. Damping 
nonlinearity can be incorporated in the analysis, though with increased algebraic effort. 
The accuracy of the nonlinear damping estimates, with the present form of response 
representation and the sample and ensemble specifications, can be discussed and 
reviewed only after such a study is made. 

The present study also attempted to implement the algorithms developed on a laboratory 
rotor-rig. The study was carried out with the existing facilities and constraints. The 
results were satisfactory and further strengthened the applicability of the procedures. 
However, the experimental study was restricted to ball bearings. Further experimental 
work needs to cany out on rotors in fluid film bearings, in order to fully investigate the 
utility of the procedures developed during the course of this study. 



REFERENCES 


Amorocho, J., and Orlob, G., 1961, ‘‘"Nonlinear Analysis of Hydrologic Systems", Water 
Resour. Cent. Contrib. 40, Univ. of Calif, Los Ang. 

Ariaratnam, S. T., 1960, "‘‘Random Vibration of Non-linear Suspensions", Journal 
Mechanical Engineering Science 2(3), pp. 195-201. 

Arumugam, P., Swamamani, S., and Prabhu, B. S., 1997, ""An Experimental Investigation 
on Static and Dynamic Characteristics of Journal Bearings Under the Influence 
of Twisting Misalignment", Transactions of ASME, Journal of Tribology, 119, pp. 
188-192. 

Baheti, S. K., and Kirk, R. G., 1994, ""Thermo-Hydrodynamic Solution of Floating Ring 
Seals for High Pressure Compressors using the Finite Element Method', STLE 
Tribology Transactions, 37(2), pp. 336-346. 

Baheti, S. K., and Kirk, R. G., 1994a, ""Finite-Element Thermo-Hydrodynamic Analysis of 
a Circumferencially Grooved Floating Oil Ring SeaF, STLE 94-AM-8D-2. 

Barret, J. F., 1963, “The Use of Functionals in the Analysis of Nonlinear Physical 
Systems", J. Electronic Control, 15, pp. 567-615. 

Bedrosian, E., and Rice, S. O., 1971, ""The Output Properties of Volterra Systems 
(Nonlinear System With Memory) Driven by Harmonic and Gaussian Inputs , 
Proceedings of the IEEE, 59(12), pp. 1688-1707. 

Bedrosian, E., and Rice, S., 1975, ""Application of Volterra-System Analysis. Proceedings 
of 1st Symposium on Test and Identification of Nonlinear Systems , pp. 15-27. 

Bendat, J. S., 1983, ""Statistical Errors for Nonlinear System Measurements Involving 
Square-Law operations". Journal of Sound and Vibration, 90(2), p. 275. 

Bendat, J. S., 1990, ""Nonlinear System Analysis and Identification from Random Data", 
New York: John Wiley. 

Bendat, J. S., 1998, ""Nonlinear System Techniques and Applications", New' York; John 
Wiley. 

Bendat, J. S., and Piersol, A. G., 1982, ""Spectral Analysis of Nonlinear Systems Involving 
Square Law Operations", Journal of Sound and Vibration, 81(2), pp. 199-213. 



Bendat, J. S., and Piersol, A. G., 1986, Random Data: Analysis and Measurement 
Procedures”, 2"^ Edition, Wiley-Interscience, New-York. 

Bendat, J. S., and Piersol, A. G., 1993, ^''Engineering Applications of Correlation and 
Spectral Analysis”, New York: John Wiley. 

Billings, S. A., 1980, '''Identification of Nonlinear Systems-A Survey”, Proceedings of 
lEE, 127(6), pp. 272-285. 

Billings, S. A., and Fakhouri, S. Y., 1978, "Theory of Separable Processes With 
Applications to the Identification of Nonlinear Systems”, Proceedings of lEE, 
125(9), pp. 1051-1058. 

Blok, H., 1 965, "Toplogical Aspects and the Impulse/Whirl Angle Method in the Orbital 
Hydrodynamics of Dynamically Loaded Journal Bearing”, Lecture Notes 
(Condensed English Version), Delft, The Netherlands. 

Bolotin, V. V., 1979, "Random Vibration in Elastic Systems”, Nauka, Moscow (in 
Russian). 

Booker, J., 1965, "Dynamically Loaded Journal Bearings: Mobility Methods of 
Solution”, Journal of Basic Engineering, pp. 537-546. 

Boyd, S., Tang, Y. S., and Chua, L. O., 1983, "Measuring Volterra Kernels”, IEEE 
Transactions on Circuits and Systems, CAS-30(8), pp. 571-577. 

Brilliant, M. B., 1958, "Theory of the Analysis of the Nonlinear systems”, MIT Res. Lab. 
of Elect., Technical report, 304. 

Brillinger, D. R., 1970, "The Identification of Polynomial Systems by means of Higher 
order Spectra”, Journal of Sound and Vibration 12(3), pp. 301-313. 

Burrows, C. R, Sayed-Esfahani, R., and Stanway, R., 1981, "A Comparison of 
Multifrequency Techniques for Measuring the Dynamics of Squeeze Film 
Bearings”, Transactions of ASME, Journal of Lubrication Technology, 103, pp. 
137-143. 

Burrows, C. R., and Sahinkaya, M. N., 1982, "Frequency Domain Estimation of 
Linearised Oil Film Coefficients”, Transactions of ASME, Journal of Lubrication 
Technology, 104, 210. 



Cai, G. Q. and Lin, Y. K., 1988, “A New Approximate Solution Technique for Randomly 
Excited Non-Linear Oscillators”, International Journal of Non-linear Mechanics, 
23, pp. 409-420. 

Caughey, T. K., 1963, Derivation and application of the Fokker-Planck Equation to 
Discrete Non-Linear Dynamic Systems Subjected to White Random Excitation”, 
Journal of Acoustical Society of America, 35(1 1), pp. 1683-1692. 

Caughey, T. K., 1971, ‘‘‘'Nonlinear Theory of Random Vibrations”, Advances in Applied 
Mechanics, 11, pp. 209-253, Academic Press, New York. 

Caughey, T. K., 1986, the Response of Non-Linear Oscillators to Stochastic 
Excitation”, Probabilistic Engineering Mechanics, 1, pp. 2-4. 

Caughey, T. K. and Ma, F., 1983, “The Exact Steady State Solution of a Class of 
Nonlinear Stochastic Systems”, International Journal of Nonlinear Mechanics, 17, 
pp. 137-142. 

Chen, J. H., and Lee, A. C., 1997, “Identification of Linearised Dynamic Characteristics 
of Rolling Element Bearings”, Transactions of ASME, Journal of Vibration and 
Acoustics, 119, pp. 60-69. 

Childs, D., 1993, “Turbomachinery Rotordynamics: Phenomena, Modeling and 
Analysis”, New York: A Wiley-Interscience Publication. 

Childs, D., Moes, H., and Van Leeuwen, H., 1977, “Journal Bearing Impedance 
Description for Rotor Dynamic Applications”, Transactions of ASME, Journal of 
Lubrication Technology, pp. 198-214. 

Choi, D., Miksad, R. W, and Powers, E. J, 1985, “Application of Digital Cross-Bispectral 
Analysis Techniques to Model the Nonlinear Response of a Moored Vessel System 
in Random Seas”, Journal of Sound and Vibration, 99, pp. 309-326. 

Choy, F. K., Braun, M. J. and Hu Y., 1992, “Nonlinear Transient and Frequency 
Response Analysis of a Hydrodynamic Journal Bearing', Transactions of the 
ASME, Journal of Tribology 114, pp. 448-454. 

Chu., C. S., Wood, K. L., and Busch-Vishniac, I. J., 1998, “A Nonlinear Dynamic Model 
With Confidence Bounds for Hydrodynamic Bearings”, Transactions of ASME , 
Journal of Tribology, 120, pp. 595-604. 



Chua, L. O., and Ng., C. Y., 1979, ^Frequency Domain Analysis of Nonlinear Systems: 

General Theory', lEE Journal on Electronic Circuits and Sys., 3(4), pp. 165-185. 
Chua, L. O., and Ng., C. Y., 1979a, ''Frequency Domain Analysis of Nonlinear Systems: 
Formulation of Transfer Functions”, lEE Journal on Electronic Circuits and 
Systems, 3(6), pp. 257-267. 

Crandall, S. H., 1963, "Perturbation Techniques for Random Vibrations of Nonlinear 
Systems”, Journal of Acoustical Society of America, 35(14), pp. 1700-1705. 
Crandall, S. H., 1966, "Applied Mechanics Surveys”, Random vibration, Abramson, H. 
N., Liebowitz, H., Crowley, J. M., and Juhasz, S., eds., Sparton Books, 
Washington, D. C., pp. 681-689. 

Crandall, S. H., 1985, "Non-Gaussian Closure Techniques for Stationary Random 
Vibration”, International Journal of Non-linear Mechanics, 20, pp. 1-8. 

Crum, L., 1975, "Simultaneous Reduction and Expansion of Multidimensional Laplace 
Transform Kernels”, Proceedings of 1st Symposium on Test and Identification of 
Nonlinear Systems, pp. 89-105. 

DeFigueiredo, R., and Dwyer, T., 1980, "A Best Approximation Framework and 
Implementation for Simulation of Large Scale Nonlinear fystems”, IEEE 
Transactions on Circuits and Systems, CAS-27, pp. 1005-1014. 

Dimentberg, M. F., 1980, Nonlinear Stochastic Problems of Mechanical Vibrations, 
Nauka, Mascow (in Rassian). 

Dimentberg, M. F., 1988, "Statistical Dynamics of Nonlinear and Time-Varying 
Systems”, Research studies press limited, John Wiley and Sons Inc., New York. 
Dimentberg, M. F. and Sokolov, A. A., 1991, "Identification of Restoring Force Non- 
Linearity from a System 's Response to a White-Noise Excitation”, International 
Journal of Nonlinear Mechanics, 26(6), pp. 851-855. 

Elsermans, M., Hongerloot, M. and Snokeys, R., 1975, "Damping in Taper Rolling 
Bearings”, Proc. 16th MTDR Conference, p. 223, Manchester, 10-12 September. 
Eschmann, P., Hasbargen, I. and Weigand, K., 1985, "Ball and Roller Bearings. Theory. 

Design and Application” , John Wiley and Sons; New York. 

Ewins, D. J., 1984 "Modal Testing :Theory and Practice”, Lectchworth, England: 
Research Studies Press. 



Eykhoff, P., 1974, "‘System Identification: Parameter and State Estimation", Wiley, New 
York, pp. 94-104. 

Eykhoff, P., Grinten, P. M. Van Der, Naak, H. K., and Veltman, B. P., 1966, "System 
Modeling and Identification", Proc. Illrd IFAC Congress, London. 

Fakhouri, S. Y., 1980, "Identification of a Class of Nonlinear Systems With Gaussian 
Non-White Inputs", International Journal of Systems Science, 11(5), pp.-541-555. 

Fakhouri, S. Y., and Billings, S. A., 1980, "Identification of Nonlinear Systems Using 
Correlation Analysis and Pseudorandom Inputs", International Journal of 
Systems Science, 11(3), pp.-26 1-279. 

Fakhouri, S. Y., Billings, S. A., and Wormald, C. N., 1981, "Analysis of Estimation 
Errors in the Identification of Nonlinear Systems", International Journal of 
Systems Science, 12(2), pp. 205-225. 

Falsone G., 1992, "Stochastic Linearisation of MDOF Systems Under Parametric 
Excitations", International Journal ofNonlinear Mechanics 27(6), pp. 1025-1037. 

Flake, R. H., 1963, "Volterra Series Representation of Nonlinear Systems", IEEE 
Transactions on App. Ind., 81, pp. 330-335. 

Fokker, A. P., 1914, "Die mittlere energie rotierender electricher dipole im 
stahlungsfeld'. Annals Phys., 43, pp. 810-815. 

Fredriksson, B., 1980, "Three-Dimentional Roller-Raceway Contact Stress Analysis", 
Advanced Engineering Corp. Report, Sweden. 

French, A. S., and Butz E. G., 1973, "Measuring the Wiener Kernels of a Nonlinear 
System using the Fast Fourier Transform Algorithm", International Journal of 
Control, 17(3), pp. 529-539. 

French, A. S., and Butz, E. G., 1974, "The use of Walsh Functions in the Wiener Analysis 
of Nonlinear Systems", IEEE Transactions on Computers, C-23 (3), pp. 225-232. 

Fretchet, M., 1910, "Sur les fonctionnels continues", Ann. Ec. Norm. Sup., 27, pp. 193- 
219. 

Fuller, A. T., 1969, "Analysis ofNonlinear Stochastic Systems by means of the Fokker- 
Planck Equation", International Journal of Control, 9 (6), pp. 603-655. 

Fuller, D. D, 1984, "Theory and Practice of Lubrication for Engineers", second edition. 
New York; A Wiley-Interscience Publication. 



George, D. 1959, '"Continuous Nonlinear Systems'', Res. Lab. Electron., M.I.T., 
Cambridge, Mass., Tech. Rep. 355. 

Gifford, S. J., and Tomlinson, G. R., 1989, "Recent Advances in the Application of 
Functional Series to Nonlinear Structures", Journal of Sound and Vibration, 
135(2), pp. 289-317. 

Glienicke, J., 1966, "Experimental Investigation of the Stiffness and Damping 
Coefficients of Turbine Bearings and their Application to Instability Prediction", 
Proceedings of Institute of Mechanical Engineers 181 (3b), p. 116. 

Goodwin, M. J., 1991, "Experimental Techniques for Bearing Impedance Measurement, 
Transaction of ASME, Journal of Engineering for Industry, 113, pp. 335-342. 

Goodwin, M. J., Ogrodnik, P. J., Roach, M. P., and Fang, Y., 1997, "Calculation and 
Measurement of the Stiffness and Damping Coefficients for a Low Impedance 
Hydrodynamic Bearing", Trans, of ASME, Journal of Tribology, 119, pp. 57-63. 

Gross, W., 1962, "Gas Film Lubrication", New York: Wiley. 

Gustavsson, O. and Tallian, T., 1962, "Detection of Damage in Assembled Rolling 
Element Bearings", Trans, of American Society of Locomotive Engineers 5. 

Hagg, A. C., 1946, "The Influence of Oil Film Journal Bearings on the Stability of 
Rotating Machines", Transactions of ASME, Journal of Applied Mechanics, 13, 
pp. A21 1-A220. 

Halme, A., Orava, J., and Blomberg, H., 1971, "Polynomial Operators in Nonlinear 
Systems Theory", International Journal of System Sciences, 2(1), pp. 25-47. 

Harris, G. H., and Lapidus, L., 1967, "The Identification of Nonlinear Systems", 
Industrial and Engineering Chemistry, 59(6), pp. 66-81. 

Harris, T. A., 1984, "Rolling Bearing Analysis", New York: Wiley. 

Hertz, H., 1896, "On the Contact of Rigid Elastic Solids and on Hardness", 
Miscellaneous Papers, Macmillan, London, pp. 163-183. 

Hess, D. P., Soom, A. and Kim, C. H., 1992, "Normal Vibrations and Friction at a 
Hertzian Contact Under Random Excitation: Theory and Experiments", Journal 
of Sound and Vibration 153 (3), pp. 491-508. 



Ho, D., and Stark, L, 1973, '"Identification of Biological Systems by Volterra Series”, 
Proceedings of 4*’’ Symposium on Nonlinear Estimation Theoiy and its 
Application. 

Honrath, K., 1960, "Concerning the Stiffness of Machine Tool Spindles and Their 
Bearings”, Thesis submitted for Dicktor Ingenieur at the Technical University of 
Aachen. 

Hummel, C., 1926, "Kristische Drehzahlen Als Folge Der Nachgiebigkeit Des 
Schmiermittels Im Larger”, VDI-Forschungsheft, 287. 

Hung, G., and Stark, L., 1977, "Introductory Review: The Kernel Identification Method: 
Review of Theory, Calculation, Application, and Interpretation,” Mathematical 
Biosciences, 37, pp. 135-190. 

Hung, G., Stark, L., and Eykhoff, P., 1977, "On the Interpretation of Kernels. I 
Computer Simulation of Responses to Impulse Pairs”, Annals of Biomedical 
Engineering, 5, pp. 130-143. 

Ibrahim, R. A. and Roberts, J. W., 1978, "Parametric Vibrations, Part 5: Stochastic 
Problems”, Shock and Vibration Digest, 10, pp. 17-38. 

Jacoby, S., 1966, "A Mathematical Model for Nonlinear Hydrologic Systems”, J. 
Geophys. Rev., 71, pp. 4811-4824. 

Kalker, K., 1982, "Numerical Calculation of the Elastic Field in a Half-Space due to an 
Arbitrary Load Distributed Over a Bounded Region of the Surface”, SKF 
Engineering and Research Center report NL82D002 (Appendix). 

Katzenelson J., Gould, L., 1962, "The Design of Nonlinear Filters and Control Systems”, 
Part I, Inf. Control, 5, pp. 108-143. 

Khasminskii, R. Z., 1966, "A Limit Theorem for the Solutions of Differential Equations 
With Random Right-Hand Sides”, Theory of Probability and Applications, 11, pp. 
390-405. 

Kim, K. I., and Powers, E. J., 1988, "A Digital Method for Modeling Quadratically 
Nonlinear Systems with a General Random InpuC, IEEE Transactions on 
Acoustics, Speech and Signal Processing, ASSP-36, pp. 1758-1769. 

Kirk, R. G., and Gunter E. J., 1970, "Transient Journal Bearing Analysis”, NASA CR 


1599. 



Kirk, R. G., and Gunter E. J., 1975, “Short Bearing Analysis Applied to Rotor Dynamics, 
Part I: Theory", ASLE-ASME Joint Lubrication Conference 75-Lub-30. 

Kirk, R. G., and Gunter E. J., 1975a, “Short Bearing Analysis Applied to Rotor 
Dynamics, Part II: Results of Journal Bearing Response", ASLE-ASME Joint 
Lubrication Conference 75-Lub-31. 

Klein, S., and Yasui, S., 1976, “Nonlinear Systems Analysis with General Non-Gaussian 
Stimuir, Tech. Rep. Div. Of Biol., Calif, Inst. OfTechnol. 

Koh, T., and Powers, E. J., 1985, “Second Order Volterra Filtering and its Application to 
Nonlinear System Identification”, IEEE Transactions on Acoustics, Speech and 
Signal Processing, ASSP-33, pp. 1445-1455. 

Kolmogorov A., 1931, “Uber Die Analytischen Methoden in 

Wahrsheinlicheitsrechnung", Mathematische Annalen, 104, pp. 415-458. 

Kononenko, V. 0. and Plakhtienko, N. P., 1970, “Determination of Nonlinear Vibration 
System Characteristics from Analysis of Vibrations", Prikladnaya Mekhanika, 6, 
Vyp. 9. 

Korenberg, M. J., 1973, “Identification of Nonlinear Differential Systems", JACC, pp. 
597-603. 

Korenberg, M. J., and Hunter, 1. W., 1990, “The Identification of Nonlinear Biological 
Systems: Wiener Kernel Approaches" , Annals of Biomedical Engineering, 18, pp. 
629-654. 

Kraus, J., Blech, J. J. and Braun, S. G., 1987, “In Situ Determination of Rolling Bearing 
Stiffness and Damping by Modal Analysis", Transaction of ASME, Journal of 
Vibration, Acoustics, Stress, and Reliability in Design, 109, pp. 235-240. 

Krausz, H. I., 1975, “Identification of Nonlinear Systems using Random Impulse Train 
Input", Biological Cybernetics 19, pp. 217-230. 

Krodkiewski, J. M., and Ding, J., 1993, “Theory and Experiment on a Method for On-Site 
Identification of Configurations of Multi-Bearing Rotor Systems", Journal of 
Sound and Vibration, 164(2), pp. 281-293. 

Krodkiewski, J. M., and Ding, J., 1993a, “Inclusion of Static Indetermination in the 
Mathematical Model for Nonlinear Dynamic Analyses of Multi-Bearing Rotor 
System", Journal of Sound and Vibration, 164(2), pp. 267-280. 



Kunert K., 1961, '’‘‘Spannungsverteilung im halbraum bei elliptischer 

flchenpressungsverteilung ber einer rechteckigen druckfi„che’\ Forsch. Geb. 
Ingerieurwes 27 (6), pp. 165-174. 

Lee, Y. W., Levinson, N., and Martin, W.T., 1964, ^'‘Selected Papers ofNorbert Wiener", 
Published Jointly by Society for Industrial and Applied Mathematics and the M. I. 
T. Press, Cambridge, Massachusetts. 

Lee, Y. W., and Schetzen, M., 1965, "Measurement of Wiener Kernels of Nonlinear 
System by Crosscorrelation", International Journal of Control, 2(3), pp. 237-254. 
Leon, B. J., and Schaefer, D. J., 1978, "Volterra Series and Picard Iteration for 
Nonlinear Circuits and Systems," IEEE Transactions on Circuits and Systems, 
CAS-25,pp. 789-793. 

Lesiak, C., and Krener, A. J., 1978, “Existence and Uniqueness of Volterra Series," IEEE 
Transactions on Automatic Control, AC-23, pp. 1090-1095. 

Levy, E. C., 1959, “Complex-Curve Fitting ” IRE Transactions on Automatic Control, 
AC-4, pp. 37-44. 

Lim, T. C. and Singh, R., 1990a, “Vibration Transmission Through Rolling Element 
Bearings, Part I: Bearing Stiffness Formulation" , Journal of Sound and Vibration 
139 (2), pp. 179-199. 

Lim, T. C. and Singh, R., 1990b, “Vibration Transmission Through Rolling Element 
Bearings, Part 11: System Studies", Journal of Sound and Vibration 139 (2), pp. 
201-225. 

Lim, T. C. and Singh, R., 1991, “Vibration Transmission Through Rolling Element 
Bearings, Part III: Geared Rotor System Studies", Journal of Sound and Vibration 
151 (l),pp. 31-54. 

Lim, T. C. and Singh, R., 1992, “Vibration Transmission Through Rolling Element 
Bearings, Part IV: Statistical Energy- Analysis", Journal of Soimd and Vibration 

153(l),pp. 37-50. 

Lim, T. C. and Singh, R., 1994, “Vibration Transmission Through Rolling Element 
Bearings, Part V: Effect of Distributed Contact Load on Roller Bearing Stiffness 
Matrix", Journal of Sound and Vibration 169 (4), pp. 547-553. 



Lin, Y. K. and Cai, G. Q., 1995, Probabilistic Structural Dynamics. Advanced Theory 
and Applications. New York: McGraw-Hill. 

Lehman, G., 1953, ^^Untersuchung des Laufgerausches von Walzlgern'\ Konstruktion 5. 

Lubbock, J. K. and Bansal V. S., 1969, ‘‘'’Multidimensional Laplace Transforms for 
Solution of Nonlinear Equations'", Proceedings of lEE, 116(12), pp. 2075-2082. 

Lund, J. W., 1968, ‘‘‘’Calculation of Stiffness and Damping Properties of a Gas Bearings," 
Transactions of ASME, Journal of Lubrication Technology, pp. 793-803. 

Lund, J. W., 1987, ‘‘‘’Review of the Concept of Dynamic Coefficients for Fluid Film 
Bearings," Transactions of ASME, Journal of Tribology, 109, pp. 37-41. 

Lund, J. W., and Stemlicht, B., 1962, “Rotor Bearing Dynamics with Emphasis on 
Attenuation", Transactions of ASME, Journal of Basic Engineering, pp. 491-502. 

Marmarelis, P., and Naka, K. I., 1973, “Nonlinear Analysis and Synthesis of Receptive 
Field Responses in the Cat-Fish Retina-1: Horizontal Cell Ganglion Cell Chain", 

J. NeurophysioL, 36, pp. 605-618. 

Marmarelis, P., and Naka, K. I., 1973b, “Nonlinear Analysis and Synthesis of Receptive 
Field Responses in the Cat-Fish Retina-ll: One Input White Noise Analysis", J. 
Neurophysiol, 36, pp. 619-633. 

Marmarelis, P., and Naka, K. I., 1973c, “Nonlinear Analysis and Synthesis of Receptive 
Field Responses in the Cat-Fish Retina-III: Two Input White-Noise Analysis", J. 
Neurophysiol, 36, pp. 634-648. 

McFadden, P. D. and Smith, J. D., 1984, “Vibration Monitoring of Rolling Element 
Bearings by the High-Frequency Resonance Technique - A Review", International 
Journal of Tribology, 17, pp. 3-10. 

McFadden, P. D. and Smith, J. D., 1985, “The Vibration Produced by Multiple Point 
Defects in a Rolling Element Bearing", Journal of Sound and Vibration, 98 (2), 
263-273. 

Morton, P. G., 1971, “Measurements of Dynamic Characteristics of a Large Sleeve 
Bearings", Transactions of ASME, Journal of Lubrication Technology, 143. 



Muszynska, A. and Bently, D. E., 1990, ""Frequency-Swept Rotating Input Perturbation 
Techniques and Identification of the Fluid Force Models in Rotor/Bearing/Seal 
Systems and Fluid Handling Machines'", Journal of Sound and Vibration, 143 (1), 
pp. 103-124. 

Myrick, T. Jr., and Rylander H. G., 1975, ""Analysis of Flexible Rotor", ASME Design 
Engineering Technology Conference, Washington D.C. 75-Det-39. 

Nam, S. W., Kim, S.B., and Powers, E. J., 1990, ""Nonlinear System Identification With 
Random Excitation Using Discrete Third Order Volterra Series", Conference 
Proceedings of 8*’’ International Modal Analysis, pp. 1278-1283. 

Newkirk, B. L., 1925, ""Shaft Whippin^\ General Electric Review, 169. 

Newkirk, B. L., and Taylor, H. D., 1924, ""Shaft Whipping Due to Oil Action in Journal 
Bearings", General Electric Review, pp. 559-568. 

Nikias, C. L., and Raghuveer, M. R., 1987, ""Bispectral Estimation: Digital Signal 
Processing FrameworF', Proceedings of IEEE, 75, pp. 869-891. 

Nordman, R., and Schollhom, K., 1980, ""Identification of Stiffness and Damping 
Coefficients of Journal Bearings by Means of Impact Method’', Institution of 
Mechanical Engineers Conference, Vibration of Rotating Machinery, Cambridge. 

Ocvirk, F., 1952, ""Short Bearing Approximations for Full Journal Bearings", NASA 
Tech. Note 2808. 

Odiari, E. A., and Ewins, D. J., 1992, ""Parameter Identification for Nonlinear Rotor- 
Stator Systems- the Volterr a/Wiener Based Approach", Conference Proceedings 
of Institution of Mechanical Engineers (IMechE), Vibrations in Rotating 
Machineries, pp. 193-202. 

Orabi, I. I., and Ahmadi, G., 1987, ""A Functional Series Expansion Method For 
Response Analysis of Nonlinear Systems Subjected to Random Excitations", 
International Journal of Nonlinear Mechanics, 22(6), pp. 451-465. 

Orcutt, F. K., Arw-as, E. B., 1967, ""The Steady State and Dynamic Characteristics of Full 
Circular Bearing and a Partial Arc Bearing in Laminar and Turbulent Flow 
Regimes", Transactions of ASME, Journal of Lubrication Technology 89 Series 
r(2),pp. 143-153. 

Palmgren, A., 1959, ""Ball and Roller Bearing Engineering", Third Edition, Burbank. 



Papanicolaou, G. C., and Kohler, W., 1974, ‘‘Asymptotic Theory of Mixing Stochastic 
Ordinary Differential Equations", Comm, on Pure and Applied Mathematics, 27, 
pp. 614-668. 

Planck, M., 1917, ‘‘‘Uber einen satz der statischen tfynamik und seine erweiterung in der 
quanten-theorie", Sitzungber, Preuss. Akadamie Weiss, pp. 324-341. 

Powers, E. J., and Miksad, R. W., 1987, “Polyspectral Measurement and Analysis of 
Nonlinear Wave Interactions in Fluids", ASME Winter Annual Meeting, Boston. 

Pradlwarter, H. J., Bucher, C. G. and Schueller G. I., 1991, “Structural Dynamics, Recent 
Advances", Editor Schueller, G. /., New York: Springer - Verlag. 

Raemer, H. R., 1969, “Statistical Communication Theory and Applications", Englewood 
Cliffs, New Jersey: Prentice-Hall. 

Ragulskis, K. M., Jurkauskas, A. Yu., Atstupenas, V. V., Vitkute, A. Yu., and Kulvec, A. 
P., 1974, “Vibration of Bearings", Vilnyus: Mintis Publishers. 

Rao, J. S., 1998, “Rotor Dynamics", Third Edition, New Age International (P) Limited 
Publishers. 

Reusner, H., 1977, “Drucl<fl„chenbelastung und Overfluchenverschiebung in 

W„lzkontakt von Rot.Jionskurpern", Dissertation, Schweinfurt, West Germany. 

Rice, H. J., and Fitzpatrick, J. A., 1988, “A Generalised Technique For Spectral Analysis 
of Nonlinear Systems", Mechanical Systems and Signal Processing 2(2), pp. 195- 
207. 

Rice, H. J., and Fitzpatrick, J. A., 1991, “A Procedure for the Identification of Linear and 
Nonlinear Multi-Degree of Freedom Systems", Journal of Sound and Vibration, 
149(3), pp. 397-411. 

Roberts, J. B., 1981, “Response of Nonlinear Mechanical Systems to Random Excitation, 
Part 1: Markov Method", Shock and Vibration Digest, 13 (4), pp. 17-28. 

Roberts, J. B., 1983, “Energy Method for Nonlinear Systems With Non-White Excitation , 
Proceedings lUTAM Symposium Random Vibrations and Reliability, K. Hennig, 
Ed, Akademie, Berlin, pp. 285-294. 

Robertson, D., 1946, “Whirling of a Journal in a Sleeve Bearing", Philosphical 
Magazine, 15(7), pp. 113-130. 



Rouvas C. and Childs D. W., 1993. Parameter Identification Method for the 
Rotordynamic Coefficients of a High Reynolds Number Hydrostatic Bearing', 
Transactions ASME, Journal of Vibration and Acoustics, 115, pp. 264-270. 

Roy, R. V., and Spanos, P. D., 1990, “fViener-Hermite Functional Representation of 
Non-Linear Stochastic Systems", Journal of Structural safety. 

Rubinstein, R. Y., 1981, “Simulation and the Monte Carlo Method', Wiley, New York. 

Rugh, W. J., 1981, “Nonlinear System Theory: The Volterra/Wiener Approach" , Johns 
Hopkins University Press, Baltimore, Maryland 

Sahinkaya, M. N., and Burrows, C. R., 1984, “Estimation of Linearised Oil Film 
Parameters from Out of Balance Response", Proceedings of Institute of 
Mechanical Engineers 198c(8), 131. 

Sandberg, I. W., 1982, “Expansions for Nonlinear Systems", Bell System Tech. J., 61 pp. 
159-200. 

Sayler, R., deSilva, G., Leather, J., Anderson, J. and MacPherson, P., 1981, “Elastic 
Conformity in Hertzian Contacts", Tribol. Int., pp. 315-322. 

Schetzen, M., 1965, “Measurement of the Kernels of a Nonlinear System of Finite 
Order," International Journal of Control, 1(3), pp. 251-263. 

Schetzen, M., 1965a, “Synthesis of a Class of Nonlinear Systems," International Journal 
of Control, 1(5), pp. 401-414 

Schetzen, M., 1974, “A Theory of Nonlinear System Identification^ International Journal 
of Control, 20(4), pp. 577-592. 

Schetzen, M., 1980, “The Volterra and Wiener Theories of Nonlinear Systems", New 
York; John Wiley and Sons. 

Shinozuka, M., 1972, “Monte Carlo Solution of Structural Dynamics", Computers 
Structures, 2, pp. 855-874. 

Simandri, S., and Hahn, E. J., 1975, “Effect of Pressurization in the Vibration of Isolation 
Capability of Squeeze Film Bearings", ASME Design Engineering Technology 
Conference, Washington D. C., 75-Lub-70. 

Soong, T. T. and Grigoriu, M., 1993, Random Vibration of Mechanical and Structural 
Systems. New Jersy: Prentice-Hall. 



Spanos, P. D., 1981, ''‘'Monte Carlo Simulations of Responses ofNon-Symmetric Dynamic 
Systems to Random Excitations", Computers Structures, 13, pp. 371-376. 

Spanos, P. D., 1983, "Approximate Analysis of Random Vibration Problems Through 
Stochastic Averaging", Proceedings lUTAM Symposium Random Vibrations and 
Reliability, K. Hennig, Ed, Akademie, Berlin, pp. 327-337. 

Spanos P. D. and Donley M. G., 1991, "Equivalent Statistical Quadratization for Non- 
Linear Systems", Journal of Engineering. Mechanics, 117(6), pp. 1289-1310. 

Spanos P. D. and Donley M. G., 1992, "Non-Linear Multi-Degree-of-Freedom-System 
Random Vibration by Equivalent Statistical Quadratization", International 
Journal of Non-Linear Mechanics, 27(5), pp. 735-748. 

Spanos, P. D. and Mignolet, M. D., 1989, "ARMA Monte Carlo simulation in 
probabilistic structural analysis". Shock Vibration Digest, 21, pp. 3-14. 

Stanway, R., Burrows, C. R., and Holmes, R., 1979, 1. Mechanical Engineering Science, 
Institute of Mechanical Engineers, 21(6). 

Stodola, A., 1925, "Kristiche Wellenstorong Infolge Der Nachgiebigkeit Des Oelpolsters 
Im Larger", Schweizerische Bauzeiting, 85, pp. 265-266. 

Stoker, J. J., 1950, "Non-linear Vibrations", Interscience Publishers. 

Stolarski, T. A., 1990, "Tribology in Machine Design", Oxford: Heinemann Newnes. 

Stratonovitch, R. L., 1967, "Topics in the Theory of Random Noise", 2, Gordon and 
Breach, New York. 

Su, Y. T., Lin, M. H., and Lee M. S., 1993, "The effects of Surface Irregularities on 
Roller Bearing Vibrations", Journal of Sound and Vibration, 165 (3), pp. 455-466. 

Sunnesjo, C. S., 1985, "Rolling Bearing Vibrations - The Effects of Geometrical 
Imperfections and Wear", Journal of Sound and Vibration, 98 (4), pp. 455-474. 

Swanson, E. E., and Kirk, R. G., 1997, "Survey of Experimental Data for Fixed Geometry 
Hydrodynamic Journal Bearings", Transactions of ASME, Journal of Tnbology, 

119, pp. 704-710. 

Tick, L., 1961, "The Estimation of Transfer Functions of Quadratic Systems", 
Technometric 3(4) pp. 563-567. 



Tiwari, R. and Vyas, N. S., 1995, “Estimation of Nonlinear Stiffness Parameters of 
Rolling Element Bearings from Random Response of Rotor Bearing Systems", 
Journal of Sound and Vibration, 187(2,) pp.229-239. 

Tiwari, R. and Vyas, N. S., 1997, “Nonlinear Bearing Stiffness Parameter Extraction 
from Random Response in Flexible Rotor-Bearing Systems", Journal of Sound 
and Vibration, 203(3), pp. 389-408. 

Tiwari, R. and Vyas, N. S., 1997a, “Parameter Estimation in Imbalanced Nonlinear 
Rotor-Bearing Systems from Random Response", Journal of Sound and Vibration, 
208(1), pp. 1-14. 

Tiwari, R. and Vyas, N. S., 1998, “Stiffness Estimation from Random Response in Multi- 
Mass Rotor-Bearing Systems", Probabilistic Engineering, Mechanics, 13(4), pp. 
255-268. 

Tonneson, J., 1975, “Experimental Parametric Study of a Squeeze Film Bearing', ASLE- 
ASME Joint Lubrication Conference, Miami Beach, Fla., 75-Lub-42. 

Vance, J. M., and Kearton, A. J., 1975, “Experimental Measurement of the Dynamic 
Force Response of a Squeeze Film Damper", ASME Design Engineering 
Technology Conference, Washington D.C., 75-DET-39. 

Vandiver. J. K., and Jong, J. Y., 1986, “The Relation Between in Line and Cross Flow 
Vortex Induced Vibration of Cylinders", Proceedings of 1st OMAE Symposium 
on Offshore and Arctic Frontiers, ASME, pp. 299-309. 

Vasta M, 1995, “Exact Stationary Solution for a Class of Non-Linear Systems Driven by 
a Non-Normal Delta-Correlated Process", International Journal of Non-Linear 
Mechanics, 30(4), pp. 407-418. 

Volterra, V., 1889, “Surr une generlisation de la theorie des functions dune variable 
imaginaire. Ler Memoire", Acta Math., Stockh., 12, pp. 233-286. 

Volterra, V. 1930, “Theory of Functionals", Blackie and Sons, Glasgow. 

Volterra. V., 1959, “Theory of Funcliomls <md of Imegral and Integro-Differenlial 

Equalions’’, Dover Publications, Inc., New York. 

Walford. T. L. H. and Stone B. J., 1980, “The Measurement of the Radial Stiffness of 
Rolling Element Bearings Under Oscillation Condition^, Journal of Mechanical 
Engineering Science 22 (4), pp. 175-181. 



Warner. P. C., 1963, “Static and Dynamic Properties of Partial Journal Bearings", 
Transactions of ASME, Journal of Basic Engineering, 85(D), pp. 247-257. 
Watanabe, A., 1975, “The Volterra Series Expansion of Nonlinear Functionals Which are 
Defined on the Product of two Finite- Dimensional Vector Spaces, and its 
Application to Characterisation of Two Input Systems", Summer Comput. 
Simulation Conference. 

Watanabe, A., and Stark, L., 1975, “Kernel Method for Nonlinear Analysis: Identification 
of a Biological Control System", Math. Biosci., 27, pp. 99-108. 

While, M. F., 1979, “Rolling Element Bearing Vibration Transfer Characteristics: Effect 
of Stiffness", Journal of Applied Mechanics, Trans, of ASME, 46, pp. 677-684. 
Wiener, N., 1958, “Nonlinear Problems in Random Theory", The Massachusetts Institute 
of Technology and John Wiley and Sons ,Inc., New York. 

Worden, K., Manson, G. and Tomlinson, G. R., 1997, “A Harmonic Probing Algorithm 
for the Multi-Input Volterra Series", Journal of Sound and Vibration, 201(1), pp. 
67-84. 

Yamamouchi, Y., 1974, “Ship ’s Behaviour in Waves as a Stochastic Process", Dynamics 
of Marine Vehicles and Structures in Waves (Ed. R.E.D. Bishop and W.G. Price) 
London: Institution of Mechanical Engineers, 

Zhou R. S. and Hashimoto F. 1995, “A New Rolling Contact Surface and “No Run-In- 
Performance Bearings", Transactions of the ASME, Journal of Tribology, 117, 

pp. 166-170. 

Zhu, W. Q. and Yu, J.S., 1989, “The Equivalent Nonlinear System Methods", Journal of 
Sound and Vibration, 129, pp. 385-395. 

Zhu, W. Q., 1983, “Stochastic Averaging of the Energy Envelope of Nearly Lyapunov 
Systems", Proceedings lUTAM Symposium Random Vibrations and ReUabihty, K. 

Henning, Ed, Akademie, Berlin, pp. 347-357. 

Zw-ilein, 0. and Schlicht, 1980, “Werkstoffanstrengung bei walzbeanspruchung - einfluss 

von reibung und eigenspannungen", Werkstoffiech, 11, pp- 1-14- 



LIST OF PAPERS PUBLISHED/COMMUNICATED 
FROM THE PRESENT STUDY 


Khan, A. A., and Vyas, N. S., 1998, “Nonlinear Parameter Estimation Using Volterra 
and Wiener Theories'”, (in print). Journal of Sound and Vibration. 

Khan, A. A., and Vyas, N. S., 1998, “Parameter Estimation in a Nonlinear Coupled 
System Through Volterra and Wiener Kernels”, commxinicated to Journal of 
Sound and Vibration. 

Khan, A. A., and Vyas, N. S., 1999, “Nonlinear Bearing Stiffness Parameter Estimation 
in Flexible Rotor-Bearing Systems Using Volterra and Wiener Approach”, 
communicated to Probabilistic Engineering Mechanics. 

Khan, A. A., and Vyas, N. S., 1999, “Analytical and Experimental Investigations for 
Nonlinear Rotor Bearing Stiffness Estimation Through Volterra and Wiener 
Kernels”, communicated to Mechanism and Machine Theory. 




